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Chapter 1 
Maxwell's electrodynamics 



1.1 Dynamical variables of electromagnetism 

The classical theory of the electromagnetic field, as formulated by Maxwell, involves 
the vector fields 

E(t, x) = electric field at position x and time t (1.1.1) 

and 

Bit, x) = magnetic field at position x and time t. {1.1.2) 

We have two vectors at each position of space and at each moment of time. The 
dynamical system is therefore much more complicated than in mechanics, in which 
there is a finite number of degrees of freedom. Here the number of degrees of freedom 
is infinite. 

The electric and magnetic fields are produced by charges and currents. In a 
classical theory these are best described in terms of a fluid picture in which the 
charge and current distributions are imagined to be continuous (and not made of 
pointlike charge carriers). Although this is not a true picture of reality, these con- 
tinuous distributions fit naturally within a classical treatment of electrodynamics. 
This will be our point of view here, but we shall see that the formalism is robust 
enough to also allow for a description in terms of point particles. 

An element of charge is a macroscopically small (but microscopically large) 
portion of matter that contains a net charge. An element of charge is located 
at position x and has a volume dV. It moves with a velocity v that depends on 
time and on position. The volume of a charge element must be sufficiently large 
that it contains a macroscopic number of elementary charges, but sufficiently small 
that the density of charge within the volume is uniform to a high degree of accuracy. 

In the mathematical description, an element of charge at position x is idealized as 
the point x itself. The quantities that describe the charge and current distributions 
are 

p(t,x) = density of charge at position x and time t, (1.1.3) 
v(t,x) = velocity of an element of charge at position x and time t,(1.1.4) 
j(t,x) = current density at position x and time t. (1.1.5) 

We shall now establish the important relation 

j = pv. (1.1.6) 

The current density j is defined by the statement 

j ■ da = current crossing an element of area da, (1.1 .7) 
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so that j is the current per unit area. The current flowing across a surface S is then 

j ■ da. 



s 

On the other hand, current is defined as the quantity of charge crossing a surface 
per unit time. Let us then calculate the current associated with a distribution of 
charge with density p and velocity v. 

Figure 1.1 shows a number of charged particles approaching with speed v a small 
surface of area Aa; there is an angle 8 between the direction of the velocity vector 
v and the normal n to the surface. The particles within the dashed box will all hit 
the surface within a time At. The volume of this box is (vAt cos 6) Aa = v ■ nAtAa. 
The current crossing the surface is then calculated as 

current = (charge crossing the surface)/ At 

( volume of box \ 1 

= (total charge) — 

v 6 ; V total volume / At 

= (density) (volume of box) /At 

= p(v ■ nAtAa)/ At 

= pv ■ (nAo) 

= pv ■ Aa. 

Comparing this with the expression given previously, current = j ■ Ax, we find that 
indeed, pv is the current density. 

As was mentioned previously, the fluid description still allows for the existence 
of point charges. For these, the density is zero everywhere except at the position 
of a charge, where it is infinite. This situation can be described by a 5-function. 
Suppose that we have a point charge q at a position r. Its charge density can be 
written as 

p(x) = qS(x - r), 

where S(x) = 5(x)S(y)S(z) is a three-dimensional (^-function. If the charge is moving 
with velocity v, then 

j(x) — qvS(x — r) 



1.2 Maxwell's equations and the Lorentz force 
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is the current density. 

More generally, let us have a collection of point charges gu at positions rA(t), 
moving with velocities VA(t) — drA/dt. Then the charge and current densities of 
the charge distribution are given by 

p(t,x) = ^2qA6{x-r A (t)) (1.1.8) 

A 

and 

j{t,x) = J2<lAVA(t)5(x - r A {t)). (1.1.9) 

A 

Each qa is the integral of p(t, x) over a volume V A that encloses this charge but no 
other: 



= / p(t,x)d 3 x, (1.1.10) 

JVa 



qA 

JVa 

where d 3 x = dxdydz. The total charge of the distribution is obtained by integrating 
the density over all space: 

Q= f P (t lX )d 3 x = J2lA- (1.1-11) 

^ A 

1.2 Maxwell's equations and the Lorentz force 

The four Maxwell equations determine the electromagnetic field once the charge 
and current distributions are specified. They are given by 

V -E = -p, (1.2.1) 

V B = 0, (1.2.2) 

VxE = - — , (1.2.3) 
dE 

VxB = Moj + eoMo^-- (1-2.4) 

Here, eo and are constants, and V = (V x ,V y ,V z ) is the gradient operator 
of vector calculus (with the obvious notation V x / = df/dx for any function /). 
The Maxwell equations state that the electric field is produced by charges and 
time-varying magnetic fields, while the magnetic field is produced by currents and 
time-varying electric fields. Maxwell's equations can also be presented in integral 
form, by invoking the Gauss and Stokes theorems of vector calculus. 

The Lorentz-force law determines the motion of the charges once the electro- 
magnetic field is specified. Let 

f(t, x) = force density at position x and time t, (1.2.5) 

where the force density is defined to be the net force acting on an element of charge 
at x divided by the volume of this charge element. The statement of the Lorentz- 
force law is then 

/ = pE + j x B = p(E + v x B). (1.2.6) 

The net force F acting on a volume V of the charge distribution is the integral of 
the force density over this volume: 

F(t,V) = f f(t,x)d 3 x. (1.2.7) 
Jv 
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For a single point charge we have p = qS(x — r),j = qvS(x — r), and the total force 
becomes 

F = q[E(r)+v x B(r)]. 

Here the fields are evaluated at the charge's position. This is the usual expression 
for the Lorentz force, but the definition of Eq. (1.2.6) is more general. 

Taken together, the Maxwell equations and the Lorentz-force law determine the 
behaviour of the charge and current distributions, and the evolution of the electric 
and magnetic fields. Those five equations summarize the complete theory of classical 
electrodynamics. Every conceivable phenomenon involving electromagnetism can be 
predicted from them. 



1.3 Conservation of charge 

One of the most fundamental consequences of Maxwell's equations is that charge is 
locally conserved: charge can move around but it cannot be destroyed nor created. 

Consider a volume V bounded by a two-dimensional surface S. Charge conser- 
vation means that the rate of decrease of charge within V must be equal to the 
total current flowing across S: 

_d_ 
~~dt 



pd s x = f j ■ da. (1.3.1) 
v Js 



Local charge conservation means that this statement must be true for any volume 
V, however small or large. We may use Gauss' theorem to turn Eq. (1.3.1) into a 
differential statement. For any smooth vector field u within V we have 



/ V • u d 3 x = <p u ■ da. 
Jv Js 



The right-hand side of Eq. (1.3.1) can thus be written as J v V • j d 3 x, while the 
left-hand side is equal to — J v (dp/dt) d 3 x. Equality of the two sides for arbitrary 
volumes V implies 

||+V-j = 0. (1.3.2) 

This is the differential statement of local charge conservation, and we would like to 
prove that this comes as a consequence of Maxwell's equations. 

To establish this we first differentiate Eq. (1.2.1) with respect to time to obtain 

dp _ dE 

where we have exchanged the order with which we take derivatives of the electric 
field. If we now use Eq. (1.2.4) to eliminate the electric field, we get 

^ = -V-(VxB)-V-j. 

at pa 

The first term on the right-hand side vanishes identically (the divergence of a curl is 
always zero) and we arrive at Eq. (1.3.2). We have therefore established that local 
charge conservation is a consequence of Maxwell's equations. 



1.4 Conservation of energy 



[The material presented in this section is also covered in Sec. 6.7 of Jackson's text.] 



1.4 Conservation of energy 
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Other conservation statements follow from Maxwell's equations and the Lorentz- 
force law. In this section we formulate and derive a statement of energy conser- 
vation. In the next section we will consider the conservation of linear momentum. 
It is also possible to prove conservation of angular momentum, but we shall not 
pursue this here. 

Conservation of energy is one of the most fundamental principle of physics, and 
Maxwell's electrodynamics must be compatible with it. A statement of local energy 
conservation can be patterned after our previous statement of charge conservation. 
Let 

e(t,x) = electromagnetic field energy density, (1.4.1) 
S(t,x) ■ da = field energy crossing an element of 

surface da per unit time. (1-4-2) 

So e is analogous to charge density, and S (which is known as Pointing's vector), 
is analogous to current density. 

If electromagnetic field energy were locally conserved, we would write the state- 
ment 

— — / ed 3 x— (f> S-da, 
dt Jv Js 

which is analogous to Eq. (1.3.1). But we should not expect field energy to be 
conserved, and this statement is not correct. The reason is that the field does work 
on the charge distribution, and this takes energy away from the field. This energy 
goes to the charge distribution, and total energy is conserved. A correct statement 
of energy conservation would therefore be (field energy leaving V per unit time) = 
(field energy crossing S per unit time) -I- (work done on charges within V per unit 
time). To figure out what the work term should be, consider an element of charge 
within V. It moves with a velocity v and the net force acting on it is 



fdV = p(E + vx B)dV. 



While the element undergoes a displacement dx, the force docs a quantity of work 
equal to f ■ dx dV. The work done per unit time is then / • v dV = pE ■ v dV = 
j ■ E dV. Integrating this over V gives the total work done on the charges contained 
in V, per unit time. The correct statement of energy conservation must therefore 
have the form 

/ ed 3 x = I S-da + [ j-Ed 3 x. (1.4.3) 
dt J v J s J v 

In differential form, this is 

de 

— + V-S = -j-E. (1.4.4) 

This last equation must be derivable from Maxwell's equations, which have not 
yet been involved. And indeed, the derivation should provide expressions for the 
quantities e and S. What we have at this stage is an educated guess for a correct 
statement of energy conservation, but Eq. (1.4.4) has not yet been derived nor the 
the quantities e and S properly defined in terms of field variables. Now the hard 
work begins. 

We start with the right-hand side of Eq. (1.4.4) and eliminate j in favour of field 
variables using Eq. (1.2.4). This gives 

1 dE 

-j -E = E ■ (V x B) + e E • — . 

Ho dt 
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We use a vector-calculus identity to replace E-(VxB) with B-(VxE) — V-(ExB), 
and in this we replace Vx£ with —dB/dt, using Eq. (1.2.3). All of this gives us 

1 OB dE 1 

-j-E = — B ■ — — + cqE ■ — — H V ■ (E x B) 

dt dt fi 

We have arrived at an equation of the same form as Eq. (1.4.4), and this shows that 
we do indeed have energy conservation as a consequence of Maxwell's equations and 
the Lorentz-force law. 

The calculation also gives us definitions for the field's energy density, 

and for the field's energy flux (Poynting) vector, 

S=—ExB. (1.4.6) 
Mo 

It is an important consequence of Maxwell's theory that the electromagnetic field 
carries its own energy. 



1.5 Conservation of momentum 

[The material presented in this section is also covered in Sec. 6.7 of Jackson's text.] 
We should expect a statement of momentum conservation to take a form similar 
to our previous statement of energy conservation. In that case we had a scalar 
quantity e representing the density of energy, a vectorial quantity S a representing 
the flux of energy, and the conservation statement took the form of 

ds 

ttt + V a 5 a = work term. 
ot 

We have introduced an explicit component notation for vectors, and summation 
over a repeated index is understood. 

For momentum conservation we will need a vectorial quantity e a to represent 
the density of momentum, and a tensorial quantity T a b to represent the flux of 
momentum (one index for the momentum component, another index for the flux 
direction). So let 

e a (t, x) = density of a-component of field momentum, (1.5.1) 
— T a b(t, x) dab = a-component of field momentum crossing 

an element of surface da per unit time. (1.5.2) 

The minus sign in front of T a b is introduced by convention. 

The statement of momentum conservation is (field momentum leaving V per unit 
time) = (field momentum crossing S per unit time) + (rate at which momentum is 
communicated to the charges). Because a rate of change of momentum is a force, 
we have the integral statement 

~Jt S v £a = ~ f s Tah dah + Jv fa d * X ' (L5 ' 3) 
or the equivalent differential statement 

~~dt + VbTab = /a = pEa + (i x B)a - (L5 ' 4) 



1.6 Junction conditions 
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We now would like to derive this from Maxwell's equations, and discover the iden- 
tities of the vector e a and the tensor T ab . 

We first use Eqs. (1.2.1) and (1-2.4) to eliminate p and j from Eq. (1.5.4): 

(1 BE \ 

-VxB-e„- xB 
p dt ) 

= e (V-E)E+—(V-B)B + —(VxB)xB 

Mo Mo 

B ,„ „, „ dB 

-e -(ExB) + e Ex—, 

where to go from the first to the second line we have inserted a term involving 
V • B = and allowed the time derivative to operate on E x B. The next step is 
to use Eq. (1.2.3) to replace dB/dt with —V x E. Collecting terms, this gives 

/ = -e (E x B) + 6 [(V • E)E - E x (V x E)] + — [(V • B)B -Bx(VxB)]. 
at p 

To proceed we invoke the vector-calculus identity \V{u-u) = (u- V)m+mx (Vxu), 
where u stands for any vector field. We use this to clean up the quantities within 
square brackets. For example, 

(V • E)E - E x (V x E) = (V • E)E + (E ■ V)E - ■ E). 

In components, the right-hand side reads 

(V b E b )E a + (E b V b )E a - \^aE 2 = E a V b E b + E b V b E a - h ab V b E 2 

= Vfc ^E a E b - -S ab E 2 

Doing the same for the bracketed terms involving the magnetic field, we arrive at 

fa = -^ t ( E x B )a + e V b (^E a E b - \s ab E 2 ^j + ±-V b (^B a B b - ^S ab B : 

This has the form of Eq. (1.5.4) and we conclude that momentum conservation does 
indeed follow from Maxwell's equations and the Lorentz-force law. 

Our calculation also tells us that the field's momentum density is given by 

e a = e Q (E x B) a , (1.5.5) 

and that the field's momentum flux tensor is 

T ab = e (E a E b - \s ab E 2 ^j + i- (l3 a B b - \s ab B 2 ^j . (1.5.6) 

Notice that the momentum density is proportional to the Poynting vector: e — 
eqPqS. This means that the field's momentum points in the same direction as the 
flow of energy and that these quantities are related by a constant factor of 6qPq. 



1.6 Junction conditions 

[The material presented in this section is also covered in Sec. 1.5 of Jackson's text.] 
Maxwell's equations determine how the electric and magnetic fields must be 
joined at an interface. The interface might be just a mathematical boundary (in 
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which case nothing special should happen), or it might support a surface layer 
of charge and/or current. Our task in this section is to formulate these junction 
conditions. 

To begin, we keep things simple by supposing that the interface is the x-y plane 
at z — — a nice, flat surface. If the interface supports a charge distribution, then 
the charge density must have the form 

p(t, x) = p+(t, x)6{z) + p_(t, x)6(-z) + a(t, x, y)S(z), (1.6.1) 

where p + is the charge density above the interface (in the region z > 0), p_ the 
charge density below the interface (in the region z < 0), and a is the surface density 
of charge (charge per unit surface area) supported by the interface. 

We have introduced the Heaviside step function 9(z), which is equal to 1 if 
z > and otherwise, and the Dirac (^-function, which is such that J f(z)S(z) dz — 
/(0) for any smooth function f(z). These distributions (also known as generalized 
functions) are related by the identities 

±e(z) = 5(z), ±e { - z ) = S(z). (1.6.2) 

As for any distributional identity, these relations can be established by integrating 
both sides against a test function f(z). (Test functions are required to be smooth 
and to fall off sufficiently rapidly as z — > ±oo.) For example, 

-f(z)dz = - / 0(z)— — dz 



dz ' .l—r*z dz 

df 

/(0), 



o 



and this shows that d0(z)/dz is indeed distributionally equal to S(z). You can show 
similarly that f(z)6(z) = f{0)S{z) and f(z)S'(z) = f{°) S '( z ) - f'(°) S ( z ) are valid 
distributional identities (a prime indicates differentiation with respect to z). 

If the interface at z = also supports a current distribution, then the current 
density must have the form 

j(t, x) = j+(t, x)6(z) + j (t, x)0{-z) + K(t, x, y)S(z), (1.6.3) 

where j + is the current density above the interface, j_ the current density below 
the interface, and K is the surface density of current (charge per unit time and unit 
length) supported by the interface. 

We now would like to determine how E and B behave across the interface. We 
can express the fields as 

E(t,x) = E+(t,x)6(z) + E_(t,x)6(-z), 

(1.6.4) 

B(t,x) = B+(t,x)6(z) + B_(t,x)6(-z), 

in an obvious notation; for example, E + is the electric field above the interface. We 
will see that these fields are solutions to Maxwell's equations with sources given by 
Eqs. (1.6.1) and (1.6.3). To prove this we need simply substitute Eq. (1.6.4) into 
Eqs. (1.2.1)-(1.2.4). Differentiating the coefficients of 0(±z) is straightforward, but 
we must also differentiate the step functions. For this we use Eq. (1.6.2) and write 



V9(±z) =±S(z)z, 



1.6 Junction conditions 
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where £ is a unit vector that points in the z direction. 
A straightforward calculation gives 

V -E = (V-E + )6(z) + (V-E_)6{-z) + z-(E + -E_)6(z), 

V x E = (V x E+)6{z) + (V x E_)6{-z) + z x (E + - E_)S{z), 

V B = (V -B + )6{z) + (V -B_)9(-z) + z- (B + -B_)5{z), 
VxB = (V x B+)9(z) + (V x B-)6(-z) + zx (B+ - B_)6{z). 

Substituting Eq. (1.6.1) and the expression for V • E into the first of Maxwell's 
equations, V • E = p/eo, reveals that E± is produced by p± and that the surface 
charge distribution creates a discontinuity in the normal component of the electric 
field: 

z-{E + -E-)\ z=0 = ±v. (1.6.5) 

Substituting the expression for V-JB into the second of Maxwell's equations, V £? = 
0, reveals that B± satisfies this equation on both sides of the interface and that the 
normal component of the magnetic field must be continuous: 

z-(B + -B_)\ z=0 = 0. (1.6.6) 

Substituting Eq. (1.6.4) and the expression for V x E into the third of Maxwell's 
equations, Vx£= —dB/dt, reveals that E± and B± satisfy this equation on both 
sides of the interface and that the tangential (x and y) components of the electric 
must also be continuous: 



zx (E + -E_)\ z=a = 0. (1.6.7) 

Finally, substituting Eqs. (1.6.3), (1.6.4), and the expression for V x B into the 
fourth of Maxwell's equations, V x B = p j + e podE/dt, reveals that B± is 
produced by j± and that the surface current creates a discontinuity in the tangential 
components of the magnetic field: 

zx (B + -B_)\ z=0 = p K. (1.6.8) 

Equations (1.6.5)-(1.6.8) tell us how the fields E± and B± are to be joined at a 
planar interface. 

It is not difficult to generalize this discussion to an interface of arbitrary shape. 
All that is required is to change the argument of the step and S functions from z 
to s(x), where s is the distance from the interface in the direction normal to the 
interface; this is positive above the interface and negative below the interface. Then 
Vs = n, the interface's unit normal, replaces z in the preceding equations. If E + , 
B + are the fields just above the interface, and E_, B the fields just below the 
interface, then the general junction conditions are 



h-(E+-E) = — cr, (1.6.9) 

n-(S+-B_) = 0, (1.6.10) 

nx (E+-E-) = 0, (1.6.11) 

nx(B + -B_) = p K. (1.6.12) 



Here, a is the surface charge density on the general interface, and K is the surface 
current density. It is understood that the normal vector points from the minus side 
to the plus side of the interface. 
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1.7 Potentials 

[The material presented in this section is also covered in Sees. 6.2 and 6.3 of Jack- 
son's text.] 

The Maxwell equation V • B — and the mathematical identity V • (V x A) = 
imply that the magnetic field can be expressed as the curl of a vector potential 
A(t,x): 

B = VxA. (1.7.1) 

Suppose that B is expressed in this form. Then the other inhomogeneous Maxwell 
equation, V x E + dB/dt = 0, can be cast in the form 

Vx( E+ ^).0. 

This, together with the mathematical identity V x (V$) = 0, imply that E+dA/dt 
can be expressed as the divergence of a scalar potential $(t,cc): 

8A 

E = -— -V$; (1.7.2) 

dt K ' 

the minus sign in front of V$ is conventional. 

Introducing the scalar and vector potentials eliminates two of the four Maxwell 
equations. The remaining two equations will give us equations to be satisfied by the 
potentials. Solving these is often much simpler than solving the original equations. 
An issue that arises is whether the scalar and vector potentials are unique: While 
E and B can both be obtained uniquely from $ and A, is the converse true? The 
answer is in the negative — the potentials are not unique. 

Consider the following gauge transformation of the potentials: 

$ - $ocw = $-|, (1.7.3) 
A - A new = A + Vf, (1.7.4) 

in which f(t,x) is an arbitrary function of space and time. We wish to show that 
this transformation leaves the fields unchanged: 

^new = E, B ncw = B. (1.7.5) 

This implies that the potentials are not unique: they can be redefined at will by a 
gauge transformation. The new electric field is given by 

j? _ <9A ncw 

^ncw — v *new 



_!( A + v,)-v(.-§£) 



dt J+ dt' 

the last two terms cancel out and we have established the invariance of the electric 
field. On the other hand, the new magnetic field is given by 

B ncw V x A ncw 

= V x (A + V/) 
= B + Vx(V/); 

the last term is identically zero and we have established the invariance of the mag- 
netic field. We will use the gauge invariance of the electromagnetic field to simplify 
the equations to be satisfied by the scalar and vector potentials. 



1.7 Potentials 
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To derive these we start with Eq. (1.2.1) and cast it in the form 
I,-V.«-V.(-£-V*)~|v.A-V*, 

or 

-P = eo,^ -WH--[V-A + £0M0 - j , (1.7.6) 

where, for reasons that will become clear, we have added and subtracted a term 
e ^d 2 ^/dt 2 . Similarly, we re-express Eq. (1.2.4) in terms of the potentials: 



dE 

Hoj = VxB-eoHo — 



d ( dA \ 
= Vx(VxA)- £0M0 -^ ¥ -V$j 

d 2 A <9<I> 
= V(V-A)-V 2 A + e Aio^-+eoMoV — , 

where we have used a vector-calculus identity to eliminate V x (V x A). We have 
obtained 

d 2 A ( <9<I>\ 

Mo 3 = eoMo - V 2 A + V ( V • A + e oMo — ). (1.7.7) 

Equations (1.7.6) and (1.7.7) govern the behaviour of the potentials; they are equiv- 
alent to the two Maxwell equations that remain after imposing Eqs. (1.7.1) and 
(1.7.2). 

These equations would be much simplified if we could demand that the potentials 
satisfy the supplementary condition 

V-A + e Mo^=0. (1.7.8) 

This is known as the Lorenz gauge condition, and as we shall see below, it can 
always be imposed by redefining the potentials according to Eqs. (1.7.3) and (1-7.4) 
with a specific choice of function f(t,x). When Eq. (1.7.8) holds we observe that 
the equations for $ and A decouple from one another, and that they both take the 
form of a wave equation: 

1_ cP_ 
'^dt 2 



where 



V 2 )A(t,x) = -fx j(t,x), (1.7.10) 

(1.7.11) 



1 

^/mM) 

is the speed with which the waves propagate; this is numerically equal to the speed 
of light in vacuum. 

To see that the Lorenz gauge condition can always be imposed, imagine that 
we are given potentials <J> id and A id that do not satisfy the gauge condition. We 
know that we can transform them according to $ Q id — > $ = ^oid — df/dt and 
A \a — > A = A id + V/, and we ask whether it is possible to find a function f(t, x) 
such that the new potentials $ and A will satisfy Eq. (1.7.8). The answer is in the 
affirmative: If Eq. (1.7.8) is true then 

= V-A + e Q ^ — 



12 



Maxwell's electrodynamics 



= V • (Am + V/) + e fi — I $ old - -± j 

= V • A old + e Mo— ^ e °^ ^ + V 

and we see that the Lorenz gauge condition is satisfied if f(t, x) is a solution to the 
wave equation 

i a 2 2 \ 9$ old 

e2^2+V J/ = -V.A old -eoMo^ r . 

Such an equation always admits a solution, and we conclude that the Lorenz gauge 
condition can always be imposed. (Notice that we do not need to solve the wave 
equation for /; all we need to know is that solutions exist.) 

To summarize, we have found that the original set of four Maxwell equations for 
the fields E and B can be reduced to the two wave equations (1.7.9) and (1.7.10) for 
the potentials $ and A, supplemented by the Lorenz gauge condition of Eq. (1.7.8). 
Once solutions to these equations have been found, the fields can be constructed 
with Eqs. (1.7.1) and (1.7.2). Introducing the potentials has therefore dramatically 
reduced the complexity of the equations, and correspondingly increased the ease of 
finding solutions. 

In time-independent situations, the fields and potentials no longer depend on t, 
and the foregoing equations reduce to 

V 2 $(cc) = -— p(x), (1.7.12) 
V 2 A{x) = -noj(x), (1.7.13) 

as well as E = — V$ and B = V x A. The Lorenz gauge condition reduces to 
V • A = (also known as the Coulomb gauge condition), and the potentials now 
satisfy Poisson's equation. 



1.8 Green's function for Poisson's equation 

[The material presented in this section is also covered in Sec. 1.7 of Jackson's text.] 
We have seen in the preceding section that the task of solving Maxwell's equa- 
tions can be reduced to the simpler task of solving two wave equations. For time- 
independent situations, the wave equation becomes Poisson's equation. In this and 
the next section we will develop some of the mathematical tools needed to solve 
these equations. 

We begin in this section with the mathematical problem of solving a generic 
Poisson equation of the form 

V 2 V(x) = -4tt/(x), (1.8.1) 

where i/j(x) is the potential, f(x) a prescribed source function, and where a factor 
of 4-7T was inserted for later convenience. 

To construct the general solution to this equation we shall first find a Green's 
function G(x 7 x') that satisfies a specialized form of Poisson's equation: 

V 2 G(x,x') = -4tv5(x-x'), (1.8.2) 

where S(x — x') is a three-dimensional Dirac 5-function; the source point x' is 
arbitrary. It is easy to see that if we have such a Green's function at our disposal, 
then the general solution to Eq. (1.8.1) can be expressed as 



4>(x)=Mx)+ / G(x,x')f(x')d 3 x', 



(1.8.3) 



1.8 Green's function for Poisson's equation 
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where ipo{ x ) is a solution to the inhomogcneous (Laplace) equation, 

V 2 ip (x)=0. (1.8.4) 

This assertion is proved by substituting Eq. (1.8.3) into the left-hand side of Eq. (1.8.1), 
and using Eqs. (1.8.2) and (1.8.4) to show that the result is indeed equal to — Airf(x). 

In Eq. (1.8.3), the role of the integral is to account for the source term in 
Eq. (1.8.1). The role of ipoi x ) is to enforce boundary conditions that we might wish 
to impose on the potential ip(x). In the absence of such boundary conditions, ipo 
can simply be set equal to zero. Methods to solve boundary-value problems will be 
introduced in the next chapter. 

To construct the Green's function we first argue that G{x,x') can depend only 
on the difference x — x'- this follows from the fact that the source term depends only 
on x — x' , and the requirement that the Green's function must be invariant under 
a translation of the coordinate system. We then express G(x,x') as the Fourier 
transform of a function G(k): 

G(x,x') = -J--r { G(k)e lk < x - x 'U 3 k, (1.8.5) 
(2ir) 6 J 

where A; is a vector in reciprocal space; this expression incorporates our assumption 
that G(x,x') depends only on x — x'. Recalling that the three-dimensional d- 
function can be represented as 

5{x - x') = J e lk - (x - x "> d 3 fc, (1.8.6) 

we see that Eq. (1.8.2) implies k 2 G(k) = An, so that Eq. (1.8.6) becomes 

G^) = ^J—^-^k. (1.8.7) 

We must now evaluate this integral. 

To do this it is convenient to switch to spherical coordinates in reciprocal space, 
defining a radius k and angles 9 and <f> by the relations k x = k sin 9 cos <j>, k y = 
k sin 9 sin (j>, and k z — kcos9. In these coordinates the volume element is d 3 k = 
k 2 sin 9 dkd9d(f>. For convenience we orient the coordinate system so that the k z 
axis points in the direction of R = x — x' . We then have k R = kRcos9, where 
R = \R\ is the length of the vector R. All this gives 

G(x,x') = — / dk d9 d(j>e ikRcosf> sm6. 
2t Jo Jo Jo 

Integration over d<j> is trivial, and the 9 integral can easily be evaluated by using 
/i = cos 9 as an integration variable: 



J e ^Rcose s - mdd& ^ J e ikR» dfl = 



2 sin kR 
kR ■ 



We now have 



7T J Q kR 7T i? J £ 

where we have switched to £ = kR. The remaining integral can easily be evaluated 
by contour integration (or simply by looking it up in tables of integrals) ; it is equal 
to 7r/2. 
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Our final result is therefore 

= d.8.8) 

the Green's function for Poisson's equation is simply the reciprocal of the distance 
between x and x'. The general solution to Poisson's equation is then 

i>(x) = M*) + [ /^X d 3 x'. (1.8.9) 
J \x x I 

The meaning of this equation is that apart from the term ?p , the potential at x is 
built by summing over source contributions from all relevant points x' and dividing 
by the distance between x and x' . Notice that in Eq. (1.8.8) we have recovered a 
familiar result: the potential of a point charge at x' goes like 1/R. 



1.9 Green's function for the wave equation 

[The material presented in this section is also covered in Sec. 6.4 of Jackson's text.] 
We now turn to the mathematical problem of solving the wave equation, 

OI>{t,x) = -4*f(t,x), (1.9.1) 

for a time-dependent potential ip produced by a prescribed source /; we have intro- 
duced the wave, or d'Alembertian, differential operator 

□ (1.9.2) 

For this purpose we seek a Green's function G(t, x; t', x') that satisfies 

UG{t, x; t', x') = -^5{t - t')8(x - x'). (1.9.3) 
In terms of this the general solution to Eq. (1.9.2) can be expressed as 

i>(t,x) = i> (t,x) + J j ' G(t,x;t' \x')f(t' ',x')dt'd 3 x', (1.9.4) 

where ipo(t, x) is a solution to the homogeneous wave equation, 

ni> (t,x) = Q. (1.9.5) 

That the potential of Eq. (1.9.4) does indeed solve Eq. (1.9.1) can be verified by 
direct substitution. 

To construct the Green's function we Fourier transform it with respect to time, 

G(t, x; t', x') = ^- { G{uo; x, x^e'^-^ dcu, (1.9.6) 
2tt J 

and we represent the time (5-function as 

5{t -t') = ^- [ e -^(*- 4 ') duo. (1.9.7) 

27T J 

Substituting these expressions into Eq. (1.9.3) yields 

[V 2 + (lu/c) 2 ] G(u>; x, x') = -4tt5(x - x'), (1.9.8) 

which is a generalized form of Eq. (1.8.2). From this comparison we learn that 
G(0;x,x') = l/\x - x'\. 



1.9 Green's function for the wave equation 
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At this stage we might proceed as in Sec. 1.8 and Fourier transform G(u>; x, x') 
with respect to the spatial variables. This would eventually lead to 

47T f e ik i x - x ') 



An f e ll * yx-x ) 
G(.;x,x) = — 3 j k2 _ {uj/c)2 dk 



which is a generalized form of Eq. (1.8.7). This integral, however, is not defined, 
because of the singularity at k 2 = (lu/c) 2 . We shall therefore reject this method of 
solution. (There are ways of regularizing the integral so as to obtain a meaningful 
answer. One method involves deforming the contour of the u> integration so as to 
avoid the poles. This is described in Sec. 12.11 of Jackson's text.) 
We can anticipate that for u) ^ 0, G will be of the form 

G(u;x,x') = ^l*-f) , (L9 . 9) 
\x - x'\ 

with g representing a function that stays nonsingular when the second argument, 
R = \x — x'\, approaches zero. That G should depend on the spatial variables 
through R only can be motivated on the grounds that three-dimensional space is 
both homogeneous (so that G can only depend on the vector R = x — x') and 
isotropic (so that only the length of the vector matters, and not its direction). That 
G should behave as 1/R when R is small is justified by the following discussion. 

Take Eq. (1.9.8) and integrate both sides over a sphere of small radius e centered 
at x' . Since V 2 G = V • VG, we can use Gauss' theorem to get 

I (VG-R)da+(uj/c) 2 [ Gd 3 x = -An, 

JR=e JR<s 

where R = R/R. In this equation, the volume integral is of order Ge 3 and it 
contributes nothing in the limit e — ► 0, unless G happens to be as singular as 1/e 3 . 
The surface integral, on the other hand, is equal to 



A-KE 2 — 

dR 



R=s 



If G were to behave as 1/e 3 , then dG/dR would be of order 1/e 4 , the surface integral 
would contribute a term of order 1/e 2 , and the left-hand side would never give rise 
to the required —An. So we conclude that G cannot be so singular, and that the left- 
hand side is dominated by the surface integral. This implies that G must be of order 
1/e, as was expressed by Eq. (1.9.9). Setting G = g/R returns —Ang(ui, e) + 0(e) 
for the surface integral, and this gives us the condition 

ff(w,0) = l. (1.9.10) 

We also recall that g(0, R) = 1. 

We may now safely take R ^ and substitute Eq. (1.9.9) into Eq. (1.9.8), taking 
its right-hand side to be zero. Since G depends on x only through R, the Laplacian 
operator becomes 

R 2 dR dR' 

Acting with this on G = g/R yields g" / R and Eq. (1.9.8) becomes 

g" + {u/cfg = 0, (1.9.11) 

with a prime indicating differentiation with respect to R. With the boundary con- 
dition of Eq. (1.9.10), two possible solutions to this equation are 

g±(w,R) = e ±l ^/^ R . (1.9.12) 



16 



Maxwell's electrodynamics 



Substituting this into Eq. (1.9.9), and that into Eq. (1.9.6), we obtain 

G ± (t,x;t',x') = ±J e±< ^e- i "^cko = ^ J e -M*-*'=F*/=) ^ 

or 

, /^ S(t-t'T\x-x'\/c) 

G±(t, x; t', x') = -± 1 ' 7 . U } . 1.9.13 

\x — x'\ 

These are the two fundamental solutions to Eq. (1.9.3). The function G + (t, x; t', x'), 
which is nonzero when t — t' = +R/c, is known as the retarded Green's function; 
the function G_ (t, x; t', x'), which is nonzero when t — t' = —R/c, is known as the 
advanced Green's function. 

Substituting the Green's functions into Eq. (1.9.4) gives 

1>±(t,x)=Mt,*) + J / ^'^["'^ fit'^dt'dZx'. 
The time integration can be immediately carried out, and we obtain 

, s / s f f(t T \ x - x'\/c,x') o , 
V±(i, x) = Mt, x)+J n \ X _J( ' d 3 x' (1.9.14) 

for the fundamental solutions to the wave equation. Except for the shifted time 
dependence, Eq. (1.9.14) looks very similar to Eq. (1.8.9), which gives the general 
solution to Poisson's equation. The time translation, however, is very important. 
It means that the potential at time t depends on the conditions at the source at a 
shifted time 

t' = t =F | a; - x'\/c. 

The second term on the right-hand side is the time required by light to travel the 
distance between the source point x' and the field point x; it encodes the property 
that information about the source travels to x at the speed of light. The "+" 
solution, ip+(t,x), depends on the behaviour of the source at an earlier time t — R / c 
— there is a delay between the time the information leaves the source and the time 
it reaches the observer. This solution to the wave equation is known as the retarded 
solution, and it properly enforces causality: the cause (source) precedes the effect 
(potential). On the other hand, ip-(t, x) depends on the behaviour of the source at 
a later time t + R/c; here the behaviour is anti-causal, and this unphysical solution 
to the wave equation is known as the advanced solution. 

Although both solutions to the wave equation are mathematically acceptable, 
causality clearly dictates that we should adopt the retarded solution ip+(t, x) as the 
only physically acceptable solution. So we take 

^{t, x) = Mt, x) + [ fit-fc-vyc'*') d 3 x , (1 9 15) 

J \x — x'\ 

to be the solution to the wave equation (1.9.1). Recall that tpo(t, x) is a solution to 
the homogeneous equation, Dtpo = 0. 



1.10 Summary 

The scalar potential $ and the vector potential A satisfy the wave equations 

(-^ + V 2 )^,x) = -lp(^) (1.10.1) 



1.11 Problems 
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and 

1 d 2 



. (/2 +W 2 jA(t,x) = -p j{t,x) (1.10.2) 
when the Lorenz gauge condition 

1 <9$ 

V-A + -- = (1.10.3) 

is imposed. The speed of propagation of the waves, c = l/y / eo/Io, is numerically 
equal to the speed of light in vacuum. 

The retarded solutions to the wave equations arc 

1 f p(t — \x — x'\/c,x') , , , 
sb) = $o(t, as) + i / Tr"^ d x ( L10 - 4 ) 

47T£o J 1 35 — X'\ 

and 

A(t, x) = A (t, x) + ^( ^ (L10 . 5) 

47r J |aj — x'\ 

In static situations, the time dependence of the sources and potentials can be 
dropped. 

Once the scalar and vector potentials have been obtained, the electric and mag- 
netic fields are recovered by straightforward differential operations: 

BA 

E = -— -V$, B = VxA. (1.10.6) 

This efficient reformulation of the equations of electrodynamics is completely equiva- 
lent to the original presentation of Maxwell's equations. It will be the starting point 
of most of our subsequent investigations. 



1.11 Problems 

1. We have seen that the charge density of a point charge q located at r(t) is 
given by 

p(x) = q8(x - r(t)), 

and that its current density is 

j(x) = qv(t)S(x - r(t)), 

where v = dr/dt is the charge's velocity. Prove that p and j satisfy the 
statement of local charge conservation, 

2. In this problem we explore the formulation of electrodynamics in the Coulomb 
gauge, an alternative to the Lorenz gauge adopted in the text. The Coulomb 
gauge is especially useful in the formulation of a quantum theory of electro- 
dynamics. 

a) Prove that the Coulomb gauge condition, W A = 0, can always be imposed 
on the vector potential. Then show that in this gauge, the scalar potential 
satisfies Poisson's equation, V 2 $ = —p/e , so that it can be expressed as 

P{t,x') 



47re J \x — x'\ 



Show also that the vector potential satisfies the wave equation 
OA = -p j + e a p V — . 
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b) Define the "longitudinal current" by 

. , \ 1 „ f V-j(t,x') ,, , 
4tt J \x — x'\ 

and show that this satisfies V x j\ = 0. Then prove that the wave 
equation for the vector potential can be rewritten as DA = — Hojt, where 
j t = j — j\ is the "transverse current" . 

c) Prove that the transverse current can also be defined by 

*M-i T "( T «/^ W )' 

and that it satisfies V • jt — 0. 

d) The labels "longitudinal" and "transverse" can be made more meaningful 

by formulating the results of part b) and c) in reciprocal space instead of 
ordinary space. For this purpose, introduce the Fourier transform j(t, k) 
of the current density, such that 

j(t,x) = ^ J j(t,k)e ik *d 3 k. 

Then prove that the Fourier transform of the longitudinal current is 

ji(t,k) = [k-j(t,k)]k, 

where k = k/\k\ is a unit vector aligned with the reciprocal position 
vector k; thus, the longitudinal current has a component in the direction 
of k only. Similarly, prove that the Fourier transform of the transverse 
current is 

jt(t,k) = [fe x j{t,k)] x k; 

thus, the transverse current has components in the directions orthogo- 
nal to k only. For a concrete illustration, assume that k points in the 
direction of the z axis. Then show that j\ — j z z and j t = j x x + j y y. 



3. Prove that 



is a solution to 



G(x,x') = ^ e ik\x-x'\ 



d^_ 

dx 2 



k 2 ^jG{x,x')=5(x-x'), 



where k is a constant. This shows that G(x, x') is a Green's function for the 
inhomogeneous Helmholtz equation in one dimension, (d 2 /dx 2 + k 2 )ip(x) = 
f(x), where tp is the potential and / the source. 

4. The differential equation 

$+^ = /«) 

governs the motion of a simple harmonic oscillator of unit mass and natural 
frequency lo driven by an arbitrary external force f(t). It is supposed that 
the force starts acting at t = 0, and that prior to t = the oscillator was at 
rest, so that x(0) = x(0) = 0, with an overdot indicating differentiation with 
respect to t. 
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Find the retarded Green's function G(t, t') for this differential equation, which 
must be a solution to 

Then find the solution x(t) to the differential equation (in the form of an 
integral) that satisfies the stated initial conditions. 

Check your results by verifying that if f(t) = 6(t) cos(wi), so that the oscillator 
is driven at resonance, then x(t) = (2u))~ 1 t sinwt. 

5. In this problem (adapted from Jackson's problem 6.1) we construct solutions 
to the inhomogeneous wave equation Dip = —iirf in a few simple situations. 
A particular solution to this equation is 

ip{t,x)= J G(t,x;t',x')f{t',x')dt'd 3 x', 

where 

G(w ^=^- t ;- | " ,|/c) 

\x - x'\ 

is the retarded Green's function. 

a) Take the source of the wave to be a momentary point source, described 

by f(f,x') = 5(x')S(y')6(z')6(t'). Show that the solution to the wave 
equation in this case is given by 

, 5(t-r/c) 

ip(t,x,y,z) = , 

r 

where r = \/x 2 +y 2 + z 2 . The wave is therefore a concentrated wave- 
front expanding outward at the speed of light. 

b) Now take the source of the wave to be a line source described by f(t',x') = 

5(x')8(y')8{t'). (Notice that this also describes a point source in a space 
of two dimensions.) Show that the solution to the wave equation in this 
case is given by 

2c6(t - p/c) 



ip(t,x,y) 



^{ct) 2 -p 21 



where p = \J x 2 + y 2 and 9(s) is the Heaviside step function. Notice that 
the wave is no longer well localized, but that the wavefront still travels 
outward at the speed of light. 

c) Finally, take the source to be a sheet source described by f(t', x') = 
S(x')5(t'). (Notice that this also describes a point source in a space 
of a single dimension.) Show that the solution to the wave equation in 
this case is given by 

ip(t,x) = 2irc9(t- \x\/c). 

Notice that apart from a sharp cutoff at \x\ — ct, the wave is now com- 
pletely uniform. 



6. In the quantum version of Maxwell's theory, the electromagnetic interaction 
is mediated by a massless photon. In this problem we consider a modified 
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classical theory of electromagnetism that would lead, upon quantization, to a 
massive photon. It is based on the following set of field equations: 

V -E + k 2 <5> 
V B 
V x E 

V x B + k 2 A 

where k is a new constant related to the photon mass to by k = mc/h, and 
where the fields E and B arc related in the usual way to the potentials $ and 
A. The field equations are supplemented by the same Lorentz-force law, 

f = P E + j x B, 

as in the original theory. 

a) Prove that the modified theory enforces charge conservation provided that 

the potentials are linked by 

^ d$ n 

V- A + c^q— =0. 

The Lorenz gauge condition must therefore always be imposed in the 
modified theory. 

b) Find the modified wave equations that are satisfied by the potentials $ 

and A. 

c) Verify that in modified electrostatics, the scalar potential outside a point 

charge q at x — is given by the Yukawa form 



47reo r 

d) Prove that the modified theory enforces energy conservation by deriving 
an equation of the form 

% + V.S = -j.E. 

Find the new expressions for e and S; they should reduce to the old 
expressions when k — > 0. In the original theory the potentials have no 
direct physical meaning; is this true also in the modified theory? 



= 0, 

= Moj + eoMo 



dE 

at' 



Chapter 2 
Electrostatics 



2.1 Equations of electrostatics 

For time-independent situations, the equations of electromagnetism decouple into 
a set of equations for the electric field alone — electrostatics — and another set for 
the magnetic field alone — magnetostatics. The equations of magnetostatics will 
be considered in Chapter 3. The topic of this chapter is electrostatics. 

The equations of electrostatics are Poisson's equation for the scalar potential, 

W 2 <f(x) = --p(x), (2.1.1) 

and the relation between potential and field, 

E = -V<f>. (2.1.2) 
The general solution to Poisson's equation is 

* W = toW + J_/_^V, (2.1.3, 

where &o( x ) satisfies Laplace's equation, 

V 2 $oOz) = 0. (2.1.4) 

The role of &o( x ) is to enforce boundary conditions that we might wish to impose 
on the scalar potential. 



2.2 Point charge 

The simplest situation involves a point charge q located at a fixed point b. The 
charge density is 

p(x) = q5(x-b), (2.2.1) 
and in the absence of boundaries, Eq. (2.1.3) gives 

®( x ) = T i —r JL -u- ( 2 - 2 - 2 ) 

v ; 47re \x-b\ y ' 

To compute the electric field we first calculate the gradient of R = \x — b\, the 
distance between the field point x and the charge. We have R 2 = (x — b x ) 2 + 
(y — by) 2 + (z — b z ) 2 and differentiating both sides with respect to, say, x gives 
2RdR/dx = 2(x — b x ), or dR/dx = (x — b x )/R. Derivatives of R with respect 
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Figure 2.1: Geometry of a dipole. 



to y and z can be computed in a similar way, and we have established the useful 
vectorial relation 

Notice that Vi? is a unit vector that points in the direction of R = x — b. 
The calculation of the electric field involves Vi? -1 = — R~ 2 WR, or 



1 



\x-b\ \x-b\ 3 ' 

This gives 



(2.2.4) 



E{x) = -^—, (2.2.5) 

47re \x - b| 3 

The electric field goes as q/R 2 and points in the direction of the vector R = x — b. 
This well-known result is known as Coulomb 's law. 



2.3 Dipole 

Another elementary situation is that of two equal charges, one positive, the other 
negative, separated by a distance d = 2e. We align the charges along the direction 
of the unit vector ft, and place the origin of the coordinate system at the middle 
point; see Fig. 2.1. 

The charge density of this distribution is given by 



p{x) = qS(x — eft) — qS(x + eft), (2.3.1) 

and the potential is 

®( x ) = ir-(- l — ^r-i — ( 2 - 3 - 2 ) 
47re \\x-ep\ \x + ep\J 

This is an exact expression. We can simplify it if we take the observation point x 
to be at a large distance away from the dipole, 



r = \x\ > e. 



(2.3.3) 
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We can then approximate 
\xT£p\ 



or 



1 



\fr 2 =F 2ep • x + e 2 
1 =F ep ■ x/r 2 + 0(e 2 /r 2 ) 

l±ep-x/r 2 + 0(e 2 /r 2 ) 



\x T ep\ 

Keeping terms of order e only, the potential becomes 

1 (2qep) ■ x 



47re 



The vector p = 2qep is the dipole moment of the charge distribution. In general, 
for a collection of several charges, this is defined as 



p = ^2qAX A , 



(2.3.4) 



where xa is the position vector of the charge In the present situation the 
definition implies p = (+q)(ep) + (—q)(—sp) — 2qep, as was stated previously. The 
potential of a dipole is then 



$(sc) 



1 p ■ X 

Aneo r 3 



(2.3.5) 



when r satisfies Eq. (2.3.3). It should be noticed that here, the potential falls off as 
1/r 2 , faster than in the case of a single point charge. This has to do with the fact 
that here, the charge distribution has a vanishing total charge. 

It is instructive to rederive this result by directly expanding the charge density 
in powers of s. For any function / of a vector x we have f(x + ep) = f(x) + ep ■ 
V/(cc) + 0(e 2 ), and we can extend this result to Dirac's distribution: 

5(x ± ep) = 5{x) ± ep ■ V5(x) + 0(e 2 ). 

The charge density of Eq. (2.3.1) then becomes p{x) = —2qep ■ \ 7 S(x) + 0(e 2 ), or 

p(x) = -p-VS(x). (2.3.6) 

This is the charge density of a point dipole. 

Substituting this into Eq. (2.1.3) and setting $ = yields 



<f>(x) 



= Pa f 
47re J 



a x , 



IX — X' 



where summation over the vectorial index a is understood. The standard procedure 
when dealing with derivatives of a (^-function is to integrate by parts. The integral 
is then 



.3 ' 



and we have recovered Eq. (2.3.5). Notice that in these manipulations we have 
inserted a result analogous to Eq. (2.2.4); the minus sign does not appear because 
we are now differentiating with respect to the primed variables. 
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Figure 2.2: An inner boundary S{ n on which $ is specified, an outer boundary 
S ou t on which <I> is specified, and the region V in between that contains a charge 
distribution p. The outward unit normal to both boundaries is denoted n. 



From the potential of Eq. (2.3.5) we can calculate the electric field of a point 
dipole. Taking a derivative with respect to x gives 



1 



47Te 



V x (p-x) 3{p-x) 



We have V x (p ■ x) = p x and according to Eq. (2.2.3), V^r = x/r. So 



1 

4ire 



p x 3(p • x)x 



and a similar calculation can be carried out for the y and z components of V$. 
Introducing the unit vector 



r = x/r = (x/r,y/r, z/r), 
we find that the electric field is given by 

1 3(p • f)f — p 



E(x) 



47Te 



(2.3.7) 



(2.3. 



This is an exact expression for the electric field of a point dipole, for which the 
charge density is given by Eq. (2.3.6). For a physical dipole of finite size, this 
expression is only an approximation of the actual electric field; the approximation 
is good for r>£. 



2.4 Boundary- value problems: Green's theorem 

[The material presented in this section is also covered in Sees. 1.8 and 1.10 of 
Jackson's text.] 

After the warmup exercise of the preceding two sections we are ready to face 
the more serious challenge of solving boundary-value problems. A typical situation 
in electrostatics features a charge distribution p(x) between two boundaries on 
which the potential $(x) is specified (see Fig. 2.2) — the potential is then said to 
satisfy Dirichlet boundary conditions. A concrete situation might involve an inner 
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boundary S- ln that is also a grounded conducting surface (on which $ = 0) and 
an outer boundary S ou t that is pushed to infinity (so that $ = there also). We 
would like to derive an expression for <!>(a;) in the region between the boundaries; 
this should account for both the charge density and the boundary data. To arrive 
at this expression we will need to introduce a new type of Green's function, which 
we will call a Dirichlet Green's function and denote G^{x, x'); its properties will be 
identified along the way. To get there we will make use of Green's identity, which 
is essentially an application of Gauss' theorem. 

Gauss' theorem states that for a volume V enclosed by a surface S, 

V • bd 3 x = f b- da, 
v Js 

where b is any vector field defined within V and da is an outward surface element 
on S. The theorem, as stated, says nothing about a region V bounded by two 
boundaries Sj n and S ou t, the situation that interests us. But we can make our 
situation fit the formulation of the theorem by digging a narrow "channel" from 
Sout to Si n , as depicted in Fig. 2.3. The union of S- m , S ou t, and the channel is a 
closed surface, and integrating b over both sides of the channel produces a zero 
result because da points in opposite directions. We therefore have 

V • bd 3 x = b - da+ b ■ da, 

and noting the opposite orientations of da and n on Si n (see Fig. 2.3), we write this 
as 

V-b<f 3 .r = / b-nda — J b ■ nda 
b ■ hda. 

It is important to notice that the unit vector n points out of both Sj n and S ou t ■ 
We now choose the vector b to be 
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where <fi(x) and ip(x) are two arbitrary functions. We have V • b = <f>W 2 tp — ipV 2 <f) 
and Gauss' theorem gives 

[ (ct)V 2 iP-ipV 2 (t>)d 3 x= j (4>Vip - ipV(t>) ■ nda. (2.4.1) 

JV JS ou t — S\ n 

This is Green's identity. For convenience we write it in the equivalent form 

f (4>V' 2 ip - tpV /2 <j)) d 3 x' = f (0W - VV» • hda', 

J V J S ou t — s\ n 

where all fields are now expressed in terms of the new variables x' . 
We now pick 

4>(x') = $(2/) = scalar potential 

and 

4>(x') = Gb(x,x') = Dirichlet Green's function, 

where the properties of the "Dirichlet Green's function" will be described in detail 
below. For now we assume that it satisfies the equation 

V ,2 G d (x,x') = -4ttS(x-x'), (2.4.2) 

while the scalar potential satisfies Poisson's equation, V' 2 $(a;') = —p(x')/eo. Mak- 
ing these substitutions in the equation following Eq. (2.4.1) leads to 

-47r$(a;) + - f G D (x,x')p(x')d 3 x' = [ h(x')V'G D (x,x') 

£ JV J S out — Sin L 

- G D (x,x')V^(x')} ■ nda', 



or 



$(cc) = — — / G D (x,x')p(x')d 3 x' 
47re Jv 

--?-/" \l>(x')V'G B (x,x')-G D (x,x')V'<l>(x') 

47r Js out -s ia L 



• h da . 



Apart from a remaining simplification, this is the kind of expression we were seeking. 
The volume integral takes care of the charge distribution within V, and the two 
surface integrals take care of the boundary conditions. But there is one problem: 
While the value of the scalar potential is specified on the boundaries, so that the 
functions $(cc') are known inside the surface integrals, we are given no information 
about n ■ V' ( E > (a; / ), its normal derivative, which also appears within the surface 
integrals. It is not clear, therefore, how we might go about evaluating these integrals: 
the problem does not seem to be well posed. 

To sidestep this problem we demand that the Dirichlet Green's function satisfy 
the condition 

Gr>{x, x') = when x' is on the boundaries. (2.4.3) 

Then the boundary terms involving n ■ V' ( E ) (a; / ) simply go away, because they are 
multiplied by Gr,(x 7 x') which vanishes on the boundaries. With this property we 
arrive at 

= — *— / G D (x,x')p(x')d 3 x' 
47re J v 

<P(x')h-V'G D (x,x')da'. (2.4.4) 

47T / q ^ _ q. 

■' Omit Om 
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This is our final expression for the scalar potential in the region V. The problem 
of finding this potential has been reduced to that of finding a Dirichlet Green's 
function that satisfies Eqs. (2.4.2) and (2.4.3). Notice that this is a specialized form 
of our original problem: Gu(x, x') is the potential produced by a point charge of 
strength Aire^ located at x, and its boundary values on Si n and S ou t arc zero. If we 
can solve this problem and obtain the Dirichlet Green's function, then Eq. (2.4.4) 
gives us the means to calculate &(x) in very general circumstances. 

We already know that Gd(x,x') satisfies Eq. (2.4.2) and vanishes when x' lies 
on Si n and S out . We now use Green's identity to show that the Dirichlet Green's 
function is symmetric in its arguments: 

G u (x',x) = G u (x,x'). (2.4.5) 

This implies that it is also a solution to the standard equation satisfied by a Green's 
function, 

V 2 G D (a;, x') = -4tt5(x - x'), (2.4.6) 

which is identical to Eq. (1.8.2). To establish Eq. (2.4.5) we take y to be the 
integration variables in Green's identity, 

f {<pV 2 y 4> - d 3 y = { (^V y V - ^V„0) • nda y , 

JV "'Sout -Sin 

and we pick <j>{y) = Gr>(x, y) and ip(y) = Gu(x', y). Then according to Eq. (2.4.2) 
we have V 2 ,^ = —Att5(x — y) and V 2 ^ = —An5(x' — y). The left-hand side of 
Green's identity gives — 4tt[Gd(x, x') — Gd(x',x)] after integration over d 3 y. The 
right-hand side, on the other hand, gives zero because both </> and are zero when 
y is on Si n and S ou t- This produces Eq. (2.4.5). 

To summarize, we have found that the Dirichlet Green's function satisfies Eq. (2.4.6), 

V 2 G D (a:, x') = -4ttS(x - x'), 

possesses the symmetry property of Eq. (2.4.5), 

Gt>(x',x) = Gb(x,x'), 

and satisfies the boundary conditions of Eq. (2.4.3), 

Gb{x, x') = when x or x' is on the boundaries. 

Notice that the boundary conditions have been generalized in accordance to Eq. (2.4.5): 
the Dirichlet Green's function vanishes whenever x or x' happens to lie on the 
boundaries. Once the Dirichlet Green's function has been found, the potential 
produced by a charge distribution p{x) in a region V between two boundaries S- m 
and Sout on which $(cc) is specified can be obtained by evaluating the integrals 
of Eq. (2.4.4). The task of finding a suitable Dirichlet Green's function can only 
be completed once the boundaries S- ln and S out are fully specified; the form of the 
Green's function will depend on the shapes and locations of these boundaries. In 
the absence of boundaries, Gd(x,x') reduces to the Green's function constructed 
in Sec. 1.8, G(x,x') = \x-x'\-\ 

In the following sections we will consider a restricted class of boundary-value 
problems, for which 

• there is only an inner boundary S ln (S ou t is pushed to infinity and does not 
need to be considered); 

• the inner boundary is spherical (Si„ is described by the statement \x'\ — R 1 
where R is the surface's radius); 
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• the inner boundary is the surface of a grounded conductor ($ vanishes on 

Sin). 

With these restrictions (introduced only for simplicity), Eq. (2.4.4) reduces to 

<P(x) = - i —[ G D (x,a;>(a;')(iV. (2.4.7) 
47r£ o J\ X '\>R 

In these situations, the electric field vanishes inside the conductor and there is an 
induced distribution of charge on the surface. The surface charge density is given 
by Eq. (1.6.9), 

a = e E-h, (2.4.8) 

where E is the electric field just above 5; n . This is given by E = — V$(|cc| = R). To 
proceed we will need to find a concrete expression for the Dirichlet Green's function. 
For the spherical boundary considered here, this is done by solving Eq. (2.4.6) in 
spherical coordinates. 



2.5 Laplace's equation in spherical coordinates 

[The material presented in this section is also covered in Sees. 3.1, 3.2, and 3.5 of 
Jackson's text.] 

When formulated in spherical coordinates (r, 9, </>), Laplace's equation takes the 
form 

2 1 8 ( 2 8ip\ ( 1 8 f^ Q dtp\ , 1 <9 2 V> 



V>=^^ +^— — sing^ + , a- — § = 0. (2.5.1) 

r 2 8r\ or J r 2 sin 9 89 \ 89 J r 2 sin 9 8cj> 2 

To solve this equation we separate the variables according to 

i){r,6,4>) = R{r)Y{6,4>), (2.5.2) 

and we obtain the decoupled equations 



d f 2 dR \ 



r 



dr\ dr ) 

and 



'+1)R (2.5.3) 



8 ( 8Y\ 1 8 2 Y 

1 sm6"-) +^-"-^ =-£(£+ l)Y, (2.5.4) 



sin 6»<96>V d9 J sin 2 9 8 '</> 2 

where 1(1 +1) is a separation constant. 

The solutions to the angular equation (2.5.4) are the spherical harmonics Yt m (9, (f>] 
which are labeled by the two integers I and m. While £ ranges from zero to infinity, 
m is limited to the values —£, —£+ 1, ■■■,£— 1,£. For m > 0, 



Y im {e,4>) - fr^ pr(cos0)e "" ' (2 - 5 - 5) 

where Pp(cos9) are the associated Legendre polynomials. For m < we use 

Y tm = (-l) m Y e *_ m , (2.5.6) 

where an asterisk indicates complex conjugation. Particular spherical harmonics 
are 

1 

Ynn = 



47T 
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Y w 
Y 22 
Y 2 i 
Y20 



sin ye" 



4tt 
1 fib 
4V 2^ 
fl5 
V 8tt 
~5 



cos 0, 



1 



2 V 4tt 



sin 2 0e 2i * 



sin6»cos6>e # , 
(3cos 2 6>- 1). 



Two fundamental properties of the spherical harmonics are that they are orthonor- 
mal functions, and that they form a complete set of functions of the angles 9 and 
4>. The statement of orthonormality is 



/ 



V* 
1 01 m l 



)Y,, 



)dCl = ■ 



(2.5.7) 



where dSl = sin 9 d9d(f> is an element of solid angle. The statement of completeness 
is that any function f(9, (f>) can be represented as a sum over spherical harmonics: 



(2.5.8) 



e=o m=-t 

for some coefficients ft m . By virtue of Eq. (2.5.7), these can in fact be calculated 



as 



(2.5.9) 



Equation (2.5.8) means that the spherical harmonics form a complete set of basis 
functions on the sphere. 

It is interesting to see what happens when Eq. (2.5.9) is substituted into Eq. (2.5.8). 
To avoid confusion we change the variables of integration to 9' and <f)': 

= EE J f(o',<P')Y; m (9',4>')dn' 



-I 



/(#',</>') 



EE r ^'>')w^) 



e m 



dfl'. 



The quantity within the large square brackets is such that when it is multiplied 
by f(6',<p') and integrated over the primed angles, it returns f(9,<ft). This must 
therefore be a product of two 5-functions, one for 9 and the other for <fi. More 
precisely stated, 



E E Y LW,<t>')Ye m {0,4>) = 



5(6 - 6')5{(j> - <f/) 



i=0 m=- 



sm( 



(2.5.10) 



where the factor of 1/ sin 9 was inserted to compensate for the factor of sin 9' in 
dQ' (the (5-function is enforcing the condition 9' = 9). Equation (2.5.10) is known 
as the completeness relation for the spherical harmonics. This is analogous to a 
well-known identity, 
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in which the integral over dk replaces the discrete summation over i and m; the 
basis functions [2i:)~ 1 /' 2 e lkx are then analogous to the spherical harmonics. 

The solutions to the radial equation (2.5.3) are power laws, R oc r e or R oc 
r -(^+!)_ The general solution is 

R im (r) = a lr y + b im r-^ +1 \ (2.5.11) 

where ae m and bi m are constants. 

The general solution to Laplace's equation is obtained by combining Eqs. (2.5.2), 
(2.5.5), (2.5.11) and summing over all possible values of I and to: 

oo I 

iP(r,6,<f>) = J2 E (a lm r l + b lm r-^y im {e,4>)- (2.5.12) 

£=Q m =-i 

The coefficients ae m and be m are determined by the boundary conditions imposed 
on ip(r, 6*, (/)). If, for example, ip is to be well behaved at r — 0, then bi m = 0. If, on 
the other hand, tp is to vanish at r = oo, then a em = 0. From these observations 
we deduce that the only solution to Laplace's equation that is well behaved at the 
origin and vanishes at infinity is the trivial solution -0 = 0. 



2.6 Green's function in the absence of boundaries 

[The material presented in this section is also covered in Sees. 3.6, and 3.9 of Jack- 
son's text.] 

Keeping in mind that our ultimate goal is to obtain the Dirichlct Green's function 
Gu(x,x') for a spherical inner boundary, in this section we attempt something 
simpler and solve 

V 2 G(x, x') = -4n6(x - x') (2.6.1) 

in the absence of boundaries. We already know the answer: In Sec. 1.8 we found 
that 

• (2-6.2) 

We will obtain an alternative expression for this, one that can be generalized to 
account for the presence of spherical boundaries. 

The 5- function on the right-hand side of Eq. (2.6.1) can be represented as 

= 5{r-r')8{6-e')5^ 
r 2 sin 9 

S(r - r') 



J2Ye* m (0'A')Ye m (0,<P), (2.6.3) 



fm 



where we have involved Eq. (2.5.10). The factor of l/(r 2 sin#) was inserted be- 
cause in spherical coordinates, the volume element is d 3 x — r 2 sin 9 drd9d(j), and the 
Jacobian factor of r 2 sin# must be compensated for. We then have, for example, 

f{x') = J fix)5ix-x')d 3 x 

= [fir M ) 5{r 7' )5{e - e ' )5{ t^\^^9drd9d ( p 
J r z sin 9 

= J fir, 9, 4>)8ir - r')<5(6» - 0')<J(<£ - 0') drd9d(j) 
= /(r',0',A 
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as we should. 

Inspired by Eq. (2.6.3) we expand the Green's function as 



g(x,x') = J29Ury)Y; m (e , ,<p , )Y em (e,^, 



(2.6.4) 



using the spherical harmonics Y^ m (9, <f>) for angular functions and gt m (r, r')Y^ m (9' , <j>') 
for radial functions. For convenience we have made the "radial function" depend 
on x' through the variables r', 9', and 0'; these are treated as constant parameters 
when substituting G(x,x') into Eq. (2.6.1). Applying the Laplacian operator of 
Eq. (2.5.1) on the Green's function yields 



V 2 G(x,x')=J2 



LA 

r 2 dr 



! dge„ 
dr 



£{£ + !) 



~9tn 



Y; m (6',4>')Y em (0,^, 



and setting this equal to — AnS(x — x') as expressed in Eq. (2.6.3) produces an 
ordinary differential equation for ge m {r, r'): 



dr 



> dge„ 
dr 



-£(£ + l)g. 



im 



-47T<5(r — r'). 



(2.6.5) 



This equation implies that the radial function depends only on £, and not on m: 
gi m (r,r') = gz(r,r'). To simplify the notation we will now omit the label "£" on 
the radial function. 

To solve Eq. (2.6.5) we first observe that when r ^ r', the right-hand side of the 
equation is zero and g is a solution to Eq. (2.5.3). We let #<(r, r') be the solution to 
the left of r' (for r < r') and g>(r,r') be the solution to the right of r' (for r > r'). 
We then assume that g{r 7 r') can be obtained by joining the solutions at r = r'\ 



g(r, r') = g< (r, r')9(r' - r) + g> (r, r')9(r ■ 



(2.6.6) 



where 9(r' — r) and 9(r — r') are step functions. Differentiating this equation with 
respect to r gives 

g'(r,r') = g' < (r,r')9(r' - r) + g' > (r,r')9(r - r') + [g>(r',r') - g < (r' ,r')]S(r - r'), 

where a prime on g indicates differentiation. We have involved the distributional 
identities 9'(r — r') = 5(r — r'), 9'(r' — r) = —5(r — r'), and f(r)S(r — r') — 
f(r')5(r — r'); these were first encountered in Sec. 1.6. Multiplying by r 2 yields 

r 2 g\r,r') = r 2 g' < (r,r')9(r' - r) + r 2 g' > (r,r')9(r - r') 
+ r' 2 [g > (r',r')-g < (r\r')]S(r-r') 

and taking another derivative brings 

[r 2 g'(r,r')]' = [r'g'^r, r')]' 6{r' - r) + [r'g'^r, r')]' 6{r - r 1 ) 
+ r' 2 [g' > (r',r')-g' < (r' 1 r')]S(r-r') 
+ r' 2 [g > (r',r')-g < (r>,r>)]6'(r-r'). 

Substituting this and Eq. (2.6.6) back into Eq. (2.6.5), we obtain 



-47T<5(r - 



A 

dr 
+ 

+ r 
+ r 



>.dg< 



dr 



dr \ 

12 



( 2 d92 



dr 



£(£ + l) 9< 
-£(£ + l) 9> 



dg> i i i \ dg < i , 
(r ,r ) - — (r ,r) 



dr 



dr 



9{r' - r) 
9{r - r') 
5{r-r') 



9>{r',r') - g<(r',r') —5{r-r'). 
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For the right-hand side to match the left-hand side we must demand that g<(r, r') 
and g>(r, r') be solutions to the homogeneous version of Eq. (2.6.5), as was already 
understood, and that these solutions be matched according to 

9>(r',r')-g < (r',r , )=0 (2.6.7) 

and 

^</y)-^</y) = -^ (2.6. 8 ) 

The function g(r,r') is therefore continuous at r = r' , but its first derivative must 
be discontinuous to account for the (^-function on the right-hand side of Eq. (2.6.5). 

Because g< is a solution to Eq. (2.5.3), according to Eq. (2.5.11) it must be a 
linear combination of terms proportional to r E and l/r l+1 . But only the first term 
is well behaved at r = 0, and we set 

9<(r,r') = a(r')r e . 

For g> we select instead 

g> (r,r')=b(r')/r e+ \ 

which is well behaved at r = oo. It does not matter that <?< diverges at r = oo 
because this function is restricted to the domain r < r' < oo; similarly, it does not 
matter that <?> diverges at r = because it is restricted to the domain r > r' > 0. To 
determine the "constants" a and b we use the matching conditions. First, Eq. (2.6.7) 
gives b = ar' 2e+1 . Second, Eq. (2.6.8) gives a = 4irr'- ( -< ! + 1 '> /(2£ + 1). So 

9<M = ^^i (2-6.9) 

and 

We observe a nice symmetry between these two results. Combining Eqs. (2.6.6), 
(2.6.9), and (2.6.10), we find that the radial function can be expressed in the com- 
pact form 

^y) = ^-|r, (2.6.H) 

where r< denotes the lesser of r and r', while r> is the greater of r and r'. Thus, 
if r < r' then r< = r, r> = r', and Eq. (2.6.11) gives g(r,r') = g<(r,r'), as it 
should. Similarly, when r > r' we have r< = r', r> = r, and Eq. (2.6.11) gives 
g(r,r') = g>(r,r'). 

Substituting Eq. (2.6.11) into Eq. (2.6.4) and taking into account Eq. (2.6.2), 
we arrive at 

=E E ^^(^^(M). (2-6.12) 

1 1 e=0 m=-i ' > 

This very useful identity is known as the addition theorem for spherical harmonics. 

At first sight the usefulness of this identity might seem doubtful. After all, 
we have turned something simple like an inverse distance into something horrible 
involving special functions and an infinite double sum. But we will see that when 
l/\x — x'\ appears inside an integral, the representation of Eq. (2.6.12) can be 
very useful indeed. Not relying on this identity would mean having to evaluate an 
integral that involves 



\x-x'\ ^r 2 - 2rr'[sin 6 sin 6' cos(</> -(/>')+ cos 9 cos ¥] + r' 2 ' 
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and such integrals typically cannot be evaluated directly: they are just too compli- 
cated. 

Let us consider an example to illustrate the power of the addition theorem. We 
want to calculate the electrostatic potential both inside and outside a spherical 
distribution of charge with density p(x') = p(r'). There are no boundaries in this 
problem, but the distribution is confined to a sphere of radius R. After substituting 
Eq. (2.6.12) into Eq. (2.1.3) and writing d 3 x' = r' 2 dr'dfl', we obtain 



= 4^ E Yem{6 > ® I ^mP {r '> n dr ' J Y ^(°'> #) dn '- 

To evaluate the angular integral we multiply and divide by (4-7r) -1 / 2 = Y ao (9', 4>'): 

J Y; m (e',<p')dn' = V^J Y; m (e',<p')Y oo (0',cp')dn' = V^s t , s m , , 

where we have involved Eq. (2.5.7). At this stage we happily realize that the infinite 
sum over £ and m involves but a single term. We therefore have 

1 f R r' 2 1 f R r' 2 

&(x) = - 4ttV4^Yoo(9,(I>) / —p(r')dr' = - 4tt / —p{r')dr', 

47re Jo r > 47r£ o Jo r > 

and as we might have expected from the spherical symmetry of the problem, the 
potential depends on r only. 

To evaluate the radial integral we must be careful with the meaning of r> = 
max(r,r'). Suppose first that x is outside the charge distribution, so that r > R. 
Then r > r' , r> = r, and we have 

4tt / —p(r')dr' = -4tt / p(r')r' 2 dr' = ^, 
Jo r > r Jo r 

where Q = J p(x') d 3 x' = 4ir J a R p(r')r' 2 dr' is the total charge of the distribution. 
In this case we have 

*outO«0 = T^--. (2-6.13) 
47re r 

the same result as in Eq. (2.2.2). Suppose next that x is inside the charge distri- 
bution, so that r < R. Then part of the integration covers the interval < r' < r 
and the remaining part covers r < r' < R. In the first part r' is smaller than r and 
r> = r; in the second part r' is larger than r and r> = r'. So 

f R r' 2 1 r r R 

4?r / — p(r')dr' = -4?r / p(r')r' 2 dr' + 4tt / p(r')r' dr' . 

Jo r > r Jo Jr 

In the first term we recognize 

q(r)= [ p(x') d 3 x' = 4tt f p(r')r' 2 dr', (2.6.14) 

J\x'\<r Jo 

the charge enclosed by a sphere of radius r; if r > R then q(r) = Q. The internal 
potential can thus be expressed as 



47re 



^1+4. 
r 



/ p(r')r' dr' 

J r 



(2.6.15) 



While the expression (2.6.15) for the internal potential depends on the details 
of the charge distribution — the particular form of the function p(r') — a nicer 
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conclusion applies to the electric field. The field is obtained by differentiating $; 
with respect to r, and 



d<Z> m 
dr 



1 

47re 



q_ Idq 
r 2 r dr 



4-rrpr 



But since dq/dr = iirpr 2 according to Eq. (2.6.14), we see that the last two terms 
cancel out. The electric field is then 



E(x) = 



1 



-r, 



47reo r 



(2.6.16) 



both inside and outside the charge distribution; recall that when r > R, q{r) be- 
comes Q, the total charge. We have recovered the well-known result that for any 
spherical distribution of charge, the electric field at a radius r depends only on the 
charge q(r) enclosed by a sphere of that radius; and the field is the same as if all 
this charge were concentrated in a point at r = 0. 



2.7 Dirichlet Green's function for a spherical 

inner boundary 

[The material presented in this section is also covered in Sec. 3.9 of Jackson's text.] 
We now return to the task of finding the Dirichlet Green's function Gd(x, x') for 
the specific situation described near the end of Sec. 2.4, in which we have a spherical 
inner boundary at \x\ = R and no outer boundary. It will be a simple matter to 
modify the treatment presented in Sec. 2.6 to account for the inner boundary. 
As in Eq. (2.6.4), the Dirichlet Green's function can be expanded as 

G D (x,x') =Y,9t(rJ)Yt m (#,<l>')Y lm {e,4>)- (2-7.1) 

As in Eq. (2.6.6), the radial function can be expressed as 

9t(r, r') = 9< (r, r')0(r' - r) + 9> (r, r')6{r - r'), (2.7.2) 

where g< and g> are solutions to Eq. (2.5.3) that must also satisfy the matching 
conditions of Eqs. (2.6.7) and (2.6.8), 

9>(r',r')-g < (r',r , )=0 (2.7.3) 

and 



-(r>y)-^(r',r>) = -^. (2.7 .4) 



dg> , , dg< , , 4tt 

^r (r ' r) ^ (r ' r) = -^ 

The only difference with respect to what was done in Sec. 2.6 is that now, according 
to Eq. (2.4.3), the Green's function must vanish on the inner boundary. This gives 
the additional conditions g e (r = R, r') = g e (r, r' = R) = 0, or 

9< (r = R,r') = 0, 9> (r, r' = R) = 0; (2.7.5) 

recall that R is now the smallest possible value of both r and r'. 

Since g<(r, r') is no longer required to be well behaved at r = 0, we must set it 
equal to 

g< ( r , r') = a(r', R)r l + a{r', R)/r l+ \ 
the general solution to Eq. (2.5.3). For g>(r,r') we choose instead 

g> (r,r') = b(r',R)/r e +\ 
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which is well behaved at r = oo. Solving Eqs. (2.7.3)-(2.7.5) for the "constants" 
produces a = 4 7 rr / -^+ 1 )/(2f + 1), a = -aR 2e +\ and b = a{r' 2e+1 - R 2l+1 ). Col- 
lecting these results gives 

9<M = ^ ri ^-%-) (2-7.6) 

and 

Combining Eqs. (2.7.2), (2.7.6), and (2.7.7), we find that the radial function can be 
expressed in the compact form 

, m ^ 1 / / R 2i+1 \ ^ rl i+1 -R 2i+1 , nno , 



2^ + 1 ri +1 V < / 2^ + 1 (rr') £+1 

where r< denotes the lesser of r and r', while r> is the greater of r and r' (notice 
that r<r> = rr'). We see that ge(r,r') vanishes on the inner boundary: When 
r = R it must be that r < r', so that r< = r = i?, and the terms within the 
brackets vanish; on the other hand, when r' = R it must be that r' < r, so that 
r< = r' = _R, and the same conclusion applies. 

The Dirichlet Green's function is obtained by substituting Eq. (2.7.8) into Eq. (2.7.1): 

A 1 / R 2i+1 \ 
G D (x, X ') = J2^J^Ti{ r < - 7 +r)Y? m (0 , ,<t>')Y tm (6,<i>). (2.7.9) 

This is slightly more complicated than Eq. (2.6.12), but still manageable. Notice 
that the Dirichlet Green's function reduces to l/\x — x'\ when R = 0. It would be 
straightforward to generalize this derivation to account also for the presence of an 
outer boundary; we shall not pursue this here. 



2.8 Point charge outside a grounded, spherical 

conductor 

[The material presented in this section is also covered in Sec. 2.2 of Jackson's text.] 
The electrostatic potential outside a grounded, spherical conductor of radius R 
is given by Eq. (2.4.7), 

$(x) = — *— / G B {x,x')p(x')d 3 x', (2.8.1) 

47re J\x'\>R 

and the Dirichlet Green's function was obtained as 

a 1 / p2£+l \ 

Gn(x,x') =J2 27TT^r< ~ - 7 +r) Y i m (0'^')Y em (0A) (2.8.2) 

in the preceding section. The charge density induced on the surface of the conductor 
is given by Eq. (2.4.8), a = e$E ■ n, where E is the electric field just above the 
conductor and n = r is the surface's unit normal. Because E = — V$, this works 
out to be 

a = -e - 



dr 



(2.8.3) 

r=R 



We place a point charge q on the z axis, at a distance zq > R from the centre of 
the conducting sphere. The charge's position vector is 2 = z z and its spherical 
coordinates are r = z , 9 n = 0, and 0o is undetermined. The charge density is 

p(x') = qS(x' - z ). (2.8.4) 
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Substituting this into Eq. (2.8.1) gives 



*(se) 



47re 



47re 



E 



47T 



2^+1 Vrl +1 



R 2£+1 

e+i e+i j^mtMo)^™ (#,</>), 
r < r > 



where r< now stands for the lesser of r and z , while r> stands for the greater of 
the two. It is a property of the spherical harmonics that 



Y em (0,<h) 



2£+l 

47T 



J m,0, 



so that the potential is actually independent of <j)Q and involves a single sum over 
I. It is another property of the spherical harmonics that for m = 0, they reduce to 
ordinary Legendre polynomials, 



y*,o(M) = 



2^ + 1 

47T 



P e (cos9). 



After cleaning up the algebra, the final expression for the potential is 



47re ^ 



£=0 v ' > 



^+1 



P e (cos9). 



(2.8.5) 



We see that the potential depends on r and 6* only; the fact that it does not depend 
on cf> reflects the axial symmetry of the problem. 

Before moving on we pause and recall the main properties of the Legendre 
polynomials. These are ordinary polynomials of order £• if £ is even the function 
Pi{x) contains only even terms, while if I is odd it contains only odd terms. For 
example, 



Po(x) 


= 1, 




Pi(x) 


= x, 




P2(X) 


- 


1), 


Ps(x) 


- 2^- 


3a;), 


P 4 {x) 


= l(35x 4 - 


- 30x 2 + 3) 



The Legendre polynomials are orthogonal functions, 



P i ,(x)P l (x)dx=-^- l < 



and they possess the special values 
,„(2n-l)!! 



P2n{Q) = (-If 



(2n)\\ 



P 2 n+i(0)=0, P t (l) = l, 



(2.8.6) 

(2.8.7) 
• (2) or 



where the double factorial notation means (2n)H = (2n)(2n — 2)(2n — 4) 
(2n - 1)!! = (2n - l)(2n - 3)(2n - 5) ■ ■ ■ (1). 

The potential of Eq. (2.8.5) is expressed as an infinite sum over i. We will 
evaluate this sum in a moment, but for now we calculate the charge density a (9) 
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Figure 2.4: Surface charge induced by a point charge above a spherical conductor. 
Plotted is ATrR 2 a/q' as a function of 9 for zq/R = 1.3 (upper curve), zq/R = 1.6 
(middle curve), and z /R = 1.9 (lower curve). The charge density is peaked at 
9 = and its width increases with increasing values of zq /R: the closer the charge 
is to the surface the narrower is the charge distribution. 

on the surface of the conductor. For r close to R we have r < zq, r < = r, r> = zq, 
and Eq. (2.8.5) becomes 

Let /(r) = r e — H™* 1 /r +1 be the function of r that appears inside the brackets. 
Its derivative is f'(r) = Ir 1 ^ 1 + (£ + l)R 2t+1 /r e+2 and evaluating this at r = R 
gives f'(R) = (2£+l)R e - 1 . From Eq. (2.8.3) we then obtain 

£ Z 

or 

j oo 

77 1=0 

where 

q' = -q(R/z ) (2.8.9) 

will be shown presently to be the total induced charge on the surface of the con- 
ductor. Plots of the charge density are presented in Fig. 2.4 for selected values of 
R/z Q . 

That q' is truly the total surface charge can be shown by evaluating the integral 
of the charge density over the surface of the conductor; we should recover q' = 
§ a da. To evaluate this we let da = R 2 dQ = R 2 sin 9 d6d(f>, which reduces to 
da = 2irR 2 sin 9 d9 after integration over <f>. Then 

ada= -q'J2{2£ + l)(R/z o y / P^cos 0) sin 9 d9. 

2 „_n JO 
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The integral over 6 is equal to J j Pe(x)P a (x) dx — 2£^o by virtue of Eq. (2.8.6), 
and we obtain § a da = q' ', as required. 

The potential of Eq. (2.8.5) can be decomposed as 

$(a:) = $ g (aj) + $,/(aj), (2.8.10) 

where 

OC £ 

e=o > 

is the potential that would be obtained in the absence of a boundary, and 

^)-^j:%Pdcos0) (2.8.12) 

is the modification introduced by the boundary condition $(r = R) = 0; here 
q' = —q(R/zo) was defined in Eq. (2.8.9) and 

4 = R 2 /z . (2.8.13) 

Notice that z' /R = R/z < 1, so that z' < R<r. 

The sums over I can now be evaluated. For $ 9 we already know that Eq. (2.8.11) 
is equivalent to 

$,(*)= 1 , q , , (2.8.14) 
9V 47re |aj-z | 

because $ g is the potential produced by a charge g at a position z = z i in the 
absence of boundaries. By analogy, because z' Q < r we have that Eq. (2.8.12) is 
equivalent to 

j i_ 

47TC0 \X — Zq| 

where z' = z' z. This is the potential of a fictitious image charge q' at a position z' 
inside the conductor. While in reality the charge q' is distributed on the surface of 
the conductor, the potential <ly is the same as if all this charge were concentrated 

This observation is the basis for the method of images, which consists of finding 
solutions for the potential in the presence of conducting surfaces by guessing the 
correct positions of the correct number of image charges. It is a powerful, but tricky, 
method. 



M*)= ^/. . ( 2 - 8 - 15 ) 



2.9 Ring of charge outside a grounded, spherical 

conductor 

In this section we consider a uniformly charged ring of radius a > R. If the ring is 
placed in the x-y plane, its charge density is described by 

p(x') = ^_5(r'-a)5(6'-%), (2.9.1) 

where q = J p(x') d 3 x' is the ring's total charge. The 5-function in r' indicates that 
the ring is infinitely thin, while the (5-function in 6' tells us that it is placed in the 
x-y plane, as was previously stated. 

Substituting Eq. (2.9.1) into Eq. (2.8.1) along with Eq. (2.8.2) gives 
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where d£l' = sin 6' d6'd<f>, r< = min(r, r'), and r> = max(r, r'). To perform the <j>' 
integration we recall that Y£ m oc e~ lm ^ , so that integrating gives zero unless m = 0. 
So 

f -27T 

Y; m (9',<f>')d<t>' = 2nY tl0 (9')5 m , . 



I 

Jo 



la 

Integration over 9' then gives 



21 + 1 

2^,0(1)^,0 = 2ir ] J^—P e (0)5 m , . 
On the other hand, the radial integration gives 

a 



j+i r e+i r e+i 



where r< now stands for the lesser of r and a, while r> stands for the greater of 
the two. With these results the potential becomes 

•<-)-^i:^o(|r-^r)^«(0). 

Invoking once more the relation between Y«,o(0) and P^(cos#), we arrive at 

= 4^x>K# _ ) p<(cosg) - (2 - 9 - 2) 

Once more the potential is independent of 0, and is expressed as an infinite sum 
over I. Once again we will see that the sum over I can be evaluated. 

We first calculate the charge density on the surface of the conductor. For r close 
to R and therefore smaller than a, we have r< = r, r> = a, and Eq. (2.9.2) becomes 

^) = 4^oE^(°)^(cos^), 

where the function f(r) = r E — R 2e - +1 /r* +1 was previously encountered in Sec. 2.8. 
We recall that f'(R) = {21 + l)i^"\ and Eq. (2.8.3) gives 



or 



where 



= 4^R2 ^2(2£ + l)P t {O){R/aYPt(coB0), (2.9.3) 

£=0 

q' = -q(R/a) (2.9.4) 

is the total charge on the conducting surface. Plots of the surface charge density 
are presented in Fig. 2.5 for selected values of R/a. 

As in the previous section the potential of Eq. (2.9.2) can be decomposed as 

$(a:) = $ g (aj) + $,/(a:), (2.9.5) 

where <& q is the potential of a ring of total charge q and radius a, while <Ey is the 
potential of a fictitious image ring of total charge q' and radius a! = R 2 /a. (You 
should work through the details and make sure that these statements are true.) 
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Figure 2.5: Surface charge induced by a ring of charge outside a spherical conduc- 
tor. Plotted is 4nR 2 a/q' as a function of 9 for a/R — 1.3 (upper curve), a/R = 1.6 
(middle curve), and a/R= 1.9 (lower curve). The charge density is peaked at 9 = | 
and its width increases with increasing values of a/i?: the closer the ring is to the 
surface the narrower is the charge distribution. 

2.10 Multipole expansion of the electric field 

[The material presented in this section is also covered in Sec. 4.1 of Jackson's text.] 
We now leave boundary-value problems behind, and derive the fact that in the 
absence of boundaries, the electrostatic potential outside a bounded distribution of 
charge can be characterized by a (potentially infinite) number of multipole moments. 
Among them are 

Q = J p(x) d 3 x = total charge = monopole moment, (2.10.1) 

p a = J p{x)x a d 3 x = dipole moment vector, (2.10.2) 

Q a b = J p(x) (3x a xt, — r 2 5 a b) d 3 x = quadrupole moment tensor. (2. 10.3) 

Additional multipole moments can be constructed in a similar way, up to an infi- 
nite number of tensorial indices. But to proceed along this road quickly becomes 
cumbersome, and we shall find a better way of packaging the components of these 
multipole-moment tensors. Notice that the quadrupole moment is a symmetric and 
tracefree tensor: Q oa — Q a b and Q aa — (summation over the repeated index is 
understood); it therefore possesses 5 independent components. 
Consider the complex quantities 



qtm = I p(xyY e * m (9,(t))d 3 x. 

By virtue of Eq. (2.5.6), they satisfy the identity 

qt,- m = (-l) m q$m- 



(2.10.4) 



(2.10.5) 
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Because m ranges from — £ to £, for a given £ there are 2£ + 1 independent real 
quantities within the qe m s>- For £ = we have only one, and we will see that goo is 
proportional to the total charge. For I — 1 we have three independent quantities, 
and we will see that the gi m s give the components of the dipolc moment vector. For 
£ = 2 we have five independent quantities, and we will see that the (femS give the 
components of the quadrupole moment tensor. Thus, the quantities of Eq. (2.10.4) 
give a convenient packaging of the multipole moments, to all orders. 

To evaluate Eq. (2.10.4) for I = we recall that Y* = {Air)- 1 / 2 , so that 

9oo = -!=Q. (2.10.6) 

V47T 



We see that indeed, goo is directly related to the total charge. For I = 1 we have 
rift = -(3/87r) 1 / 2 (rsin6»)e-^ = -(3/8 7 r) 1 /2( x _i y ) an d r y* = (3/4^)^^ cos = 
(3/4tt) 1 / 2 z. This gives 

/ 3 / 3 

1u = -d^{p x -iPy), 1w = \I^Pz 1 (2.10.7) 

and we see that indeed, the <?i r „s are directly related to the dipolc moment vec- 
tor. For I = 2 we have r 2 Y 2 * 2 = {15/32n) 1 ^(r sin Oe-^f = (15/32w) 1 / 2 (x - iy) 2 , 
r 2 Y 2 \ = -(15/87r) 1 / 2 (rsin6»e"^)(rcos6') = -{lb/Sir) 1 ' 2 ^ - iy)z, and r 2 Y 2 \ = 
(5/167r) 1 / 2 r 2 (3cos 2 6'- 1) = {5/16tt)^ 2 (3z 2 _ r 2y Thig gives 

1 / 15 

922 = J^y-^(Qxx - Qyy ~ 2^0^), 
1 / 15 

921 = -qJ ^(Qxz-iQyz), 

920 = l\f^Q^ ( 2 - 10 - 8 ) 

and we see that indeed, the g2mS are directly related to the quadrupole moment 
tensor. For higher values of £ it is best to stick with the definition of Eq. (2.10.4). 

How are the multipole moments qg m related to the potential? The answer comes 
from Eq. (2.1.3) in which we set $o( £C ) = (because we are no longer considering 
boundaries) and substitute Eq. (2.6.12): 

Because we are looking at the potential outside the charge distribution, r is larger 
than r' and we can set r < —r',r > — r. This gives 

lm 

We recognize the integral over the primed variables as the definition of the <# m s, 
and we arrive at 

u i=Q m =-e 



This is an expansion of the potential in inverse powers of r, and the multipole 
moments play the role of expansion coefficients. 
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The leading term in Eq. (2.10.9) comes from the monopole moment qooj or total 
charge Q: 

$monopole(aO = (2.10.10) 

The next term comes from the dipole moment. Using Eq. (2.10.7) and recalling the 
explicit expressions for the spherical harmonics of degree £ = 1, we have 

1 3 

^2 9l,m^l,m = ^ (PxX + PyV + p z z) 
m— — 1 



and 



*di P oie(aO = -. (2.10.11) 



as in Eq. (2.3.5). The next term in the expansion of the electrostatic potential comes 
from the quadrupole moment. Using Eq. (2.10.8) and the explicit expressions for 
the spherical harmonics of degree £ = 2, it is easy to show that 

*quadrupole(*) = J^Q«6^- (2.10.12) 

In a situation in which Q ^ 0, the electrostatic potential far outside the charge 
distribution is well approximated by ^monopole! in this case the potential is charac- 
terized by a single quantity, the total charge Q. When the total charge vanishes, 
however, the potential is approximated by 3>di P oie, and in this case we need three 
quantities, the components of p, the dipole moment vector. When this vanishes 
also, the potential is more complicated and is approximated by $ quadruple! in this 
case the characterization of the potential involves the five independent components 
of the quadrupole moment tensor. 

It is important to be aware of the fact that the multipole moments qp m depend 
on the choice of origin for the coordinate system. For example, a point charge q 
at x = has a charge density described by p(x) = qS(x) and multipole moments 
given by 

_ f q/Vhr £ = 

qim ~ \ otherwise ' 

As expected, only the monopole moment is nonzero. But the same point charge 
placed at another position Xq possesses a very different set of multipole moments. 
In this case the density is p{x) = q5(x — x ) and Eq. (2.10.4) gives 

qtm = qr o e Y e * m (0 o ,(t>o), 

where (r , 0q,4>o) are the spherical coordinates of the point x . We have agreement 
only for <7oo> the lowest nonvanishing multipole moment. This is a special case of 
a general result: Only the lowest nonvanishing multipole moments have values that 
are independent of the choice of origin for the coordinate system; all other moments 
change under a translation of the coordinates. We shall not provide a proof of this 
statement. 



2.11 Multipolar fields 

To conclude this chapter we construct plots of electric field lines that would be 
produced by a pure multipole of order £. We shall consider only axisymmetric 
situations and set m = 0. 



2.11 Multipolar fields 



43 



The potential for such a multipole is obtained from Eq. (2.10.9), in which we 
set [4ir/(2£ + l)]q m Y m {6 , 4>) = q e P e {cos6). This gives 

and the corresponding electric field is 

E{r,6) = iTo\(t+ l)Pe ( cos 0)f + sin 0Pi (cos 0)§], (2.11.2) 

47reo r t+ L J 

in which a prime indicates differentiation with respect to the argument. 

The electric field lines are represented by parametric curves x(X) for which A is 
an arbitrary parameter. These curves are integral curves of the vector field E(x), 
which means that the tangent vector to the curves is defined to be everywhere equal 
to E. Thus, the field lines are determined by solving the differential equations 

r] qr> 

-=£(*). (2.11.3) 

In the case of an axisymmetric field the curves can depend on r and 9 only, and we 
have 

dx dx dr dx dO dr „ dd ~ 

— = 1 = — r + r — 8. 

dX dr d\ 86 dX dX dX 

This shows that dr/dX = E r and rd9/dX = Eg, and using Eq. (2.11.2) produces 
dr oc {£ + 1)P(, and rdO oc sin#P £ ', or 

r ~ sin0P;(cos0) (1 - /i 2 )P;( M ) M ' 1 ' 

where /i = cos 0. The field lines are therefore determined by integrating Eq. (2.11.4). 

To accomplish this it is useful to recall the differential equation satisfied by the 
Legendre polynomials, 

+ « + iH*-a 

We solve this algebraically for Pi and substitute into Eq. (2.11.4) to get 

dr [(l-^)Pj]' 
% (l-^)P^) d ^ 

This can be integrated at once, and we obtain 

^lnr = ln[(l - ^ 2 )P'i\ + constant, 



or 



r(9) = r Q \sm 2 0P^cos6)\ 1/e , (2.11.5) 



where ro is constant on each curve, but varies from one curve to the next. Equation 
(2.11.5) gives the mathematical description of the electric field lines of an axisym- 
metric multipole of degree I. The field lines are plotted for selected values of I in 
Fig. 2.6. 




1 
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2.12 Problems 

1. Calculate the electrostatic potential <f>(x) inside a charge distribution whose 
density is given by 

3 

p{x) = ^Po(r/r ) 2 sin 2 9 for r < r , 

where po and r are constants. For r > r the density is zero. There are no 
boundary surfaces in this problem. 

2. (Jackson's Problem 1.12) Prove Green's reciprocation theorem: If $1 is the 
potential due to a volume charge density p\ within a volume V and a surface 
charge density u\ on the conducting surface S bounding the volume V, while 
$2 is the potential due to another charge distribution p 2 and er 2 , then 



X 



Pi$ 2 d 3 x+ f cr 1 $ 2 ^a= / P2®i d 3 x + f a 2 ^ida. 
Js Jv Js 



A point charge q is brought to a distance zq > R measured from the centre 
of a conducting sphere of radius R (as in Sec. 2.8). The conductor, however, 
is not grounded, but kept in careful isolation; as it initially supports a zero 
net charge on its surface, its final net surface charge must also vanish. The 
surface of the conductor is therefore at a constant potential V ^ 0. Find V 
as a function of q, zq, and R. 

Suppose that at a spherical boundary surface S of radius R, the electrostatic 
potential is assigned the values <fr(R,6,<p) = V sin 9 cos 9 cos <f>, where V is a 
constant. Suppose that there are no charges outside this surface, and that S 
is the only boundary present in the problem. Find the potential for r > R. 

In this problem we consider the electrostatic potential due to a uniformly 
charged disk located outside a grounded, conducting sphere. The conductor 
has a radius R and the potential over its surface is zero. The disk has an inner 
radius a > R and an outer radius b > a. It is convenient to place the disk 
in the equatorial plane of the spherical coordinate system, so that its charge 
density is described by 

3g g(g-f) 
P(X)= 2^^^' ° <r<6 ' 
with q denoting the disk's total charge. 

a) Verify that q is indeed the disk's total charge. 

b) Calculate the potential in the region R < r < a. 

c) Calculate the surface charge density a(9) on the surface of the conducting 

sphere, and provide a plot of this function. (Use the values q = 1, R = 1, 
a = 1.2, and b = 2; make sure to include a sufficient number of terms in 
the infinite sum over I.) 

d) Show that the total charge on the surface of the conductor is given by 

, _ 3q (a + b)R 
q ~ ~Y a 2 +ab + b 2 ' 

We generalize the situation considered in Sec. 2.8 by placing an additional 
charge — q at position — Zq = —zqZ below the conducting sphere. The first 
charge +q is still at +z a = z a z, the conducting sphere still has a radius R, 
and the sphere is still grounded. 

Calculate the vector p, the total electric dipole moment of this charge distri- 
bution. Express your result in terms of q, z a , and R. 
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7. Suppose that the upper hemisphere of a conducting sphere is maintained at 
a potential +V, while its lower hemisphere is maintained at a potential —V. 
(The two hemispheres are separated by a thin layer of insulating material, but 
there is no need to include this into the mathematical model.) The conducting 
surface has a radius R, and there are no charges outside the conductor. 

Show that the potential outside the conducting sphere can be expressed as 

oo 

$(r,0) = Vj2^(R/rY +1 P e (cos9), 
e=i 

where the numerical coefficients are zero if £ is even, and equal to (21 + 
1) Jq 1 Pe(x) dx if £ is odd. Calculate ae explicitly for £ = (1, 3, 5). 

8. A conducting sphere of radius R is maintained at a potential V. It is sur- 
rounded by a thin spherical shell of radius a > R on which there is a surface 
charge density proportional to cos6>. The volume density of charge is given 
by 

3p 

p(r, 9) = - — - cos 9 S(r — a), 

where p is a constant — it is the dipole moment of the charge distribution on 
the thin shell. 

a) Verify that the shell's total charge is zero, and that p = pi is the dipole 

moment vector of this charge distribution. 

b) Calculate the potential in the region R < r < a. 

c) Calculate the potential in the region r > a. 

d) What is the total charge of this configuration? What is the total dipole 

moment vector? (The total charge configuration includes the free charge 
on the spherical shell and the induced charge on the surface of the con- 
ductor.) 

e) Calculate the density of charge a (6) on the surface of the conductor. What 

is the total charge on the conductor? 

9. A dielectric sphere of radius R (the inner boundary) is maintained at a po- 
tential $ = Vb cos 9. It is surrounded by a thin spherical shell of radius a > R 
with a density of charge given by 

p = ^ cos9S(r — a), 

where p is a constant. 

Find the value of Vq which makes the scalar potential $(a;) vanish for r > a. 

10. Calculate the electrostatic potential <S>(x) outside a thin, hollow sphere of 
radius R on which there is a surface density of charge proportional to 1 + 
3cos(2#). You may assume that there are no boundaries, and that the volume 
charge density is given by 

5 \ r -i 

P(r,0) - [l + 3cos(20)Jj(r-i2), 

where A is a constant. 

What are the components of the quadrupole moment tensor for this charge 
distribution? 
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3.1 Equations of magnetostatics 

For time-independent situations, the part of Maxwell's equations that involve the 
magnetic field are decoupled from those involving the electric field, and they can be 
dealt with separately. The equations of magnetostatics are Poisson's equation for 
the vector potential, 

V 2 A(x) = -n j{x), (3.1.1) 
the Lorenz (or Coulomb) gauge condition 

V-A = 0, (3.1.2) 

and the relation between potential and field, 

B = VxA. (3.1.3) 

We recall that the current density must satisfy the statement of charge conservation, 
V- j = 0. 

Because Eq. (3.1.1) is very similar to the electrostatic equation (2.1.1), the tools 
introduced in Chapter 2 to solve Poisson's equation can be taken over directly to 
magnetostatics. And because we already have gained much experience solving such 
equations, it will suffice to keep our discussion of magnetostatics quite brief. In 
particular, it will be sufficient to solve Eq. (3.1.1) in the absence of boundaries. We 
therefore write the solution as 

A(x) = y~ ( p^TffiV, (3.1.4) 
Air J \x-x'\ 

recalling the results of Sec. 1.8. 



3.2 Circular current loop 

[The material presented in this section is also covered in Sec. 5.5 of Jackson's text.] 
As our first application we calculate the magnetic field produced by a current / 
flowing inside a circular loop of radius R. We place the loop in the x-y plane, at 
z = 0. To evaluate the integral of Eq. (3.1.4) for this situation we will use spherical 
coordinates and express the Green's function as in Eq. (2.6.12), 

1 °° 1 Air r e 

MI^EE ^^(WUM), (3.2.1) 

where r< is the lesser of r and r', while r> is the greater of the two. 
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The density of the negative charge carriers inside the wire is given by 

p{x) = -^5{r-R)5{6-l), (3.2.2) 

where A is the magnitude of the linear charge density (minus the wire's total charge 
per unit length); the total charge is given by J p(x) d 3 x = —2irR\. The charge 
carriers move in the negative <f> direction so as to produce a positive current. Their 
velocity vector is 

v = — v(f> = — sin 0a; + cos<py), (3.2.3) 

where <f> is a unit vector that points in the direction of increasing <f>. Its expression in 
terms of the constant vectors x and y is obtained by noting that the position vector 
of a point in the plane is given by r = r cos <fix + r sin <j>y, so that a vector pointing 
in the (j) direction is dr/dcj) — —r sin ipx + r cos (f>y; dividing by r produces the 
normalized vector <p. From Eqs. (3.2.2) and (3.2.3) we obtain the current density, 

j(x)=pv = ^4>5(r-R)5(e-%), (3.2.4) 

where I = Xv is the current. 

To obtain the vector potential we substitute Eqs. (3.2.1) and (3.2.4) into Eq. (3.1.4). 
For concreteness we take x to be outside the loop, so that r > r' , r< = r' , and 
r> = r. We must pay attention to the fact that <p is not a constant vector: its 
expression in terms of x and y must be substituted inside the integral before inte- 
gration is attempted. It is convenient to consider the complex combination A y —iA Xl 
which is generated by 

j y - ij x = ^(cos^' + ismcf>')6(r' - R)S(6' - § ) = ^ 5(r' - R)6(6' - §); 

we have now expressed this in terms of the primed variables. We notice that the 
</>' dependence of the current density is contained entirely in the exponential factor 
e 1 ^ . After evaluating the r 1 and 9' integrals the vector potential becomes 

f,m " 

The (j>' integral evaluates to Y^ m (§,0) e - im ^' e^' dxj>' = 2nY em (^ , 0)6 mil , and 
the potential simplifies to 

the sum begins at t = 1 because there is no spherical harmonic Y" i- The spherical- 
harmonic functions can then be expressed as Y^i(6,<p) — Y( t \(6, 0)e lc ^, and this 
allows us to extract the real and imaginary parts of the right-hand side. We obtain 

00 o 2 / p\W 

^ - -£ J E ><»{") *WMW 

and 

and these combine to form the vector 

A ^ x ) = iS^E 27TT^(f>°) W Y£ ' lie,0) - (3 ' 2 ' 5) 
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We see that the vector potential is expressed as an expansion in powers of R/r. 

When r 3> R the vector potential is well approximated by the leading term of 
the expansion, given by I = 1. In this case we have Yn(0, 0) = —y/3/(8n) sin# and 
Eq. (3.2.5) becomes 

A(x) ~ g (^i? 2 /) 4 (sin 00) . (3.2.6) 

We will write this in a different form. First we define the vector 

m = (irR 2 I)z, (3.2.7) 

which points in the z direction and whose magnitude is the product of the current I 
and the area ttR 2 enclosed by the loop; this vector is known as the magnetic moment 
of the current distribution. Then we notice that z x x/r = —(y/r)x + (x/r)y = 
— sin 9 sin <fix + sin 9 cos <py = sin 9<f>. So Eq. (3.2.6) can be written as 

A(sB) -£~ ^s— ' (3 ' 2 ' 8) 

and this shall be our final expression. 

The curl of A gives the magnetic field. Apart from a common factor of (/Lto/47r)|m| 
its x component is equal to — X7 z (x/r 3 ) = 3xz/r 5 , its y component is equal to 
— V z (y/r 3 ) = 3yz/r 5 , and finally, its z component is equal to V^rr/r^+V^y/r 3 ) = 
2/r 3 — 3x 5 /r 5 — 3y 2 /r 5 = — 1/r 3 +3z 2 /r 5 . These results combine to form the vector 
3zx/r 5 — z/r 3 , and multiplying by \m\ gives 3(m • x)x/r 5 — m/r 3 . The magnetic 
field is therefore 

%), f (m 'f W (3.2.9) 
47r r 6 

where r = x/r is a unit vector that points in the radial direction. Notice that this 
expression for the magnetic field is very similar to what was obtained in Sec. 2.3 
for the electric field of a dipole: 

_ 1 3(p-f)f-p 
dipolc ~ 47reo r 3 

So the magnetic field far outside a current loop has a dipolar behaviour, with m 
playing the role of dipole moment. 



3.3 Spinning charged sphere 

As a second application we calculate the magnetic field produced by a hollow sphere 
of radius R, uniformly charged, that is rotating with an angular velocity SI = Slz. 
Here the charge density is given by 

p(x) = ^S(r - R), (3.3.1) 

where q is the sphere's total charge. To calculate the velocity vector v of a point x 
on the surface of the sphere, we consider a narrow strip at an angle 9 with respect 
to the z axis (see Fig. 3.1). A point on this strip moves in the (p direction; in a 
time dt the angle increases by d<f> = Vldt. Since this point is at a distance Rsin9 
from the z axis, the length traveled is dl — (R sin 6) d<j>. The speed is therefore 
v = dl/dt = (Rsin6)fl, and the velocity vector is v = vcf> — ClRsin94>, or 

v(x) = tlxx. (3.3.2) 

This is the velocity of the charge distribution, and the current density is given by 

3(x) = ^ 5 ^0&(r-R). (3.3.3) 
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Figure 3.1: A hollow sphere rotating at an angular velocity O around the z axis. 

We recall that <j> = — sin <px + cos <f>y. 

Again it is convenient to consider the complex combination A y — iA x , which is 
generated by 



3y ^Jx 



qfl smB'e^' S(r' -R). 



4ttR 



The angular function happens to be equal to — W8tt/3Yii(9' , </)'). Substituting this 
and Eq. (3.2.1) into the expression for A y — iA x that comes from Eq. (3.1.4), we 
find that the angular integration is nonzero if and only if t = m = 1; the double 
sum over I and m therefore contains a single term. The radial integration reduces 
to 

/•no r 

r 

R 3 /r 2 

and we obtain 



T -^S(r' - R)r' 2 dr' 



r < R 
r>R ' 



iA x — 



Ho q& 47r 



r 

3 /„2 



R /r 



The angular function within square brackets is sin^e 1 ^, and the vector potential is 

A = — — (— sin 9 sin 6x + sin 6 cos d>y) < ^J, 9 
47r3i? V ; \ R 3 /r 2 

After multiplication by r the vector within brackets becomes — yx + xy = z x x. 
If we let 

m = -qQR 2 z = -qR 2 fl, (3.3.4) 
3 3 

we find that the vector potential inside the sphere is given by 

. . /in m x x , „. 
A in (x)= f £^ r - (r<R), 

while the potential outside the sphere is 

Ho m x x 



(3.3.5) 



4tt r 3 



(r > R). 



(3.3.6) 



3.4 Multipole expansion of the magnetic field 
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Notice that these are exact results — no approximations were involved in the cal- 
culation of the vector potential. Notice also that A out takes the same form as in 
Eq. (3.2.8), except for the fact that here, the magnetic moment m is a different 
vector. From this we conclude that the magnetic field J3 ut outside the sphere has 
a dipolar behaviour, and that it takes the form of Eq. (3.2.9). 

To calculate the magnetic field inside the sphere we express the vector potential 
as A ln = (/io/47r)(m/i? 3 )(— yx + xy) and take its curl. Its x and y components are 
both zero, and its z component is proportional to \7 x x + \7 y y = 2. We arrive at 

Mo 2m 

BM = ^- (3 - 3 ' 7) 
This states that inside the sphere, the magnetic field is constant and points in the 
same direction as the magnetic moment (which corresponds to the rotation axis). 



3.4 Multipole expansion of the magnetic field 

[The material presented in this section is also covered in Sec. 5.6 of Jackson's text.] 
For both the current loop (Sec. 3.2) and the rotating sphere (Sec. 3.3) we found 
that the vector potential at a large distance from the current distribution was given 

where the vector m, known as the magnetic moment, characterizes the distribution 
of current. And in both cases we found the magnetic field to be described by 

g, 3(m-*)f-m 
47r r 6 

where r = x/r; this is a dipolar field. Here we will show that these expressions 
arc in fact very general: they are the leading terms in an expansion of the vector 
potential and the magnetic field in powers of 1/r. This multipole expansion is valid 
for points x that lie outside of the current distribution. 

Let the distribution of current be confined to a volume V bounded by a surface 
S. We take x to be outside of V, while x' is necessarily inside. We examine the 
vector potential at a large distance r » r', starting from the exact expression 

47T J \X - X'\ 

We recall that the current density satisfies V • j = 0, and we recognize that there 
is no current flowing across S; this implies that j ■ da = 0, where da is the surface 
element on S. As we shall show below, these equations imply the integral statements 



X 



j(x')d 3 x' = (3.4.3) 

IV 

and 

/ (f • x')j(x') d s x' = - \f x [ [x'x j(x')] d 3 x'. (3.4.4) 
Jv 2 Jv 

To simplify the vector potential we expand l/\x — x'\ in powers of 1/r. We could do 

this very systematically by involving the addition theorem for spherical harmonics 

(Sec. 2.6), but it will be good enough to keep only a small number of terms. For 

this we proceed directly: 





-x') 


■ (x — x') 




- 2x ■ 


x' — r' 2 



— r(l — r ■ x'/r + •••), 
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so that 

Substituting this into the exact expression for A(x) gives 

M0 ^ J j(x') d 3 x' +^J v (r- x')j(x') d 3 x' + ■ ■ 



The first term on the right-hand side vanishes by virtue of Eq. (3.4.3), and the 
second term is altered by Eq. (3.4.4). We have 

= x J v [x>xj(x>)] rfV + .-j, 

and if we define the magnetic moment of the current distribution by 

m= \ I x' x j(x') d 3 x', (3.4.5) 

then we recover Eq. (3.4.1). We see that indeed, 

v ' 4tt r 2 

is a very general expression, and that the magnetic moment of an arbitrary distri- 
bution of current is given by Eq. (3.4.5). It is a straightforward exercise to show 
that for the current loop of Sec. 3.2, Eq. (3.4.5) reproduces Eq. (3.2.7). Similarly, 
Eq. (3.4.5) reduces to Eq. (3.3.4) for the rotating sphere of Sec. 3.3. 

We still have to establish the identities of Eqs. (3.4.3) and (3.4.4). We begin 
with the z component of Eq. (3.4.3), in which we drop the primes on the integration 
variables. Consider 

V • {zj) = z(V ■ j) +j-Vz= j z , 
which follows because V- j = and Vz = z. Integrating this equation over V gives 

/ j z d 3 x = / V ■ (zj) d 3 x = <p zj ■ da, 
Jv Jv ~ Js 

where we have invoked Gauss' theorem. But the right-hand side is zero because no 
current crosses the surface S, and we have J v j z d 3 x = 0. Proceeding similarly with 
the x and y components, we arrive at Eq. (3.4.3). 

To establish Eq. (3.4.4) requires a bit more work, but the idea is the same. We 
consider the z component of this equation, in which we drop the primes on the 
integration variables. (At the same time we replace the vector f, which is constant 
with respect to the primed variables, by an arbitrary constant vector e.) We have 



/ (e • x)j z d 3 x = e x xj z d 3 x + e y / yj z d 3 x + e z / zj z d 
Jv Jv Jv Jv 

for the left-hand side of the identity. Consider now 

V • (xzj) = xz(V ■ j) + j ■ V(xz) = xj z + zj x , 
which integrates to 

/ (xj z + zj x ) d 3 x = f xzj -da = 0, 
Jv Js 



3 x 



3.5 Problems 



53 



giving rise to the integral identity 

/ xj z d 3 x = - / zj x d 3 x. 
Jv Jv 

Similarly, 

/ Viz d 3 x = - / zj y tfx. 
Jv Jv 

On the other hand, 

V • {z 2 j) = z 2 (V • j) + j ■ Vz 2 = 2zj z 

and this yields 



/ zjz 
Jv 



, d 3 x = 0. 

Using these results, we go back to the left-hand side of Eq. (3.4.4) and write 
J (e-x)j z d 3 x = \ x (^J xj z d 3 x - J zj x d 3 x^j 
+ \ e v [ I VJzd 3 x- j zj y d 3 x 



v 



\fx J {zj x - xj z ) d 3 x + ^e y J (yj z - zj v ) d 3 x. 



Inside the first integral we recognize the y component of the vector x x j, while the 
second integral contains its x component. So 



J (e • x)j z d 3 x = ~^e x J {x x j) y d 3 x - e y J (x x j) x d ; 
= -^e x J^(x x j)d 3 x^j , 
and we have arrived at the right-hand side of Eq. (3.4.4). 

3.5 Problems 

1. Show that for the current loop of Sec. 3.2, Eq. (3.4.5) reproduces Eq. (3.2.7). 
And show that Eq. (3.4.5) reduces to Eq. (3.3.4) for the rotating sphere of 
Sec. 3.3. 

2. A uniformly charged sphere of radius R is rotating around the z axis with 
an angular velocity f2. Take the density of charge inside the sphere to be 
3g/(47ri? 3 ) (for r < R), where q is the sphere's total charge. Calculate m, the 
sphere's magnetic moment vector. 

3. Calculate the vector potential A(x) inside the spinning sphere of the preceding 
problem. Express it in terms of the sphere's magnetic moment vector m. 

4. Prove that for a current loop of arbitrary shape that lies in a plane, Eq. (3.4.5) 
gives 

\m\ — (current) (area enclosed by loop), 
and that the direction of m is normal to the plane. 
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Chapter 4 
Electromagnetic waves in 

MATTER 



4.1 Macroscopic form of Maxwell's equations 

[The material presented in this section is also covered in Sec. 6.6 of Jackson's text.] 
Our goal in this chapter is to describe how electromagnetic waves propagate in 
dielectric materials. Our first task will be to reformulate Maxwell's equations in a 
way that conveniently incorporates the electromagnetic response of the medium. 

4-1.1 Microscopic and macroscopic quantities 

Maxwell's equations, as they were written down in Chapter 1, are the fundamental 
equations of classical electrodynamics; in principle they require no modification 
even in the presence of a nontrivial medium. But it is convenient to write them 
in a different form, in which the charges that are free to move within the medium 
are distinguished from the charges that are tied to the medium's molecules, and in 
which the microscopic fluctuations of all the variables involved have been smoothed 
out. The resulting form of Maxwell's equations involves only macroscopic densities 
of free charges and currents, and macroscopic electromagnetic fields. 

To distinguish them from macroscopic quantities, the fundamental, microscopic 
variables will be assigned different symbols in this chapter. The microscopic density 
of charge (including both free and bound charges) is denoted 77, while the micro- 
scopic density of current (including both free and bound currents) is denoted J. 
The microscopic electric field is denoted e, and the microscopic magnetic field is 
denoted b. These quantities are related by the microscopic Maxwell equations. 

V-e = —T), (4.1.1) 
= 0, (4.1.2) 
= 0, (4.1.3) 

= (H> J, (4.1.4) 

which are the same as in Chapter 1. The macroscopic density of free charge will be 
denoted p, while the macroscopic density of free current will be denoted j. These 
are obtained from the microscopic quantities by removing the contribution from 
the bound charges (those that are tied to the molecules and become part of the 
medium's electromagnetic response) and by macroscopic smoothing; the precise 
relationships will be given below. The macroscopic fields will be denoted E, D, 



db 

dt 

V x b - eoMo^r 



V x e + 
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B, and H. The macroscopic electric field E is obtained from e by macroscopic 
smoothing, and D accounts for the electric polarizability of the medium; precise 
relationships will be given below. Similarly, the macroscopic magnetic held B is 
obtained from b by macroscopic smoothing, and H accounts for the magnetization 
of the medium; precise relationships will also be given below. 



4-1.2 Macroscopic smoothing 

Microscopic quantities fluctuate widely over atomic distances, and we would like 
to smooth out these fluctuations to dchne corresponding macroscopic quantities. 
For example, the charge density r\ switches from zero to infinity whenever a point 
charge is encountered, and we would like to replace this by a quantity that varies 
smoothly to reflect only macroscopic changes of density. 

One simple-minded way of smoothing out a microscopic quantity / is to average 
it over a volume V(x) centered at some position x; this volume is imagined to be 
microscopically large (so that all microscopic fluctuations will be averaged over) 
but macroscopically small (so that variations over relevant macroscopic scales are 
not accidentally discarded). The averaging operation is an integration of / over the 
volume V(x): (f(x)) = V^ 1 J f(x') d 3 x', where the domain of integration consists 
of all points x' that lie within V(x). If the volume is a sphere of radius R, then the 
condition is \x — x'\ < R, and (/(#)) can be expressed as 



X - J f(x')9(R-\x-x'\)d 3 x\ 



V 

with a now infinite domain of integration. 

In this averaging operation it is convenient to replace the sharply defined func- 
tion V~ 1 6(R — \x — x'\) by a better-behaved function w(x — x') that smoothly goes 
from V~ x when x' = x to zero when x' is far away from x. Such a function, whose 
detailed specification is not required, will be called a smoothing function] a specific 
choice would be a Gaussian function of width R centered &tx — x' = 0. The smooth- 
ing function is required to be normalized, in the sense that / w(x — x') d 3 x' = 1. 

We define the macroscopic average of a microscopic quantity f(t, x) to be 

F(t,x) = (f(t,x)) = J f(t,x')w(x-x')d 3 x', (4.1.5) 

where w is a suitable smoothing function (assumed to be normalized) and the 
integration is over all space. The property that w does not have sharp boundaries 
implies that the operations of averaging and differentiation commute. For example, 

V„(/(t,a;)) = (V„/(t,aj)) ) (4.1.6) 

and the same is true of time derivatives. This important result follows by straight- 
forward manipulations. Consider, for example, 



dx 



= J f(t,x')^-w(x-x')d 3 x' 



I 



f(t,x')-^w(x-x')d 3 x' 



= [ ^-w(x-x')d 3 x' 
J dx' 



df 

dx 



In the second line we use the fact that w depends on the difference x — x' , so that a 
derivative with respect to x is minus a derivative with respect to x' . In the third line 
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we integrate by parts, noting that the boundary terms at infinity do not contribute 
because w vanishes there; the result is by definition the average of the microscopic 
quantity df/dx. 

Following the usage introduced in Eq. (4.1.5), we define the macroscopic electric 
field to be 

E(t,x) = (e(t,x)), (4.1.7) 
and the macroscopic magnetic field to be 

B(t,x) = (b(t,x)). (4.1.8) 

Because differentiation and averaging are commuting operations, the macroscopic 
form of Maxwell's equations are 

V -E = -( V ), (4.1.9) 

V B = 0, (4.1.10) 
B B 

VxE+— = 0, (4.1.11) 
dE 

VxB-e of i — = fi (J). (4.1.12) 

The homogeneous equations (4.1.10) and (4.1.11) are the same as their microscopic 
version, but the inhomogeneous equations (4.1.9) and (4.1.12) involve a macroscopic 
averaging of the source terms. 



4-1.3 Macroscopic averaging of the charge density 

We distinguish between the charges that can move freely within the medium — the 
free charges — and the charges that are tied to the molecules — the bound charges. 
We will consider the bound charges to be part of the medium and remove them 
from the right-hand side of Maxwell's equations, which will then involve the free 
charges only. The influence of the bound charges will thus be taken to the left-hand 
side of the equations; it will be described by new macroscopic quantities, P (electric 
polarizability of the medium) and M (magnetization of the medium) . From these 
we will define the new fields D = e E + P and H = B//j, — M. 

In Fig. 4.1 we have a picture of the charge distribution in the medium. Let xa (t) 
be the position vector of the free charge qA labeled by "A" , and let VA(t) = dxA/dt 
be its velocity vector. Let x n be the position vector of the centre of mass of the 
molecule labeled by "n" ; for simplicity we assume that the molecules do not move 
within the medium. Finally, let x B (t) be the position vector of the bound charge 
q B labeled by "£?" within a given molecule, relative to the centre of mass of this 
molecule, and let VB(t) — dx B /dt be its velocity vector. 

The density of free charge is 

g bcc (t,x) = ^2q A S(x - x A {t)). (4.1.13) 

A 

The density of bound charge within the molecule labeled by "n" is 

r] n (t,x) = ^2q B S(x - x n - x B (t)). (4.1.14) 

B 

We assume that each molecule is electrically neutral, so that ^2 B q B = 0. The 
density of bound charge is then 



Abound (t, X) = y~] T] n (t, x), 
n 



(4.1.15) 
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Figure 4.1: The left panel shows the free charges (small solid disks) and molecules 
(large open circles) in the medium. The right panel shows the distribution of bound 
charges within a given molecule. 

and the total density of charge is r/ = ?jf ree + Abound- 

We first calculate the macroscopic average of the molecular charge density. We 
have 

{Vn) = J r) n (t,x')w{x - x')d 3 x' 

= / S(x' — x n — xb)w(x — x') d 3 x' 

B J 

= gBw(x - x n -x B )- 

B 

The function w(x — x n — x B ) is peaked at Because the scale R over 

which w varies is very large compared with the intra- molecular displacement xb, 
we can approximate w by Taylor-expanding about x — x n : 

w(x - x n - x B ) = w(x - x n ) - x a B V a w{x - x n ) + -XgX b B VaVbw(x - x') H , 

where to avoid confusing the notation we make the vectorial index on xb a super- 
script instead of a subscript. (Recalling that repeated indices are summed over, we 
have that x B \7 a stands for xb • V.) Substituting this into the previous expression 
yields 

{Vn) = X! ^ BW ~ X! ( lBX l B \7 a w + - ^ qBX c B x b B X7 a X7 b w H , 

B B B 

where w now stands for w{x — x n ). We simplify this equation by recalling that each 
molecule is electrically neutral, so that the first sum vanishes. In the second term 
we recognize the expression for the molecular dipole moment, p n — ^2 B Qb^b, or 

p a n {t) = Y J lBX a B {t), (4.1.16) 

B 

and in the third term we recognize the expression for the molecular quadrupole 
moment, 

Qt(t)=3j2QBX a B (t)x b B (t). (4.1.17) 

B 

Notice that this actually differs from the definition made in Sec. 2.10: here our 
quadrupole moment is not defined as a tracefree tensor. Combining these results 
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gives 

(r] n {t, x)) = -p a n {t)V a w{x - x n ) + ^Q« b {t)V a V b w(x - x n ) + ■ ■ ■ (4.1.18) 

for the molecular charge density. 

The macroscopic average of the density of bound charge is obtained by summing 
(r] n ) over all molecules. This gives 

Abound) = -^KVaW+^<2fV a V 6 W; + --- 

n n 

= -V a ^p> + iv a V 6 ^ + 

n n 

= -Va(j>n<*(* - X n)) + l^a^b(j2 Qn S ( X ~ X n)) + ■ ■ ■ , 
n n 

because w = w(x — x n ) is the macroscopic average of S(x — x n ). The quantities 
P a {t,x) = C^p®(t)5(x — x n )\ = macroscopic polarization (4.1.19) 



and 

Qab{t, x) = (^2, Qn'(t)fi( x ~ X n)^ = macroscopic quadrupole density (4.1.20) 

n 

have a straightforward physical interpretation. Consider, for example, Eq. (4.1.19). 
The quantity within the averaging brackets is ^2 n Pn$( x ~ x n), the sum over all 
molecules of the product of each molecule's dipole moment with a 5-function cen- 
tered at the molecule. This is clearly the microscopic density of molecular dipole 
moments, and P a is its macroscopic average. So the vector P(t, x), called the macro- 
scopic polarization, is nothing but the macroscopic average of the dipole-moment 
density of the medium's molecules. Similarly, Q ab (t,x) is the macroscopic average 
of the molecular quadrupole-moment density. In terms of these we have 

Abound (t, x)) = -V a P a {t, x) + ^V a V b Q ab {t, su) + ■ ■ ■ (4.1.21) 

for the macroscopic average of the density of bound charge. 
The average of the free-charge density is simply 

p(t,x) = (T)tree(t,x)) = (^q A 5{x - X A )) , (4.1.22) 

A 

and the macroscopic average of the total charge density is therefore 

(r,) =p- V ■P+\\7 a V b Q ab + ---. (4.1.23) 


Substituting this into Eq. (4.1.9) yields 

V a (e a E a + P a - ^X7 b Q ab + = p . 

The vector within the large brackets is the macroscopic field D(t, x), and we arrive 
at 

V-D = p, (4.1.24) 
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with the definition 

D a = e E a + P a -\v b Q ab + ---. (4.1.25) 

Notice that this macroscopic Maxwell equation involves only the density of free 
charges on the right-hand side, and that apart from a factor of eo, the effective 
electric field D is the electric field E augmented by the electrical response of 
the medium, which is described by the molecular dipole-moment and quadrupolc- 
moment densities. In practice the term involving Q ab in Eq. (4.1.25) is almost 
always negligible, and we can write D = e n E + P. 

4-1-4 Macroscopic averaging of the current density 
Following the definitions of Eqs. (4. 1.13) -(4. 1.15) we introduce the current densities 

Jf ICC (t,x) = y^q a v A (t)5(x - x A (tj), (4.1.26) 

A 

J n (t,x) = J2 ( lBV B (t)S(x-x n -x B (t)), (4.1.27) 

B 

Jboundit^x) = 5^J„(t,x), (4.1.28) 

n 

and the total microscopic current density is J = Jf roo + Jbound- 
The macroscopic average of the molecular current density is 

( J n) = ^V B W{X-X n -X B ) 
B 

= ^2 q B v B w — q B v B x b B V b w H , 

B B 

where w stands for w(x — x n ) in the second line. In the first term we recognize 
the time derivative of the molecular dipole moment, dp^/dt = (d/dt) J2 B q B x a B = 
J2 B 1b v b- To put the second term in a recognizable form we define the molecular 
magnetic moment by 

m n{t) = ^ ^2qb(x b x v B ) (4.1.29) 

B 

and we note that the vector Vu> x m n has components 
(Vwxm n ) a = ]^q B (x a B v b B -v a B x B )\7 h w 

B 

= i^qB{x B v B +v%x b B )^ b w- y £ j q B v , B x B V b w 

B B 

= l^Q^bW -Y,lBV B x b B V b w, 

recalling the definition = 3^2 B q B x%x b B for the molecular quadrupole moment. 
Combining these results gives 

(Jn) = | (p> ~ \QtV b w^ + (Vw x m n ) a + • • • , (4.1.30) 
and summing over n yields 

Abound) = MY,P>-l Wb Y,^ bw ) + (vx^m nU ;) a + ..., 
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+ V X (^TTln^I - X n ) 



Introducing 

M(t,x) = C^m n (t)S(x — x n ) \ = macroscopic magnetization (4.1.31) 



as the macroscopic average of the magnetic-moment density of the medium's molecules, 
we finally arrive at the expression 

(Abound) = | (Pa g-Qafc) + (V X M)„ + ■ ■ ■ (4.1.32) 

for the macroscopic average of the density of bound current. 
The average of the free-current density is simply 

j(t,x) = (Jhee(t,x)) = (y^q A v A S{x - x A )\ (4.1.33) 



A 

and the macroscopic average of the total current density is therefore 



. d 

3 + ^1 



(d-cqE} + V x M + ---, (4.1.34) 



dt' 

having used Eq. (4.1.25). Substituting Eq. (4.1.34) into Eq. (4.1.12) gives 

dD 

V x (£?/mo - Af) - — = j. 

The vector within brackets is the macroscopic field H(t, x), and we arrive at 

dD 

VxH- — =j (4.1.35) 

with the definition 

H = B/hq- M. (4.1.36) 

Notice that this macroscopic Maxwell equation involves only the density of free 
currents on the right-hand side, and that apart from a factor of /xo, the effective 
magnetic field H is the magnetic field B diminished by the magnetic response of 
the medium, which is described by the molecular magnetic-moment density. 

4-1-5 Summary — macroscopic Maxwell equations 
The macroscopic Maxwell equations are 

V D = p (4.1.37) 

V B = 0, (4.1.38) 

OB 

VxE+— = 0, (4.1.39) 
dD 

VxH- ¥ = j. (4.1.40) 

Here p is the macroscopic free-charge density, and j is the macroscopic free-current 
density. The fields D and H are defined by 

D = e Q E + P (4.1.41) 
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and 

H = —B - M, (4.1.42) 

where 

P = (^Pn{t)5{x-x n )) (4.1.43) 

n 

is the macroscopic polarization (p n = ^2b Qb x b is the molecular dipole moment) 
and 

M= (J2m n (t)5(x - x n )) (4.1.44) 

n 

the macroscopic magnetization (m„ = \ J2b QbXb x «b is the molecular magnetic 
moment). 

4.2 Maxwell's equations in the frequency domain 

[The material presented in this section is also covered in Sec. 7.1 of Jackson's text.] 
We now begin our discussion of electromagnetic waves propagating in a dielectric 
medium. We suppose that the medium is not magnetized, so that we can set M = 
is our equations. We further suppose that the waves are propagating in the absence 
of sources, so that we can also set p = and j = 0. Maxwell's equations therefore 
reduce to 

V D = (4.2.1) 

V B = 0, (4.2.2) 

VxE+— = 0, (4.2.3) 
dD 

VxB-^o— = 0, (4.2.4) 

where D = e E + P. 

We want to form plane-wave solutions to these equations, and examine how 
the waves interact with the medium. We shall assume that the medium is uniform, 
isotropic, and unbounded, and leave unexplored issues of reflection and transmission 
at boundaries. To begin we will look for monochromatic waves oscillating with an 
angular frequency u). We therefore set 

(4.2.5) 
(4.2.6) 
(4.2.7) 
(4.2.8) 

with the understanding that only the real parts of these complex fields have physical 
meaning. 

The restriction to monochromatic waves is not too severe. Because the field 
equations are linear (assuming, as we shall verify below, that the relation between 
D and E is itself linear), solutions with different frequencies can be added to form 
wave-packet solutions. For example, we can express the electric field as 

E(t,x) = -— E{uj, x)e- wt duj, (4.2.9) 

V27T J-oo 

in terms of the frequency-domain field E(u; 7 x); this is a Fourier representation of 
the electric field. The inverse transformation is 

1 f°° 

E(w,x) = -= E{t,x)e wt dt, (4.2.10) 

V27T J-oo 



E(t,x) 


= E(x)e~ 


-iujt 


D(t,x) 


= D{x)e 


-iujt 


P(t,x) 


= P{x)e~ 


-iujt 


B{t,x) 


= B(x)e 


-itut 
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and frequency-domain fields D(u), x), P(u>,x), B(to,x) can be defined a similar 
way. 

The frequency-domain fields satisfy the equations 



V • D 


= 


(4.2.11) 


V B 


= o, 


(4.2.12) 


V x E - iuoB 


= o, 


(4.2.13) 


V x B + iujuoD 


= o, 


(4.2.14) 



and we have the relation 

D = e Q E + P. (4.2.15) 

These are the sourceless Maxwell equations in the frequency domain. They can only 
be solved once an explicit relationship (called a constitutive relation) is introduced 
between the polarization P and the electric field E. The constitutive relation 
describes the medium's response to an applied electric field, and its derivation is 
based on complicated molecular processes. While quantum mechanics is required 
for a proper treatment, in the next section we shall consider a simple classical model 
that nevertheless captures the essential features. 



4.3 Dielectric constant 

[The material presented in this section is also covered in Sec. 7.5 of Jackson's text.] 
We suppose that the medium is sufficiently dilute that the electric field felt 
by any given molecule is just the applied field E itself, the field exerted by other 
molecules being negligible. We assume that in the absence of an applied field, the 
molecular dipole moment is zero, at least on average. When an electric field is 
applied, however, some of the molecular charges (the electrons) move away from 
their equilibrium positions, and a dipole moment develops. 

The displacement of charge q B from its equilibrium position is denoted £ B (t). 
This charge undergoes a motion that is governed by the applied electric field (which 
tends to drive £b away from zero) and the intra-molecular forces (which tend to drive 
the charge back to its equilibrium position). We model these forces as harmonic 
forces, all sharing the same natural frequency loq. What we have, therefore, is a 
system of simple harmonic oscillators of natural frequency loq that are driven at 
a frequency lo by an applied force. Since the motion of the charges will produce 
electromagnetic radiation, energy will gradually be removed from the oscillators (to 
be carried off by the radiation), and the oscillations will be damped; we incorporate 
this effect through a phenomenological damping parameter 7 > 0. The equations 
of motion for the charge q B are therefore 

m B (g B +"f£ B +ujte) =q B E(uj,x n )e- i " t , (4.3.1) 

where overdots indicate differentiation with respect to t. Notice that the applied 
electric field, which is macroscopic and varies slowly over intra-molecular distances, 
is evaluated at the molecule's centre of mass. To solve the equations of motion 
in the steady-state regime we write £ B (t) = £ B e~ luJt and we substitute this into 
Eq. (4.3.1). This yields 

— or — vyu 
for the frequency-domain displacement. 
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n 

E(w,x)(y2S(x - x n ) 



The molecule's dipole moment is p n — ^ B qs{xB + £b), where Xb is the equi- 
librium position of charge qs- Because the dipole moment vanishes at equilibrium, 
this reduces to p n = ^2 B qs^B- Substituting Eq. (4.3.2) then gives 

UJq — Ld Z — 1-fUJ 

Since all the contributing charges qs are electrons, we can write YIb^bI^-b — 
Ze 2 /m, where Z is the number of electrons per molecule, e the electronic charge, 
and m the electron's mass. We therefore arrive at 

Ze 2 1 

Pn(w) = 2 2 — X n ) (4.3.3) 

for the molecular dipole moment in the frequency domain. 
The polarization vector is then 

^2p n S(x - x n ) 

n 

Ze 2 1 

m uj 2 — to 2 — ijuj 

Ze 2 1 

m to 2 id 2 iju 

The quantity within the averaging brackets is the microscopic molecular density. Its 
macroscopic average is therefore the medium's density N (the number of molecules 
per macroscopic unit volume), which we take to be uniform. We have arrived at 

Ze 2 N 1 

P(u, x) = , 3 ; — E(lu, x) (4.3.4) 

for the frequency-domain polarization. We recall that Z is the number of electrons 
per molecule, that N is the molecular number density, and that e and m are the 
electron's charge and mass, respectively. We notice that the relationship between 
P and E is linear, as was anticipated in Sec. 4.2. 

The relationship between D = e E + P and E also is linear. We express it as 

- / Ze 2 N 1 \~ 
D = e 1+ - 2 - 2 —)E, 

or as 

D(lo, x) = e(oj)E(w, x), (4.3.5) 
having introduced the (complex, and frequency-dependent) dielectric constant 

^ = 1 + 2 4 ■ (4-3.6) 
and the medium's plasma frequency 



;Ze 2 N\ 1/2 (Ao „ 
= —ZT- ■ ( 4 - 3 - 7 ) 



e m 



Equation (4.3.6) is the product of a simplistic classical model. But a proper 
quantum-mechanical treatment would produce a similar result, except for the fact 
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Figure 4.2: Real and imaginary parts of the dielectric constant, as functions of 
lu/luo. To fit both curves on the same graph we removed 1 from Re(e/eo). The real 
part of the dielectric constant goes through 1 at resonance. The imaginary part of 
the dielectric constant achieves its maximum at resonance. 

that the electrons actually oscillate with a discrete spectrum of natural frequen- 
cies uij and damping coefficients jj. A more realistic expression for the dielectric 
constant is then 



where fj is the fraction of electrons that share the same natural frequency ojj . In the 
rest of this chapter we will continue to deal with the simple expression of Eq. (4.3.6). 
It is good to keep in mind, however, that the model is a bit crude. 

In Fig. 4.2 we plot the real and imaginary parts of the dielectric constant. You 
will notice that for most frequencies, when u> is not close to Wq, Re(e/eo) increases 
with increasing w, and Im(e/eo) is very small; this typical behaviour is associated 
with normal dispersion. When u is very close to uiq, however, the electrons are 
driven at almost their natural frequency, and resonance occurs. Near resonance, 
then, Re(e/eo) decreases with increasing uj, and Im(e/eo) is no longer small; this 
behaviour is associated with anomalous dispersion. In realistic models the number 
of resonances is larger than one: there is one resonance for each natural frequency 

At low frequencies the dielectric constant behaves as e/e ~ 1 + (w p /wo) 2 , and 
is therefore a constant slightly larger than unity. At high frequencies it behaves as 
e/eo — 1 — (uj p /lS) 2 , and therefore approaches unity from below. 

4.4 Propagation of plane, monochromatic waves 

[The material presented in this section is also covered in Sec. 7.1 of Jackson's text.] 
After making the substitution of Eq. (4.3.5), the frequency-domain Maxwell 
equations (4.2.11)-(4.2.14) become 




V E = 







(4.4.1) 
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V B = 0, 
V x E - icuB = 0, 
V x B + iuje(u)fi E = 0. 



(4.4.2) 
(4.4.3) 
(4.4.4) 



We recall that the time-domain fields are reconstructed, for example, as E — Ee~ lult 
if the wave is monochromatic, or as in Eq. (4.2.9) if the wave possesses a spectrum 
of frequencies. 

It is easy to show that as a consequence of Eqs. (4.4.1)-(4.4.4), both E and B sat- 
isfy a frequency-domain wave equation. For example, taking the curl of Eq. (4.4.3) 
gives 

V x (V x E) - iujV x B = 0; 

substituting Eq. (4.4.4) and making use of the vectorial identity V x (V x E) = 
V(V • E) - V 2 E puts this in the form 



where 



(V 2 + k 2 )E(oj,x) = 0, 



k 2 = ^ 2 e H 
c 2 e 



We have introduced a complex wave number k(u>), and we recall that c : 
is the speed of light in vacuum. Similar manipulations return 

(V 2 + k 2 )B(uj,x) = 0, 



(4.4.5) 

(4.4.6) 

(6oMo)- 1/2 

(4.4.7) 



the statement that the magnetic field also satisfies a wave equation. 

To see how monochromatic electromagnetic waves propagate in a dielectric 
medium, consider any component of E or B, and call it ip. Then suppose that 
the wave is a plane wave that propagates in the z direction, so that ip = ip(z). The 
wave equation simplifies to 



dz 2 



+ k 2 )i){w,z) = 0. 



(4.4.8) 



A right-moving wave is described by the particular solution ip = e . To help 
identify the physical properties of this wave we substitute Eq. (4.3.6) into Eq. (4.4.6) 
and get 



k 2 = 



1 + 



uJq lo 2 ijuj 



(4.4.9) 



To be concrete we suppose that the second term within the large brackets is small. 
Taking the square root gives 



k 



1 + 



or 



where 



k R = -{l + 



2(w 2 — lo 2 — i-fu;) 
k = kn + ikj, 



(4.4.10) 



2[(u 



LO 



2\2 



2,. ,21 



is the real part of the complex wave number, while 
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is its imaginary part. 

The right-moving wave is therefore described by 

ip(t, z) = 4>(u, z)e- lut = e -kiz e -i("t-k R z)_ (4.4.11) 

From this expression we recognize that the wave is indeed traveling in the positive 
z direction, with a speed 

v p = r-^— = (4.4.12) 

known as the phase velocity; n = ckji/uj is the medium's index of refraction. From 
our previous expression for Ur we obtain 

1 . u 2 p {u 2 -u 2 ) 
n = 1 + 



2[(wg - lo 2 ) 2 + 7 2 w 2 ] ' 

and we see that n ~ 1 + |(ajp/aj ) 2 > 1 for low frequencies, so that v p < c. On the 
other hand, n ~ 1 — i(w p /a>) 2 < 1 for high frequencies, so that v p > c. The phase 
velocity of the wave can therefore exceed the speed of light in vacuum! This strange 
fact does not violate any cherished relativistic notion, because a monochromatic 
wave does not carry information. To form a signal one must modulate the wave 
and superpose solutions with different frequencies. As we shall see, in such situations 
the wave always travels with a speed that does not exceed c. 

Another feature of Eq. (4.4.10) is that the wave is an exponentially decreasing 
function of z provided that ki is positive, which is guaranteed to be true if 7 > 0. 
A positive damping constant means that some of the wave's energy is dissipated by 
the oscillating electrons (which emit their own electromagnetic radiation), and the 
wave must therefore decrease in amplitude as it travels through the medium. It is 
interesting to note that in the case of lasers and masers, the electronic oscillations 
actually reinforce the wave and lead to an exponential amplification instead of an 
attenuation; in such media the damping constant is negative. 

Having constructed solutions to Eqs. (4.4.5) and (4.4.7), we should examine 
what constraints the full set of equations (4.4.1)-(4.4.4) place on these solutions. 
We write 

E(u, z) = eE (cj)e lk ^ z (4.4.13) 

and 

B(u, z) = bB (cj)e ik ^ z , (4.4.14) 

where e and b are real unit vectors, and where E and B are complex amplitudes. 
Substituting these expressions into Eqs. (4.4.1) and (4.4.2) informs us that the unit 
vectors are both orthogonal to z, 

e-z = b-z = Q; (4.4.15) 

the fields are therefore transverse to the direction in which the wave propagates. 
Equation (4.4.3) then produces B n b = (k/oj)E (z x e), which implies 

B q (lu) = ^Eo(w) (4.4.16) 

and 

b = z x e. (4.4.17) 

These equations state that the magnetic field is orthogonal to the electric field, and 
that its amplitude is linked to that of the electric field. Finally, Eq. (4.4.4) gives us 
B (z x 6) = —(k/uj)E e, which merely confirms Eqs. (4.4.16) and (4.4.17). 
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4.5 Propagation of wave packets 

[The material presented in this section is also covered in Sees. 7.8 and 7.9 of Jack- 
son's text.] 

4-5.1 Description of wave packets 

In this section we shall assume, for simplicity, that e(uj) is real: we work away from 
resonances and ignore dissipative effects. We let k(ui) be the positive solution to 
Eq. (4.4.6). 

A wave packet is obtained by superposing solutions to Eq. (4.4.8) with different 
frequencies lj. Including both left- and right-moving waves, we write 
1 r°° 

' a+{uj)e lk ^ z + a-{oj)e- lkl -^ z 



(4.5.1) 



where cl±{lo) are complex amplitudes that determine the shape of the wave packet. 

Alternatively, and more conveniently for our purposes, we can instead superpose 
solutions with different wave numbers k. Defining uj(k) to be the positive solution 
to Eq. (4.4.6), we write 



i/>(t,z) 



3 lkz dk, 



(4.5.2) 



V^TT J-o 

and it is not too difficult to show that the decomposition of Eq. (4.5.2) is equivalent 
to that of Eq. (4.5.1). 

The complex amplitudes A±(k) are determined by the initial conditions we wish 
to impose on ip(t, z) and its time derivative. It follows from Eq. (4.5.2) that 

1 r°° 

il>(0,z) = -j= J [A+{k) + A-(k)]e lkz dk, 
and Fourier's inversion theorem gives 

1 f°° 

A+(k) + A_(k) = -= ^{0,z)e- lkz dz. 



On the other hand, differentiating Eq. (4.5.2) with respect to t produces 
4>(0,z) = —= / iuj(k)[A+(k) - A_{k)]e ikz dk, 

V27T J-oo 



and inversion gives 



iui{k)[A+(k) - A-(k)] = 



1 



Solving for A±(k), we finally obtain 



A±(k) = - 



iuj{k) 



ikz 



dz. 



e~ lkz dz. 



(4.5.3) 



The procedure to follow to construct a wave packet is therefore to choose an ini- 
tial configuration by specifying tp(0,z) and tp(0, z), then calculate A±(k) using 
Eq. (4.5.3), and finally, obtain ip(t, z) by evaluating the integral of Eq. (4.5.2). 

In the following discussion we will choose the initial configuration to be time- 
symmetric, in the sense that 

i>(0,z) = 0. (4.5.4) 

Equation (4.5.3) then simplifies to 

1 1 f°° 

A(k) = A ± (k) = --;= I ip{0,z)e~ lkz dz, (4.5.5) 

J — OO 



2^2^ 

and we have an equal superposition of left- and right-moving waves 
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4-5.2 Propagation without dispersion 

Let us first see what the formalism of the preceding subsection gives us when there 
is no dispersion, that is, when n = c\k\/w is a constant independent of k. We then 
have 

u(k) =v p \k\, (4.5.6) 

where v p = c/n is the wave's phase velocity. 

For time-symmetric initial data, Eq. (4.5.2) reduces to 



iP(t,z) 



1 



A(k) 



i\k\v p t _|_ g— i\k\v p t 



? kz dk. 



The factor within the square brackets is 2cos(|fc|w p t) = 2cos(kv p t), and we can 
therefore replace the previous expression with 



V>(M) = 



/27T 



A{k) 



e ik(z+v p t) _|_ e ik(z-Vpt) 



dk. 



But from Eq. (4.5.5) we have 

1 1 f'°° 
A(k) = --= mz')e~ lkz dz', 

and substituting this gives 

1 f°° 1 r°° r / / 

ll>(t,z) = - dz'tp(Q,z')— e ik(z+v p t-z ) + e ik(z-v p t-z ) 



dk 



dz'ip(0, z') S(z + v p t - z') + S(z - v p t - z') 



or 



i>{^ z ) = ^(0, z + v p t) + iv(o, z - Vpt). 



(4.5.7) 



The first term on the right-hand side is (half) the initial wave packet, ip(0,z), 
translated in the z direction by — v p t; it represents that part of the wave packet 
that travels undisturbed in the negative z direction with a speed v p . The second 
term, on the other hand, is the remaining half of the wave packet, which travels 
toward the positive z direction with the same speed v p . 

The main features of wave propagation without dispersion are that the wave 
packet travels with the phase velocity v p — w/k and that its shape stays unchanged 
during propagation. 

4-5.3 Propagation with dispersion 

Let us now return to the general case of propagation with dispersion. Here lu is a 
nonlinear function of k, and for time-symmetric initial data we have 



iP(t,z) 



1 f°° 

-7= / A ( k ) 
v2tt J -oo 



giu>(k)t _|_ ^ — iuj(k)t 



ikz 



dk. 



To be concrete, suppose that A(k) is sharply peaked around k — k , with /c some 
arbitrary wave number. [A proper treatment would have A(k) be a symmetric 
function peaked about both k = ko and k = — ko. Unlike what we shall find below, 
this would produce a real wave function tp(t,z). It is not difficult to make this 
change, but we shall not do so here.] Suppose further that ui(k) varies slowly near 
k = ko, so that it can be approximated by 



u>(k) ~ ujq + ui (k - k ) = (ui - u' k o) + u' k, 



(4.5.8) 
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1.04 



1.03 



1.02 



1.01- 




Figure 4.3: Group velocity v g in a dielectric medium, in units of c, as a function 
of (jj/u>q. 



where luq = uj(ko) and uj' = [du j 'dk)(ko) . Then the wave packet can be expressed 
as 



i 



27T J- 



A(k) 



e ik(z+uj' t) e i(uj -uj' k )t , e ik(z-u' t) e ~i(uJo-uj' k )t 

A {j.yk{zW t) dk 



dk 



or 



V2tt J- 

+ e -i{u -u k )t / A(k)e tk{z -^ l) dk, 

v27r j-oo 

^{t, z) = e 4 ^- w ofeo)t^( 0; z + + e -i(«-o-«i*o)t^( 0) 2 _ ^ (4.5.9) 



where 



duj 



(4.5.10) 



is known as the group velocity. From Eq. (4.5.9) we learn that in the presence 
of dispersion, the wave packet travels with the group velocity v g instead of the 
phase velocity v p — fco/^o- And the phase factors indicate that the shape of the 
wave suffers a distortion during propagation. These are the main features of wave 
propagation in a dispersive medium. 

If the dispersion relation is expressed as k — n(u>)u)/c, in terms of the index of 
refraction n(u>), then dk/duj — c _1 [n + u>(dn/du))] and 



1 



c n + uj(dn/duj) 
For the dielectric medium described in Sec. 4.3, 

u 2 p {ul -w 2 ) 



(4.5.11) 



n A = 1 



(cl>q — LO 2 ) 2 + 7 2 w 2 ' 



and the group velocity is plotted in Fig. 4.3. You will notice that v g is smaller 
than c for most frequencies, but that it seems to exceed c near resonance. This 
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phenomenon is entirely artificial: when w is close to ujq it is no longer true that ui 
is a slowly varying function of k, and the approximations involved in defining the 
group velocity are no longer valid. Furthermore, the imaginary part of u>, together 
with the dissipativc effects it represents, cannot be neglected near resonance. The 
bottom line is that the very notion of group velocity is no longer meaningful near a 
resonance. 



4-5.4 Gaussian wave packet 

As a concrete example of wave propagation with dispersion we consider an initial 
configuration given by 

-0(0, z) = e - z2 ' {2L ^ cos(k z); (4.5.12) 

this is an oscillating wave (with wave number k ) modulated by a Gaussian envelope 
of width L. Once more we take the initial data to be time-symmetric, and we 
impose 0(0, z) = 0. In order to be able to carry out the calculations we shall take 
the dispersion relation to be 

u{k) = v[ 1 + \o?k 2 ), (4.5.13) 



where v is a minimum frequency and a an arbitrary length scale. This dispersion 
relation is not particularly realistic, but its quadratic form will allow us to evaluate 
the integrals that will appear in the course of our analysis. These will have the form 
of Gaussian integrals, and we quote the standard result 

£V- 2+ - +c d*=^e*p(£ + c). (4.5.14) 

Here the constants a, b, and c can be complex, but Re(a) must be positive for the 
integral to converge. 

We first calculate the Fourier amplitudes associated with Eq. (4.5.12). Equation 
(4.5.5) gives 

A(k) = --= / tp(0,z)e- tkz dz 



2 V2^ 
1 1 



e -*7(2£") 



e -i(fc-fe )z _|_ e -i(k+k )z 



dz. 



For the first integral we have a = 1/(2L 2 ), b = —i(k — fc ), c = 0, and b 2 /(Aa) = 
— \{k — ko) 2 L 2 ; for the second integral we simply make the substitution k — > — fc 
in this result. Equation (4.5.14) then gives 

A (k) = | [ e -^ fe - fc °) 2L2 + e -^+M 2 L 2 ] . (4 . 5 . 15) 

We see that A(k) is a symmetric function of k, and that it is peaked at both k = fc 
and k = —k . The width of the distribution around both peaks is equal to l/L. 
The wave packet at times t ^ is given by Eq. (4.5.2), 

i f'°° i f'°° 

V> = -= / A{k)e l[kz -^ k)t] dk+^ A ( k y[kz+u(k)t\ dk ^ 



which we may also express as 

= -= / A{k)e^ kz -^ k ^ dk + — = \ A(-k)e-^ kz -^- k ^ dk, 
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Figure 4.4: Wave propagation with a high-frequency dispersion relation u>(k) = 
(ujp + c 2 k 2 y/ 2 , calculated numerically with the help of a Fast Fourier Transform 
routine. The upper panel shows the initial wave profile, described by Eq. (4.5.12). 
The lower panel shows two versions of the propagated wave. The first (lower curve) 
is obtained without dispersion by setting lu p = 0; this wave propagates with speed 
c without altering its shape. The second (higher curve — the wave function is 
translated upward for ease of visualization) is obtained with dispersion; this wave 
spreads out as it travels with speed v g < c. Notice the drop in amplitude as the 
initial wave packet splits into an equal superposition of left- and right-moving waves. 
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having changed the variable of integration to — k in the second integral. Because 
A(— k) = A(k) and uj(—k) = u(k), we have 



1 



2tt 
2 



A(k)e i[kz " (fc)tl dk + complex conjugate, 

) 

/CO 
A(k)e l[kz -^ k)t] dk. 
-co 



Substituting Eq. (4.5.15) then gives 



2V2tt 



Re 



e 2 



Taking into account Eq. (4.5.13), we see that the argument of the exponential 
function is a second-order polynomial in k, which we write as —ak 2 + bk + c, with 
a = \{L 2 + ivta 2 ), b = k L 2 + iz, and c = — (iut + \k^L 2 ). With this we have 

_ (z-ikpL 2 ) 2 
4a ~ 2L 2 {l + ivtf 2 y 

where / = a/L, and Eq. (4.5.14) yields 



i> = ^Rc<j {l + ivtf 2 )- 1/2 exp 



(z - ik L 2 ) 



2\2 



-ivt-ikiL- 1 



+ (fco 



2L 2 {l+ivtf 2 ) 

It is then a matter of straightforward algebra to bring this to its final form, 



fco)}. 



%l>{t,z) = -ReUl + wi/ 



where / = a/L and 



-1/2 



exp 



2L 2 {l + ivtf 2 



J[k a z-uj(ko)t] 



(k -» -fe ) k 



va 2 k 



dbj 
dk 



(feo) 



(4.5.16) 



(4.5.17) 



is the group velocity at fc = fco for the dispersion relation of Eq. (4.5.13). 

Equation (4.5.16) informs us that the wave packet retains its Gaussian shape as 
it propagates through the medium, but that its width increases with time: 



L(t) 



(4.5.18) 



The wave therefore spreads out, and this is a typical feature of propagation with 
dispersion. In addition, we see that while the oscillations travel at the phase velocity 
v p = u)(ko)/ko, the centre of the Gaussian peak travels at the group velocity v g ; 
this result was expected after the discussion of the preceding subsection. 

A more realistic example of wave dispersion (calculated numerically) is displayed 
in Fig. 4.4. Notice that the features discussed in the preceding paragraphs are 
present in this example as well: the wave packet spreads out as it propagates with 
a group velocity v g that is smaller than c. 



4.6 Problems 

1. (This problem is adapted from Jackson's problem 7.19.) A wave packet in 
one dimension has an initial shape given by ^(0, z) = f(z)e tk ° z , where f(z) is 
the modulation envelope; the wave packet is assumed to be time-symmetric. 
For each of the forms f(z) given below, calculate the wave-number spec- 
trum |A(fc)| 2 of the wave packet, sketch |^»(0, z)\ 2 and |A(fc)| 2 , calculate the 
root-mean-square deviations Az and Afc from the means (defined in terms of 
1^(0, z)\ 2 and |^4(fc)| 2 , respectively), and test the inequality AzAk > |. 
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a) f(z) = e -«M/2; 

b) f(z) = e-« 2 - 2 /4. 

c) f(z) = 1 — a\z\ for a\z\ < 1 and f(z) = for a\z\ > 1; 

d) f(z) = 1 for a\z\ < 1 and f(z) = for a\z\ > 1. 

2. We have seen that the relationship between the frequency-domain fields D 
and E is 

D{uo,x) = e(oj)E(u>,x), 



where 



= 1 + 



Prove that as a consequence of this, the relationship between the time-domain 
fields D and E is 



D(t,x) 



eo 



E{t,x) + -^ J e-~i T l 2 sin{vT)E{t-T,x)dT 



where z/ = y — \ r ) 2 - This reveals that £) at time t depends on the electric 
field at all previous times: the relationship is nonlocal in time. The exponen- 
tial factor, however, indicates that delays longer than I/7 do not matter very 
much: only the electric field's recent history is important. 



Chapter 5 
Electromagnetic radiation 
from slowly moving sources 



5.1 Equations of electrodynamics 

In this chapter and the next we will explore some of the mechanisms by which 
electromagnetic radiation is produced. While Chapter 6 will be concerned with 
relativistic situations, here we shall restrict our attention to charge and current 
distributions whose internal velocities are small compared with the speed of light. 

We recall the main equations of electrodynamics, as we have formulated them 
in Chapter 1. We have wave equations for the scalar and vector potentials, 

\3$(t,x) = -—p(t,x) (5.1.1) 

and 

nA(t,x) = -fi j{t,x), (5.1.2) 

where 

□ = -^+V* (5.1.3) 

is the wave operator and c = (eo/ i o)~ 1/ ' 2 is the speed of light. For Eqs. (5.1.1) 
and (5.1.2) to be equivalent to Maxwell's equations, the potentials must satisfy the 
Lorenz gauge condition, 

1 <93> 

V.A + --=0. (5.1.4) 

Maxwell's equations imply that the charge and current densities satisfy the conti- 
nuity equation, 

^+V-j=0, (5.1.5) 

which states that charge is locally conserved. 

In the absence of boundaries, the wave equations have solutions 

*"-"=4^/ p ""'r.:tr"' ^ 



and 



these are retarded solutions, which incorporate the physically appropriate causal 
relationship between sources and fields. 
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From the potentials we obtain the fields via the relations 

r)A 

E=-^--V<S> (5.1.8) 

dt v ; 

and 

B = VxA. (5.1.9) 

And finally, we recall from Sec. 1.4 that the fields carry energy. The electromagnetic 
field energy density is given by 

e(t,x) = le Q E 2 + ^-B 2 (5.1.10) 

and the energy flux vector (Poynting vector) is 



S(t, x) = —E x B. (5.1.11) 



5.2 Plane waves 



The simplest situation featuring the production of electromagnetic waves involves 
monochromatic plane waves. These are generated by a planar distribution of 
charges, all oscillating in the plane with the same frequency ui. 

We take the charge distribution to be uniform in the x-y plane, so that the 
charge density is given by 

p(t,x) = (t5(z), (5.2.1) 

where a is a constant surface density. Notice that the charge density is actually 
time independent. The resulting scalar potential will also not depend on time, and 
the electromagnetic waves will be described entirely by the vector potential. 
We take each charge to be oscillating in the x direction with a velocity 

v (t, x) = v xcosujt. (5.2.2) 

Notice that the velocity vector is actually independent of position in the plane, and 
that x is actually an arbitrary direction within the plane; we could pick another 
direction by rotating the x-y coordinates around the z axis. The current density is 

j(t,x) = <Tv 5(z)xcosLdt, (5.2.3) 

and we will calculate the associated vector potential without approximations. 

In principle we might proceed by substituting Eq. (5.2.3) into Eq. (5.1.7), but 
it is easier to start instead with the equivalent form 

A(t,x) = ^ [jit' ^-'-^-^ dM, 
4-7T J \x — x'\ 

in which we substitute Eq. (5.2.3) expressed in terms of primed variables. We have 
A(t,x) = ^-<jv x / dt'T(t — t',x) coswi', 

47T J 



where 



T(t-t',x)= I S(z') S{t - t '- lx - xVc) d^ 

\x — x'\ 



is an integral that we must now evaluate. 
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Integration over dz' is immediate: 

6{t -t' - y/(x - x') 2 + (y - y'Y + z 2 /c) 
y/(x - x') 2 + (y - y') 2 + ^ 2 



dx'dy', 



or 



7 ^V 2 +y' 2 + z 2 



after translating the origin of the (x',y') coordinates to the point (x,y); this is 
allowed, because the domain of integration is the whole x'-y' plane. To integrate 
this we make the change of variables x' — £cos9, y' — £sin#. This gives, after 
writing dx'dy' — ididO and integrating over d6, 

2 s {t -e-VWTT V c )ldt 



Jo 



jo s/W^z- 

The easiest way to integrate over the (^-function is to make its argument be the 
integration variable. So let s = V ' l 2 + z 2 /c — At, where At = t — t' . Because I is 
nonnegative, we have that s is bounded below by \z\/c— At. We also have 

idl Idl 

ds 



cV£ 2 + z 2 c 2 (s + Ai) : 
and the integral becomes 



/>oo 

r = 2ttc / 5{s)ds, 

J\z\/c-At 



l\z\/c-At 

or 



T(t - t', x) = 2nc9(t - t' - \z\/c), 

where 9(s) is the step function, which is unity when s > and zero otherwise. 

Substituting this result into our previous expression for the vector potential, we 
obtain 

A(t,x) = ^(<7v )(2irc)x [ dt'6(t- \z\/c - t') cosut' 

47r J-oo 



Ho r*-\ z \/ c 

— 2-KavQCX I cosut' dt'. 

47T 7-oo 



The dt' integral evaluates to w _1 sinw(t~ |^|/c), apart from an ambiguous constant 
associated with the integral's infinite lower bound. Since this constant will never 
appear in expressions for the electric and magnetic fields, we simply set it to zero. 
Our final expression for the vector potential is therefore 

A(t,z) = — smw(t- \z\/c)x, (5.2.4) 

where 

E o = -^Hovvoc (5.2.5) 

has the dimensions of an electric field. For z > we have that sinw(i — \z\/c) = 
sinw(t — z/c) = — sinfc(z — ct), where k = ui/c is the wave number, and we see 
that the wave travels in the positive z direction with speed c. For z < we have 
instead sinaj(t — |z|/c) = sinu(i + z/c) = sin k{z + ct), and the wave now travels in 
the negative z direction. Thus, Eq. (5.2.4) describes a plane wave that propagates 
away from z = 0, the place where it is generated; the wave has a frequency lo, a 
wave number k — ui/c, and it travels with the speed of light. 
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The time-dependent part of the electric field is obtained from Eqs. (5.1.8) and 
(5.2.4): 

E(t, z) = -E cos cj(t - \z\/c)x. (5.2.6) 

We see that the electric field behaves as a plane wave that travels away from z = 
with the speed of light, and that it points in the x direction, which is perpendicular 
to the direction of propagation. The magnetic field, on the other hand, is obtained 
by taking the curl of Eq. (5.2.4). This involves differentiating \z\ with respect to z, 
which is +1 for z > and — 1 if z < 0. Denoting this function by e(z) = 0(z)—6(—z), 
we have 

E 

B(t, z) = —e(z) cosuj(t - \z\/c)y. (5.2.7) 

c 

The magnetic field also is a plane wave that travels away from z = at the speed of 
light, and it points in the y direction, which is perpendicular to both the direction 
of propagation and the direction of the electric field. Notice that E and B arc 
both proportional to cosuj(t — \z\/c): the fields are in phase. Notice also that \B\ 
is smaller than \E\ by a factor 1/c. 

To finish off our discussion of plane electromagnetic waves, we calculate the 
field's energy density and flux vector. From Eqs. (5.1.10), (5.2.6), and (5.6.7) we 
obtain 

e(t, z) = e El cos 2 uu{t - \z\/c), (5.2.8) 
and from Eq. (5.1.11) we get 

S(t,z)=ce(z)e(t,z)z. (5.2.9) 

This shows that the electromagnetic field energy travels with the wave, in the pos- 
itive or negative z direction, away from the source, and with the speed of light. 
Let us summarize the main points: 

• The plane wave travels in the (positive or negative) z direction with speed c. 

• The electric and magnetic fields are perpendicular to z and to each other. 

• The fields are in phase. 

• Their magnitudes are related by |B|/|£7| = 1/c. 

• The electromagnetic field energy travels with the wave. 

These properties are generic: they are shared by all electromagnetic waves. The 
only caveat is that in general, the direction of propagation is not uniform: the 
constant vector £ goes into a vector n that may change from point to point. 



5.3 Spherical waves — the oscillating dipole 

[The material presented in this section is also covered in Sees. 9.1 and 9.2 of Jack- 
son's text.] 

5.3.1 Plane versus spherical waves 
Any function of the form 

W,z) = f(tTz/c), (5.3.1) 

where /(£) is an arbitrary function of £ = t =F z/c, is a solution to the wave equation 
Dtp = 0. If £ = t — z/c, the plane wave travels in the positive z direction with speed 
c; if instead £ = t + z/c, the wave travels in the negative z direction. To see that 
Eq. (5.3.1) does indeed form a valid solution to the wave equation, we note that 
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dip/dt = /'(£) while dip/dz — T /'(£,) I c , where a prime indicates differentiation with 
respect to £; substitution into the wave operator yields Dip = —cr 2 f + c~ 2 /" = 0, 
as required. 

Plane waves are not a very good approximation to realistic waves, which orig- 
inate in a localized region of space and propagate outward with a decreasing am- 
plitude. A better approximation that incorporates these features is the spherical 
wave 

i>(t,r) = I/(t- r /c), (5.3.2) 

where f(u) is an arbitrary function of u = t — r/c. This is also a solution to Dip = 0. 
Here the wavefronts are spherical instead of planar, and the wave propagates out- 
ward in the radial direction; its amplitude decays as 1/r, giving rise to an intensity 
that goes as 1/r 2 . That Eq. (5.3.2) forms a valid solution to the wave equation 
follows by direct substitution: In the absence of an angular dependence the wave 
operator reduces to 

Id 2 1 d , d 

a = 1 r 2 — . 

c 2 dt 2 r 2 dr dr 

Differentiating with respect to t yields d 2 ip/dt 2 = r -1 /, while we have r 2 dip/dr — 
—rf/c — / and 

d 2 diP 1. 1=1; 1 = 

7T r IT =—f+ ~^ r f + ~f = -^2 r -f' 
Or or c c z c tr 

where an overdot indicates differentiation with respect to u. All this gives Dip = 0, 
as required. Note that a similar calculation would reveal that ip = r~ x f{t + r/c) 
also is a solution to the wave equation; this wave travels inward instead of outward. 



5.3.2 Potentials of an oscillating dipole 

The simplest electromagnetic spherical wave is produced by an oscillating dipole 
located at the origin of the coordinate system. Here and in the following subsections 
we will calculate the potentials and fields produced by this dipole, and examine the 
nature of the corresponding electromagnetic radiation. 

In Sec. 2.3 we saw that the charge density of a point dipole could be expressed as 
p{x) = — p- V<5(x), where p was a constant dipole moment vector. Here we introduce 
a time dependence by making p oscillate with angular frequency u: p — > pcosuut. 
The charge density of an oscillating dipole is then 

p(t,x) = — p ■ "V5(x) cosut. (5.3.3) 

Notice that in this equation, p is still a constant vector, which we align with the z 
axis: p = pz. 

A time- varying charge density must be associated with a current density j, 
because 

V • j = — — = —ujp ■ V8(x) smut = V • [— ujpS(x) sinwf] . 

The simplest solution to this equation is 

j(t,x) = —ujpS(x) sinwi, (5.3.4) 



and this is the current density of an oscillating dipole. 

We first compute the vector potential produced by the current density of Eq. (5.3.4). 
Once more it is convenient to start with the expression 

A(t, x) = ^ [j(t>, x>) S{t X ' l/C} dfdV. 

4-k J \x — x'\ 
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After substituting Eq. (5.3.5) we obtain 



A(t,x) 



Mo 
Mo 

47T 



Lup / dt' sin tot 



'/ S(x') 



S(t - f - \x - x'\/c) d3x , 



IX — X' 



Lup / dt' smut 



,S(t-t' -r/c) 



or 



Mo sinw(£- r/c) 

A(t, x) — — —up . 

47T r 



(5.3.5) 



This has the structure of a spherical wave, with its expected dependence on u — 
t — r/c and its amplitude decaying as 1/r. 

We may now compute the scalar potential. This calculation is a bit more labo- 
rious, and we start with 



$(t,x) 



1 



4ire 



IX — X' 



After substituting Eq. (5.3.3) we get 

<P(t,x) = -— *— p- I dt'T(t-t',x) coswt', 
47re J 



where 



T(t-t',x) = J [V'S(x')] 



ix(~J\\ *(* ~t' -\x-x'\/c) 



IX — X' 



is a vectorial integral that we must now evaluate. 
Integration by parts allows us to write 



= - J S(x')V' 



5(t-t' - \x-x'\/c) 



d 3 x'. 



Because the dependence on x' is contained entirely in \x — x'\, and because V'|a; 
x'\ = —(x — x')/\x — x'\, evaluating the gradient gives 



S(t — t' — \x- x'\/c) 1 S^t - t' - \x - x'\/c) 



IX — X 



/|2 



IX — X' 



x — X 
\x — x' 



d 3 x' 



S(t -H - r/c) 1 5'{t -t'- r/c) 



or 



T(t-t',x) = -f 



S(t — t'- r/c) 1 S'(t -t'- r/c) 
r z c r 



where r — x/r. 

Substituting this result into our previous expression for the scalar potential, we 
obtain 



$(t,x) = - ^—p-r [ dt' 
47re J 



S(t -t' - r/c) 1 S'(t -f- r/c) 
r 2 c r 



cos ut' 



and integrating against the (5-functions gives 

1 



$(i,a:) 



47re 



p ■ r 



cosu>(t — r/c) to sincj(t — r/c) 
r 2 c r 



One must be careful with the second term: a minus sign comes from the integration 
by part needed to convert S'(t — t' — r/c) into an actual S- function, but this minus 
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sign is canceled out because the prime on the <5-function indicates differentiation 
with respect to t — t' — r/c instead of t' . Wc express our final result as 



where 



and 



$ = *noar + $wavc, (5.3.6) 

1 p • V 

$ nea , T (t,x) = 5- cos(wu) (5.3.7) 

$wavct,a: ) = — A p-f— ^ '-, 5.3.8 

47re c r 

and where u = t — r/c is retarded time. 

The scalar potential contains two terms. Of these, only $ wavo has the mathe- 
matical structure of a spherical wave. The other term, $ ne ar, has the structure of a 
static dipole field, 

1 pf 



^dipole (x) 



to which a time dependence cos(wu) is attached. Unlike $ W avc, this term does not 
represent electromagnetic radiation. Instead, $ noar describes the field of a point 
dipole, and this field follows the dipole's time dependence with the usual delay: the 
cos ut dependence of the dipole is replaced by cosw(i — r/c). 

As the name indicates, <& n car dominates when r is small, while $ waV c dominates 
when r is large. The scalar potential becomes purely radiative when r is sufficiently 
large that <& n0 ar can be neglected; the region of space for which this is true is called 
the wave zone. On the other hand, the region of space for which r is sufficiently 
small that $ W ave can be neglected in front of $ n ear is called the near zone. 

In the near zone we must have 1/r 2 3> u>/(cr), or 

c 1 A 

r <C — = 7 = — (near zone), (5.3.9) 

where k is the wave number and A the wavelength. In the near zone, u = t — r/c ~ t 
and the delay contained in the retarded time is hardly noticeable — the near-zone 
field responds virtually instantaneously to changes in the source. In the wave zone 
we have instead 1/r 2 <C uj/(cr), or 

c 1 A 

r >• — = - = — (wave zone). (5.3.10) 

LO k Z7T 

Here the time delay is important, and the scalar potential behaves as a true spherical 
wave. In the region of space where r is comparable to c/lo, between the near zone 
and the wave zone (a region that is sometimes called the induction, or intermediate, 
zone), $near and $wavc are comparable, and the scalar potential behaves neither as 
an instantaneously-changing dipole field nor as a spherical wave. 



5.3.3 Wave-zone fields of an oscillating dipole 

Let us calculate the electric and magnetic fields in the wave zone. We use the 
expressions 

4 / n Mo sin(am) . „ _ . 

Alt, x) = -^ujp — ^ — '- 5.3.11) 
47r r 

and 

$(t,x) = -- p-f—^-, 5.3.12 

47re c r 

where u — t — r/c is retarded time. Notice that Eq. (5.3.11) is exact, while 
Eq. (5.3.12) is valid in the wave zone only. 
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To do this calculation we will need to compute the gradient of 

, , , sin(om) 
tf(t,r)= — 

We have V?/> = (dtp / 'dr)V V = (dtp/dr)r and 

dtp sin (ant) o> cos(uju) 
dr r 2 c r ' 

because du/dr = 1/c. But since we are working in the wave zone, the first term 
can be neglected and we obtain the useful approximation 

v shuW) = _a; cos(anQ ._ 
r c r 

We might also need the gradient of p ■ r. To see what this is, consider its x compo- 
nent: 

/ x\ pxx p x 
V K \p- - 



r j r r 

so that 

V(p-f) = p - {p -^ . (5.3.14) 

r 

The fact that this is of order 1 jr means that the gradient of p ■ f will be negligible 
in our forthcoming calculations. 

The electric field is related to the potentials by Eq. (5.1.8), 

£ = -— - V$. 
dt 

Substituting Eqs. (5.3.11) and (5.3.12), we obtain 



_ /in 9 COs(ont) 1 U> 

E = ^lu 2 P 1 >- + 

47r r 47reo c 



sin(am)„. „. . *s_sin(ant) 
— - — - V(p • f) + (p ■ r)V — - — - 



Using Eq. (5.3.14) shows that the first term within the square brackets is of order 
1/r 2 and can be neglected in the wave zone, while using Eq. (5.3.13) gives 

_, 9 cos(om) Id 2 . „, cos(uju) „ 

E = ^u?p —^-L - - - p • f)—^-f. 

47r r 47re & r 

Collecting the terms and recalling that e /j, = 1/c 2 brings this to its final form 

£wavc(i, x) = £V [p-(p- r)r] (5.3.15) 

This has the structure of a spherical wave, and we notice that E wavc is transverse 
to the direction in which the wave propagates: 

#wavc ■ f = 0. (5.3.16) 

If, for concreteness, we let p point in the z direction, then p — pz, f = sin 9 cos <px + 
s'mdsiiKpy + cos 9z, p ■ r = pcos9, and p — (p ■ f)f = — psm9(cos9 coscpx + 
cos 8 sin (py— sin 6z) = —p sin 80, where 9 is a unit vector that points in the direction 
of increasing 9. This yields 

E wave (t,x) = -^Lu 2 psm9 C ° S ^ JU) -e (5.3.17) 
47r r 

for the electric field. 



5.3 Spherical waves 



— the oscillating dipole 



83 



The magnetic field is obtained by taking the curl of the vector potential. We 
have 

^/ x Mo _sin(o;u) 
B(t,x) = ^ W pxV^^, 
47r r 

and Eq. (5.3.13) gives us 

^ , x «n w 2 . cos(o>u) „ 
S wave (t,x = ^— rxp 5.3.18 

The magnetic field also is a spherical wave, and like the electric field it is transverse 
to the direction of propagation: 

-B wavc f = 0. (5.3.19) 

Furthermore, we see from Eqs. (5.3.15) and (5.3.18) that the fields are mutually 
orthogonal, because [p — (p • r)f] -(rxp) = p • (f x p) — (p • r)r ■ (r x p) = 0, since 
the vector r x p is perpendicular to both p and r. So the fields satisfy 

-E^wave ' -Bwavc = (5.3.20) 

in the wave zone. For p ~ pz we have r x p = — psin#(— sin 05 + cos^y) = 
— p sin #</>, and Eq. (5.3.18) becomes 

_ . . «nW 2 . .COs(wu) ~ ,„„„„n 

Bwavc t,aO = -T^— Psintf — - 0. 5.3.21 

47r c r 

5.5.^ Poynting vector of an oscillating dipole 

We now calculate S = ^ 1 E x B, the energy flux vector associated with the 
electromagnetic radiation produced by the oscillating dipole. It is easy to check 
that in the wave zone, the electric and magnetic fields are related by 

E wavc — cB wavc x r. (5.3.22) 

This follows because the magnetic field is proportional to r x p, and (f x p) x f = 
(f ■ f)p — (r ■ p)r = p — (p-r)f, which gives the direction of the electric field. Using 
once more the same vectorial identity, the Poynting vector is 

S = — (B x r) x B = — \(B ■ B)r - (B ■ r)B 

Mo Mo L 

The second term vanishes because B is orthogonal to r, and we arrive at 

S wa , ve (t, x ) = \B wave {t, x) I r. (5.3.23) 
Mo 

This shows that the energy is transported in the radial direction, the same direction 
in which the wave propagates. 

For p = pz we can substitute Eq. (5.3.21) into Eq. (5.3.23) and obtain 

c / ^ \ 2 ^ 4 P 2 sin 2 9 cos 2 (am) „ 
\47ry c 2 r 2 

n p 2 u; 4 sin 2 9 o/ \ - 
= 167r2c r 2 cob (am) r. 

The energy flux oscillates in time, and after averaging over a complete wave cycle 
we get 

ic> \ MoP 2 ^ 4 sin 2 9 „ 
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From the angular dependence of this expression we see that most of the energy is 
emitted near the equatorial plane (9 = ^), and that none of the energy propagates 
along the z axis (9 = or 9 = n). The average of any function / of retarded time u 
over a complete wave cycle is defined by (/(«)) = T^ 1 f(u) du, where T = 2-k/uj 
is the oscillation period. Using this rule we get (cos 2 (ivu)) = |, and this gives 
Eq. (5.3.24). 

The total power radiated (energy per unit time) is obtained by integrating S 
(energy per unit time per unit area) over any closed surface S. The averaged power 
is then 

(P) = j>(S)-da. (5.3.25) 

It is simplest to take S to be a sphere of constant r. Then da = fr 2 sin 9 d9d(j), and 
Eq. (5.3.25) becomes 



32tt 2 c 

The integral evaluates to 4/3, and we obtain 



(P) = Jg^V. (5.3.26) 

5.3.5 Summary — oscillating dipole 

We have seen that in the wave zone (r 3> c/u), the fields of an oscillating dipole 
are given by 

„/ n \ Mo 9 „ cos u)(t — r/c) ~ 

E(t, r, 6,(p) = - ^-lu 2 P sin 9 v — 9 

47r r 

and 

n ,x i*nw 2 . „ cosoj(t — r/c) T 

B(t,r,6,<f>) = - — —psm9 v '—^ 0. 

4ir c r 

These have the structure of a spherical wave, which travels outward at the speed of 
light. They are produced by a charge distribution whose dipole moment vector is 
given by p(t) — pcos(ut)z. The field amplitudes are related by |-B|/|.E| = 1/c, and 
both fields are in phase. The electric field is directed along 

6 = cos 9 cos (j) x + cos 9 sin <fi y — sin 9 z 
while the magnetic field is directed along 

4> = — sin <p x + cos <p V ■ 
The fields are mutually orthogonal, and they are both transverse to 

r = sin 9 cos <fi x + sin 9 sin </> y + cos 9 z, 

the direction in which the wave is traveling. This is also the direction in which 
the field energy is transported. The angular profile of the energy-flux (Poynting) 
vector is described by sin 2 9, so that most of the energy is emitted in the directions 
orthogonal to the dipole. 



5.4 Electric dipole radiation 

[The material presented in this section is also covered in Sec. 9.2 of Jackson's text.] 
We now consider a much more general situation, in which electromagnetic ra- 
diation is produced by an arbitrary distribution of charges and currents, with an 
arbitrary time dependence (not necessarily oscillating with a single frequency ui). 
Our only restrictions are that 
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• the source is confined to a bounded region V of space; 

• the charges are moving slowly. 

These conditions will allow us to formulate useful approximations for the behaviour 
of the electric and magnetic fields. Although the present context is very general, we 
will see that the results to be derived in this section are very close to those obtained 
in the preceding section. 



5.4-1 Slow-motion approximation; near and wave zones 

To make the slow-motion approximation precise, and to define near and wave zones 
in the general case, we introduce the following scaling quantities: 

r c = characteristic length scale of the charge and current distribution, 
t c = characteristic time scale over which the distribution changes, 

T 

v c = — = characteristic velocity of the source, 

tc 

u c = — = characteristic frequency of the source, 

tc 

A c = = ct c = characteristic wavelength of the radiation. 

The characteristic length scale is defined such that if p(t,x') and j(t,x') are the 
charge and current densities, respectively, then \x'\ is at most of order r c throughout 
the source. This implies that the distribution of charge and current is localized 
within a region whose volume is of the order of rj?. The characteristic time scale is 
defined such that dp/dt is of order p/t c throughout the source. If, for example, p 
oscillates with a frequency f2, then t c ~ l/Cl and lu c ~ f2. 

The slow-motion approximation means that v c — r c /t c is much smaller than the 
speed of light: 

v c < c. (5.4.1) 
This condition gives us r c = v c t c <C ct c = A c , or 

r c < A c . (5.4.2) 

The source is therefore confined to a region that is much smaller than a typical 
wavelength of the radiation. 

We can now define the near and wave zones: 

A c 

near zone: \x — x'\ <§C — - = — , (5.4.3) 

A c 

wave zone: \x — x'\ >• — - = — . (5.4.4) 

Because r c A c , the wave zone can also be defined by the condition r 3> X c /(2tt). 
We will compute the potentials and fields in these zones, using Eqs. (5.4.1)-(5.4.4) 
to simplify our expressions. 



5-4-2 Scalar potential 
We begin by calculating the scalar potential, 



47re Jv \ x ~ x \ 
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where V is the region of space occupied by the source. In the near zone we can 
treat \x — x'\/c as a small quantity and Taylor-expand the charge density about the 
current time t. Because p(t — e) = p{t) — p{t)e + ^p(t)e 2 + ■ ■ ■ we have 



p(t-\x-x'\/c) = p(t)--p(t)\x 



where an overdot indicates differentiation with respect to t. Relative to the first 
term, the second term is of order \x — x'\/(ct c ) = \x — x'\/X c , and by virtue of 
Eq. (5.4.3), this is small in the near zone. The third term would be smaller still, 
and we neglect it. We then have 

*(t,x) = -M/ f^ld 3 x'-- [ P (t,x')d 3 x> + . 
47re [Jv \x -x'\ c J 



Aire 
1 

47re 



Jv \x-x'\ cdtJv 



The second term vanishes, because it involves the time derivative of the total charge 
J p(t, x') d 3 x', which is conserved. We have obtained 



^noar(i, x) 



4ire 



p(t,x') 
\x — x'\ 



1 + 



\x — X 



A? 



dV. 



(5.4.5) 



We see that the near-zone potential takes its usual static expression, except for the 
small correction of order \x — x'\ 2 /X 2 and the fact that charge density depends on 
time. The time delay between the source and the potential has disappeared, and 
what we have is a potential that adjusts instantaneously to the changes within the 
distribution. The electric field it produces is then a "time-changing electrostatic 
field" . This near-zone field docs not behave as radiation. 

To witness radiative effects we must go to the wave zone. Here \x — x'\ is large 
and we can no longer Taylor-expand the density as we did previously. Instead we 
must introduce another approximation technique. We use the fact that in the wave 
zone, r is much larger than r', so that 

\x — x'\ = \J (x — x') ■ (x — x') 
= \J r 2 — 2x ■ x' + r' 2 

1 -2f ■ x'/r + 0{r 2 /r 2 ) 
= r 1 - r ■ x'/r + 0{r 2 c /r 2 ) 



1/2 



This gives 



where 



r — r ■ x' + 0{r 2 /r). 



p(t — \x — x'\/c) = p(t — r/c + f ■ x' /c + e/c) 
= p(u + f-x'/c + e/c), 

u = t — r/c 



(5.4.6) 

is retarded time and e is of order r 2 c jr. Let us now Taylor-expand the charge density 
about the retarded time u instead of the current time t. We have 



p(t -\x- x'\/c) = p(u) + -p(u)(r ■ x' + e) H , 

where an overdot now indicates differentiation with respect to u. We see that 
relative to r ■ x', the e term is of order r c /r <C 1 and can be neglected. Now relative 
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to p{u), the time-derivative term is of order r c /{ct c ) = v c /c, and this is small by 
virtue of Eq. (5.4.1). We therefore have 



p(t - \x - x'\/c) = p(u) + 

in the wave zone. 

Inside the integral for $ we also have 



-p(u) + 0(v 2 /c 2 ) 



(5.4.7) 



1 + f ■ x'/r + 0(r 2 /r 2 ) 



1 + 0{r c /r) 



\x — x'\ 

and it is sufficient to keep the leading term only. The wave-zone potential is then 



*(t,x) 



47T60 r Jy 

1 1 " 
47T£o T 



p(u) + 1 -^p(u) + ■ 



p(u, x') d 3 x 



c 

3^' 



d 3 x' 



r d 

C Oil Jy 



/ p(u, x')x' d 3 x' 
Jv 



In the first integral we recognize Q, the total charge of the distribution; this is 
actually independent of u by virtue of charge conservation. In the second integral 
we recognize p(u), the dipole moment vector of the charge distribution; this does 
depend on retarded time u. Our final expression for the potential is therefore 



} (t,x) 



1 Q If- p(u) 



where 



is the total charge and 



47re r 47re cr 

Q = I p(u,x')d 3 x' 
Jv 

p(u) = / p(u, x')x' d 3 x' 
Jv 



+ ■ 



(5.4.8) 
(5.4.9) 
(5.4.10) 



the dipole moment. The first term on the right-hand side of Eq. (5.4.8) is the static, 
monopole potential associated with the total charge Q. This term does not depend 
on time and is not associated with the propagation of radiation; we shall simply 
omit it in later calculations. The second term, on the other hand, is radiative: it 
depends on retarded time u = t — r/c and decays as 1/r. We see that the radiative 
part of the scalar potential is produced by a time-changing dipole moment of the 
charge distribution; it is nonzero whenever dp/dt is nonzero. You should check that 



^wave(^j 3?) 



1 f • p(u) 
4ire cr 



reduces to Eq. (5.3.8) when the dipole moment is given by p(t) = po cos(wt), where 
p is a constant vector. 

5.4-3 Vector potential 
The exact expression for the vector potential is 



A(t,x) 



Mo 

4-7T 



Jv 



(t- \x-x'\/c,x') 
\x — x'\ 



d 3 x>. 



In the near zone we approximate the current density as 



j(t-\x-x'\/c)=j(t) 



1 + 
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and we obtain 

r 

. fz—jj 1 ^( ' 



^ncar^, x) . [ , t\ 
47T Jy \X - x'\ 



d A x. (5.4.11) 



We see that here also, the vector potential takes its static form, except for the 
correction of order \x — x'\/X c and the fact that the current density depends on time. 
The potential responds virtually instantaneously to changes in the distribution, and 
there are no radiative effects in the near zone. 
In the wave zone we have instead 

j(t- \x-x'\/c) = j(u + f-x'/c) = j(u)[l + 0(v c /c)] 

and 

A wave (t,x) = ^- f j(u,x')[l + 0(v c /c)]d 3 x'. 

In static situations, the volume integral of j vanishes — recall Eq. (3.4.3). But 
here the current density depends on time, and as we shall prove presently, we have 
instead 

[ j(u,x')d 3 x' =p(u). (5.4.12) 
Jv 



The vector potential is therefore 

A wavc (t,x) = ^^+--- (5.4.13) 

47T r 

in the wave zone. This has the structure of a spherical wave, and we see that 
the radiative part of the vector potential is produced by a time-changing dipole 
moment. You may verify that Eq. (5.4.13) reduces to reduces to Eq. (5.3.5) when 
the dipole moment is given by p(t) = po cos(wi), where p is a constant vector. 

To establish the z component of Eq. (5.4.12) we start with the identity V- (zj) = 
z (^ ' j) + j ' (^ z )j m which we substitute the statement of charge conservation, 
V • j = —dp/dt. After integration over the source we have 

[ j z d 3 x — -f- f pz d 3 x = [ V • (zj) d 3 x = J) zj ■ da, 
Jv dt J v J V Jg 

where S is the two-dimensional surface bounding V. Because no current is crossing 
this surface, the right-hand side vanishes and we have 

d f ^3 dp z 
j z a x = — / pzd x = ——. 
v dt J v dt 

This is the same thing as the z component of Eq. (5.4.12), except for the fact that 
p and j are expressed in terms of the variables (t,x) instead of (u, x'). The other 
components of the equation are established with similar manipulations. 



5.4-4 Wave-zone fields 

$wave(i, 3j) — ~T~~ L ^ k + 0(v c /c)} (5.4.14) 

47re cr L J 

A wavc (t, x) = [l + 0(v c /c)} , (5.4.15) 
An r I J 

where u = t — r/c and f = x/r, are generated by time variations of the dipole 
moment vector p of the charge and current distribution. They therefore give rise 



The potentials 



and 
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to electric- dipole radiation, the leading-order contribution, in our slow-motion ap- 
proximation, to the radiation emitted by an arbitrary source. We now compute the 
electric and magnetic fields in this approximation. 

To get the electric field we keep only those terms that decay as 1/r, and neglect 
terms that decay faster. For example, when computing the gradient of the scalar 
potential we can neglect Vr _1 = — r/r 2 ; we only need V(r -p), which we calculate 
as 

~ .\ ( x y z 

Vx{r-P) = Vx[ ~Px + ~Py + ~Pz 

\r r r 

= -Px(--Vxr) + -Py(--V x r) + -p z [--V x r) +0[- 
r \ c J r \ c ) r \ c J \r 

= — (r-p)-+0 - 

c r \ r 



so that 



V(r-p) =--(r-p)r + o(-\. (5.4.16) 



c ' \r 



To obtain this result we made use of the fact that p depends on x through u = 
t — r/c, so that acting with V amounts to taking a derivative with respect to u and 
multiplying by -c^'Vr = — cr x r. The electric field is then 

E = — V$ 

dt 

1 (r ■ p)r 

47T r 47re c 2 r 

or 

^wave(^) = _g» A") -['•#")]' = ^ Nj^! , (5.4.17) 

Att r 4tt r 

Notice that the wave-zone electric field behaves as a spherical wave, and that it is 
transverse to r, the direction in which the wave propagates. 

To get the magnetic field we need to compute V x p(u). To see what this is we 
examine its x component: 

1 .. y 1 .. z 

= -~Pz- + ~Py- 

c r c r 

= --( fx pL> 



so that 



The magnetic field is then 



V x p(u) = —-f x p(u). (5.4.18) 

c 



B wavc M)^^. (5.4.19) 
47r cr 

Notice that the wave-zone magnetic field behaves as a spherical wave, and that 
it is orthogonal to both r and the electric field. Notice finally that the fields are 
in phase — they both depend on p(u) — and that their magnitudes are related 
by = 1/c. The properties listed in Sec. 5.3.5 are therefore shared by all 

electric-dipole radiation fields. 
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5.4-5 Energy radiated 

The Poynting vector is S = Hq 1 E x B, and since E = cB x f in the wave zone 
- as follows from Eqs. (5.4.17) and (5.4.19) — we have S — (c/^, )(B x f) x B = 
(c/ ' no)[\B\ 2 f — (B ■ f)B\. The last term vanishes by virtue of Eq. (5.4.19), and we 
obtain 

C i i 2 

Swavc = -Bwavc f, (5.4.20) 

MO ' ' 

the same result as in Eq. (5.3.23). The fact that the Poynting vector is directed 
along r shows that the electromagnetic field energy travels along with the wave. 

The energy crossing a sphere of radius r per unit time is given by P = § S ■ da, 
where da = rr 2 dVL and dVL = sin0d9d<j). Substituting Eqs. (5.4.19) and (5.4.20) 
yields 

To evaluate the integral we use the trick of momentarily aligning the z axis with 
the instantaneous direction of p(u) — we must do this for each particular value of 
u. Then r x p(u) = \p(u)\sm9 and 



''" l p(u)\ 2 <f> sin 2 9 dn = -p^\p{u)\ 2 f * d(f> r sa?6dB. 
J (47r) c J q Jq 



r2n 

The angular integration gives (27r)(4/3) and we arrive at 

\p(u)\ 2 . (5.4.21) 



Mo 

67TC 1 



This is the total power radiated by a slowly-moving distribution of charge and 
current. The power's angular distribution is described by 

where Q(u) is the angle between the vectors r and p(u). 

5.4-6 Summary — electric dipole radiation 

To leading order in a slow-motion approximation, the wave-zone fields of an arbi- 
trary distribution of charges and currents are given by 

B mvc (t, X ) = - ^ Xpit - r/c) \l + 0( Vc /c) 
47T cr L 

and 

B wavc — cB wavc x 7*. 

The fields are generated by time variations of the source's dipole moment vector, 

= / p(t,x)xd 3 x. 
Jv 

The fields behave as spherical waves, they are each transverse to f, and they are 
mutually orthogonal. The field energy flux vector is 

S -—IB I 2 r 

MO 

and the total power radiated is given by 
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Figure 5.1: Electric field lines of an oscillating dipole, p(t) = pocos(u>t)z, at 
selected moments of time, for u> = 1. Witness the transition between near-zone 
behaviour and wave-zone behaviour which occurs near r = A = 2ttc/lu = 1. The 
curves were obtained by following the method detailed in Sec. 2.11. 

The power's angular profile is proportional to sin 2 0(u), where 0(u) is the angle 
between the vectors f and p(u) . 

The electric field lines of an oscillating dipole, p(t) — po cos(ojt) z, are rep- 
resented in Fig. 5.1. The field lines of a rotating dipole, p(f) = po[cos(ujt)x + 
sin(cc>i) y], are represented in Fig. 5.2. 

5.5 Centre-fed linear antenna 

[The material presented in this section is also covered in Sec. 9.4 of Jackson's text.] 
As an example of a radiating system we consider a thin wire of total length 2£ 
which is fed an oscillating current through a small gap at its midpoint. The wire 
runs along the z axis, from z = — I to z = £, and the gap is located at z = (see 
Fig. 5.3.) 

For such antennas, the current typically oscillates both in time and in space, 
and it is usually represented by 

j(t lX ) = Ism[k(£- \z\)]6(x)6(y)zcosujt, (5.5.1) 

where 

k = uo/c. (5.5.2) 

The current is an even function of z (it is the same in both arms of the antenna) 
and its goes to zero at both ends (at z = ±£). The current at the gap (at z = 0) is 
Isink£, and I is the current's peak value. 

We want to calculate the total power radiated by this antenna, using the electric- 
dipole approximation. To be consistent we must be sure that v c <C c for this 
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Figure 5.2: Electric field lines of a rotating dipole, p(t) — p [cos(uit) x + sm(ojt) y], 
at selected moments of time, for ui = 1. Witness the transition between near-zone 
behaviour and wave-zone behaviour which occurs near r = A = 2ttc/lu = 1. 



coaxial feed 



antenna 



z=-l 



Figure 5.3: Centre-fed linear antenna. The oscillating current is provided by a 
coaxial feed. 
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distribution of current, or equivalently, that £ = r c <C A c = 2ir/k — refer back to 
Eq. (5.4.2). In other words, we must demand that 

M< 1, 

which means that k\z\ is small throughout the antenna. We can therefore approxi- 
mate sin[fc(^— |,z|)] by k{£ — \z\), and Eq. (5.5.1) becomes 

j(t,x) = I (l - \z\/£)5(x)5{y)zcosLut, (5.5.3) 

where 

I = Ik£ (5.5.4) 

is the value of the current at the gap. In this approximation the current no longer 
oscillates in space: it simply goes from its peak value 1$ at the gap to zero at the 
two ends of the wire. In Sec. 9.4 of his book, Jackson gives an exact treatment of 
the antenna of Eq. (5.5.1). 

To compute the power radiated by our simplified antenna we first need to cal- 
culate p(t), the second time derivative of the dipole moment vector. For this it is 
efficient to turn to Eq. (5.4.12), 

P(t) = J j(t,x)d 3 x, 

in which we substitute Eq. (5.5.3). We have 

pit) = I zcosuit / (l — \z\/£) dz, 
J -i 

and evaluating the integral gives 

pit) = Iq£z cos uit. 

Taking an additional derivative and making the substitution t —* u = t — r/c yields 

p\u) = — Iouj£z sin(wu), 

or 

p\u) = — (I c){k£)z sin(wu) (5.5.5) 

after using Eq. (5.5.2). 

Since the vector p(u) is always aligned with the z axis, we have that 0(m) = 8, 
where 8 is the angle between the vectors r and z. According to Eq. (5.4.22), 
therefore, the angular distribution of the power is given by 



dfl (4tt) 2 c 



p(u)\ sin 8 



^—{I c)\Mf sin 2 {urn) sin 2 



(4tt) 2 c 

After averaging over a complete wave cycle, this reduces to 

To obtain the total power radiated we must integrate over the angles. Using 
§ sin 2 8dQ = 87r/3, we arrive at our final result 

(P) = ^(lok£) 2 . (5.5.7) 

For a fixed frequency u), the power increases like the square of the feed current Iq. 
For a fixed current, the power increases like the square of the frequency, so long 
as the condition kl <C 1 is satisfied. From Eq. (5.5.6) we learn that most of the 
energy is radiated in the directions perpendicular to the antenna; none of the energy 
propagates along the axis. 
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5.6 Classical atom 

The classical picture of a hydrogen atom contains a proton at the centre and an 
electron on a circular orbit around it. Supposing for simplicity that the proton is 
at rest, the system's dipole moment is p(t) = —ea(t) where a(t) is the electron's 
position vector. Since the electron is accelerated by the Lorentz force F exerted by 
the proton, we have p(t) = —ea = —(e/m)F. Because this is nonzero, the atom 
emits electromagnetic radiation. The atom therefore loses energy to the radiation, 
and as a result, the orbital radius a gradually decreases; it eventually goes to zero. 
As is well known, the classical model leads to an unstable atom. In this section wc 
calculate the lifetime of the classical atom. 
The Coulomb force exerted by the proton is 

1 e 2 - 
F = -a, 



where a = a(t)/a. This gives 

p(t) = - i —-^a(t) (5.6.1) 
4ne ma 2 

for the second derivative of the dipole moment vector. According to the electric- 
dipole approximation, the total energy radiated by the atom per unit time is given 
by Eq. (5.4.21), 

D A*0 I .., n |2 



Substituting Eq. (5.6.1) then gives 



6«: (47re ) 2 TOW K ' 

Energy conservation dictates that the energy carried away by the radiation must 
come at the expense of the atom's orbital energy; we therefore have 

= -P. (5.6.3) 

The orbital energy of a classical atom is 

^ = -4^ fa ^ 

Combining the last three equations, it is easy to show that the orbital radius must 
decrease, at a rate given by 

da _ 4 1 e 2 

dt " _ 3(4^e ) 2 m 2 c3 a 2- i 5 ' 6 ' 5 ) 

To obtain the lifetime we integrate Eq. (5.6.5), starting with the initial condition 
that at t = 0, the orbital radius is equal to Bohr's radius 

47re fi 2 
«o = s-- 



(In this step we make a leap and momentarily leave the classical realm, since the 
Bohr radius is very much a consequence of quantum mechanics. But we may 
still think of a as giving an appropriate initial condition, as it sets the atomic 
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length scale.) Integration of Eq. (5.6.5) is elementary, and we obtain r = JJ" dt — 
J° o (dt/da) da as 

(47re ) 5 c 3 /i 6 



T = 

4me 10 

To put this in friendlier terms we introduce 

e 2 1 

a = — — — ~ — - = fine structure constant (5.6.6) 

(47re )fic 137 v ; 

and 

Ac = — 2.43 x 10~ 12 m = electron's Compton wavelength. (5.6.7) 
mc 



The lifetime can then be expressed as 

t=-^-~1.6x10- u s. (5.6.8) 
8Tra b c 

This result shows that under classical laws, the hydrogen atom would be extremely 
short lived. Equation (5.6.8), by the way, gives an order-of-magnitude estimate of 
the time required by an atom to make a transition from an excited state to the 
ground state. 

The results of this section, especially Eq. (5.6.5), leave us with a curious paradox. 
Ignoring the fact that the classical laws of electrodynamics do not apply to a real 
hydrogen atom, we seem to find an inconsistency in our description of the classical 
atom. On the one hand, the fact that the electron feels the Coulomb field of the 
proton guarantees that the electron's motion will be circular and that its orbital 
radius will not change with time. (The most general orbit allowed by the Coulomb 
force is actually an ellipse, but the point remains that the size of this ellipse cannot 
change with time.) On the other hand, we know that the electron is accelerated, 
that it emits electromagnetic waves, that the radiation carries energy away from 
the system, and that the orbital radius cannot possibly stay constant. The paradox 
is therefore this: The electron's equations of motion do not account for the decrease 
in orbital radius; they unambiguously predict that the radius should be constant. 
Is classical electrodynamics internally inconsistent? 

Later on, in Chapter 6, we will show that the electron's radiation reaction is 
actually contained in Maxwell's theory. The effect is there, but it is well hidden, 
and it is a subtle matter to reveal it. 



5.7 Magnetic-dipole and electric-quadrupole 

radiation 

[The material presented in this section is also covered in Sec. 9.3 of Jackson's text.] 
Electric-dipole radiation corresponds to the leading-order approximation of the 
electromagnetic field in an expansion in powers of v c /c, where v c is a typical internal 
velocity of the source. In some cases, however, the dipole moment p either vanishes 
or does not depend on time, and the leading term is actually zero. In such cases, or 
when higher accuracy is required, we need to compute the next term in the expan- 
sion. This is our task in this section. We will see that at the next-to-leading-order, 
the wave-zone fields depend on the dipole moment vector p, the magnetic moment 
vector m, and the electric dipole moment tensor Q a b- We recall the definitions of 
these objects: 



p(t) = / p(t,x')x' d 3 x', 
Jv 



(5.7.1) 
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m(t) = If x' x j(t,x')d 3 x', (5.7.2) 
2 Jv 

Qab(t) = I p(t,x')(3x' a x' b -r' 2 S ab )d 3 x', (5.7.3) 
Jv 

where V is the volume that contains the charge and current distribution; this volume 
is bounded by the surface S. 

In the foregoing manipulations we will introduce the vector Q(t, x) defined by 

Qa(t,f) = Qab(t)r b , (5.7.4) 

where r b = x b /r are the components of the unit vector r = x/r, and in which 
summation over the repeated index b is understood. 



5. 7. 1 Wave-zone fields 

To make our calculations more efficient we first derive a few useful results that 
follow from the fact that in this section, we place ourselves exclusively in the wave 
zone. So throughout this discussion we shall assume that 

Notice that this condition is independent of any assumption that concerns v c : there 
is a wave zone whether or not the source is moving slowly. But if v c /c <C 1, then we 
also have that r c <C \ c /(2n) <C r, as was shown in Sec. 5.4.1. In this subsection we 
rely on the wave-zone condition of Eq. (5.7.5), but we do not rely on the slow-motion 
assumption; this is will be incorporated at a later stage. 

We can anticipate that the wave-zone potentials will have the form of a spherical 
wave. For example, we shall verify that the vector potential is given by 

AfcO-gS^ (5.7.6) 

in the wave zone, where u = t — r/c and to is a vector that will be determined. We 
will not need an expression for the scalar potential; as we shall see, in the wave zone 
E can be obtained directly from B, which is determined by the vector potential. 
Given the form of the vector potential, its curl is 

_ . wit -r/c) . d wit -r/c) 

V x A = ^— V x — i — = V^(Vr) x — . 

47r r 47r or r 

The gradient of r is the vector r, and the derivative of w(u)/r with respect to r is 
— c~ 1 w(u)/r, up to terms of order 1/r 2 that can be neglected in the wave zone. We 
therefore have 

B{t,x) = -^ fXw{u) (5.7.7) 
4-7T cr 

for the wave-zone magnetic field. 

Later we will find that w is not just a function of retarded time u, but that it 
depends also on f: w = w(u,f). This does not affect our result for the magnetic 
field, nor any other result derived in this subsection. To see this, let us calculate 
more carefully the x derivative of the vector w: 

_ <9uj dw _ „ 

" x w = ^-V^ii + — • V x r. 
ou or 

In the first term we substitute X7 x ii = — c _1 V ' x r = —c~ 1 r x and we recover our 
previous expression. We then recognize that V x r is of order 1/r and that the 
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second term can be neglected in the wave zone. We conclude that the dependence 
of w on r is invisible to V so long as all computations are carried out in the wave 
zone. 

To get the electric field we go back to one of Maxwell's equations, 



VxB = fj, j + 



1 dE 



Because the current density vanishes outside the volume V, we have that j = in 
the wave zone, and 



dE 

~dt 



c 2 V x B 



Mo£ 

47T 



V x 



r x w(u) 



By the same argument as the one just given, we can treat r as a constant vector 
when evaluating the curl, and we can also pull the 1/r factor in front of the derivative 
operator. This yields 



V x 

and we obtain 



f x w(u) 



f dw 

r dr 



= (Vr) x 
dE Mo ^ x x w(u)] 



= — r x 
c 



r x w{u) 



dt 4tt 

Integrating this with respect to t gives 



.-./ \ Mo \ r x w(u)] x r nr 

E(t,x) = - — = cB(t,x) x r 

4-7T r 



(5.7.8) 



for the wave-zone electric field. 

Using Eqs. (5.7.7) and (5.7.8) it is easy to show that the wave-zone Poynting 
vector is given by 



S{t , x) = ^B\H = ^^p\ 



(5.7.9) 



Evaluating this on a sphere of constant r gives 

Mo 



dP = S-da = 



or 



dP 



167T 2 C 

Mo 



r x w(u)\ dCt, 



16tt 2 c 



r x w 



Ml 1 



(5.7.10) 



where dtt = s\n6d6d(j) is an element of solid angle. Equation (5.7.10) gives the 
radiation's angular distribution, and here it becomes important to specify the de- 
pendence of w on the radial vector r. Notice that the fields and the radiation's 
angular profile depend only on f x w, and that they are insensitive to an eventual 
component of w along f. We will use this observation later, and discard from w 
any longitudinal component. 

From these results we infer, once more, that irrespective of the exact expression 
for the vector w(u, r), the fields are both transverse (orthogonal to r) and mutually 
orthogonal, that their amplitudes differ by a factor of c, that they are in phase, and 
that the field energy travels in the radial direction. To be more concrete we shall 
have to calculate the vector w. 
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(5.7.11) 



5.7.2 Charge-conservation identities 
Before we get to this we need to establish the useful identities 

f j(t,x')d 3 x'=p(t) 
Jv 

and 

f j(t,x'){f -x')d 3 x' = m(t) xr+-Q(t,r) + -r^ f p(t,x')r' 2 d 3 x' , (5.7.12) 
Jv 6 6 dt J v 

where p(t) is the electric dipole moment, m(t) the magnetic dipolc moment, and 
Q(t,f) the vector related to the electric quadrupole moment which we introduced 
in Eq. (5.7.3). 

Equation (5.7.11) was already derived near the end of Sec. 5.4.3; refer back 
to Eq. (5.4.12). Equation (5.7.12) is a generalization of Eq. (3.4.4), which covers 
static situations (and which was derived in Sec. 3.4), and we derive it using similar 
methods. 

Consider the z component of Eq. (5.7.12). We write the left-hand side as 



/ jz(t,x)xd 3 x, 
Jv 



in which we have changed the integration variables from x' to x and replaced the 
constant vector r by the arbitrary constant vector e. More explicitly, I z is given by 

I z = e x / xj z d 3 x + e y yj z d 3 x + e z zj z d 3 x. 
Jv Jv Jv 

Consider now the relation 

V • (xzj) = xzV • j + xj ■ Vz + zj ■ Vx 

dp 

= -xz— +XJ Z + ZJ X . 

Integrating both sides over V and using Gauss' theorem (together with the state- 
ment that j ■ da = on the bounding surface S) gives 

f xj z d 3 x = — f zj x d 3 x + -j" / pxz d 3 x. 

JV JV M Jy 

j yjz d 3 x = - j zj y d 3 x + ^ J pyz d 3 x 



We obtain 



in a similar way. Consider now 



V • (z 2 j) = z 2 V ■ j + 2zj ■ Vz = -z 



dp 
dt 



■2zj z ; 



integrating this over V leads to 

zj z d 3 x 



ld_ 
2dt 



pz 2 d 3 x. 



Let us now substitute these results into our previous expression for I z . We have 



1 



J xj z d 3 x - J zj x d 3 x + ^ J 



pxz d x 



1 

— e„ 
2 y 



yjz d 3 x - J zjy d 3 x + ^ J pyz d 3 x 



1 d f 2 ,3 

-e z — / pz a x 

2 dtJ v H 



^e x J (x x j) y d 3 x + ^e y J (x x j) x d 3 x 

/ pxz d 3 x + e y / pyz d 3 x + e z pz 2 d 3 x 
Jv Jv Jv 



ld_ 

2dl 



v 
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In this last expression, the first two terms are equal to —e x m y + e y m x = (m x e) z . 
The bracketed terms are related to Q z : Keeping in mind that we use e as a substitute 
for r, we have 

Qz{t> ^0 Qzx^-x "f" Qzy^y "f" Qzz&z 

= e x p(3xz) d 3 x + e y p{3yz) d 3 x + e z p(3z 2 - r 2 ) d 3 x 
Jv Jv Jv 

= 3[ ]-e z [ pr 2 d 3 x, 
Jv 

where [ ] stands for the bracketed terms in our previous expression for I z . Using 
these new results we finally obtain 

I z = (m x e) z + ~ (q z + e z J pr 2 d 3 x^j , 

which is the same statement (after the substitutions x — > x' and e — > f) as 
Eq (5.7.12). 



5. 7. 3 Vector potential in the wave zone 

We are now ready to calculate the vector w. We begin with the exact expression 
for the vector potential, 

A(t,x) = ^ [ 3{t-\x-x'\/c,x') d:ixl 



47T Jy \X — X 1 

In this we substitute the approximations 

1 1 



\x — x'\ r 

and 

t — \x — x'\/c ~ u + f ■ x'/c, 

where u = t — r/c is retarded time, which are appropriate in the wave zone. This 
gives 

A{t,x) = ^-- [ j(u + f ■x'/c,x')d 3 x' (5.7.13) 
47r r Jy 

for vector potential. Notice that this expression, though restricted to the wave zone, 
is otherwise general: it is not yet limited by a slow-motion assumption. 

We now incorporate the slow-motion approximation by Taylor-expanding the 
current density: 

1 dj 

j(u + f ■ x'/c) = j(u) + -(r • x') — (u) H . 

Relative to the leading term, the second term is smaller by a factor of order 
c~ 1 r c /t c = v c /c. This term was neglected in our discussion of electric-dipole radia- 
tion, but we shall now keep it. We will, however, neglect all remaining terms, which 
contribute at order (v c /c) 2 . 

Substituting the expansion for j into Eq. (5.7.13) gives 



47r r 



j(u, x') d 3 x' + l£-J j( u , x')(f ■ x') d 3 x' + 
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According to Eq. (5.7.11), the first integral is equal to p(u). And according to 
Eq. (5.7.12), the second integral is equal to 



m(u) xr + -Q(u,r) + -r — 



p(u, x')r' 2 d 3 x'. 



The quantity within the square brackets defines the vector w, and we have 



w = p(u) + 



rh(u) x f + ^Q(u, r) 



1 . d 2 

— r 

6c du 2 



p{u,x')r' 2 d 3 x'. 



The last term is proportional to r. But we have seen in Sec. 5.7.1 that the electric 
and magnetic fields, as well as the radiated power, depend only on r x w. This 
means that any component of w along f will not contribute to the wave-zone fields. 
We can therefore discard this last term in our expression for w. 
We have obtained 



w(u, r) = p(u) + 



rh(u) x r + ^Q(u, r) 



and in terms of this the wave-zone vector potential is given by 

po w(u,f) 



A(t,x) 



An r 



(5.7.14) 



(5.7.15) 



an expression that was already anticipated in Eq. (5.7.3). The first term on the 
right-hand side of Eq. (5.7.14) gives rise to electric-dipole radiation; this is the 
leading term in the expansion of the fields in powers of v c /c. The second and 
third terms give rise to magnetic-dipole radiation and electric-quadrupole radiation, 
respectively; these contribute at order v c /c beyond the leading term. 
To calculate the radiated power we shall need 

r x w = f x p(u) + -\rh(u) - \f ■ rh(u)]f \ + —r x Q (3) (?i, f), (5.7.16) 
c I J 6c 

where = d 3 Q/du 3 and we have used the vectorial identity r x (rh x f) = 

rh — (f ■ rh)r. 



5.7.4 Radiated power (magnetic-dipole radiation) 

Let us calculate the power radiated by a distribution of charges and currents for 
which 

p(u) = Q (3) (u,r) = 0. 

Such a source will emit magnetic-dipole radiation. 

Under the conditions of Eq. (5.7.17), Eq. (5.7.16) becomes 

1 



so that 



r x w 



r x w\ 



rh — (r ■ rh)f 
2 (r • rh) 2 



ml 



(5.7.17) 



which we rewrite in component form as 

\f x w\ 2 = ^rh a m b (S ab - r a r b ), 

where r a = x a /r are the components of the vector r, and where summation over 
both repeated indices is understood. Substituting this into Eq. (5.7.10) gives 



dP 



dn 16tt 2 c 3 



Mo 



m a fh b (6 ab - r a r b ), 
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and integrating over the angles yields 



P = ^^m a m b (S ab - (r a r b )), 



where 

If 1 

( r " r b> = J r andVi = -6 ab . 



These angular integrations are worked out in Sec. 5.7.7. 
Our final result for the radiated power is therefore 

-Pmagnotic-dipolo = ^^|m(u)| 2 . (5.7.18) 

We might compare this result with Eq. (5.4.21), 

p — ••( \\ 2 

clcctric-dipolc — g^^|^V^J| ■ 

In orders of magnitude we have that the electric dipole moment is p ~ (p^cX^c) ~ 
pr A and that p ~ pr^/t 2 . On the other hand, m <~ (j^c)^ 3 ) ~ jr\ ~ P v cri, and 
m ~ pv c r A c jt\. The ratio of powers is then 

magnctic-dipolc / ^1 \ I 



^ electric-dipole \ CP y \ C 

So for slowing-moving sources, the power emitted in magnetic-dipole radiation is 
typically smaller than the power emitted in electric-dipole radiation by a factor of 
order (v c /c) 2 <C 1. 

5.7.5 Radiated power (electric-quadrupole radiation) 

Let us now calculate the power radiated by a distribution of charges and currents 
for which 

p{u) = rh(u) = 0. (5.7.19) 

Such a source will emit electric-quadrupole radiation. 

Under the conditions of Eq. (5.7.19), Eq. (5.7.16) becomes 

r x w = — f x Q^, 
6c 

so that 

l f x = x ^ (3) ) ' xC ? (3) )' 

The inner product can be expressed as 

(r • r)(Q^ ■ Q< 3 >) - (f • Q (3) ) (r • Q (3) ) = |Q (3) | 2 - (r • Q (3) ) 2 , 
and we have 

\r^w\ 2 = ^- 2 Q^Qf\5 ab -r a r b ). 

But the vector depends on r, because by the definition of Eq. (5.7.4), Q a — 
Q ac r c - Writing also Q b = Q bd r dl we arrive at 



r x w\ 



36c 2 

Substituting this into Eq. (5.7.10) gives 



= ^QacQbd( 6 a-br c rd ~ r a r b r c r d ). 



dn ~ 16n 2 c 36c 2 Q ^ Q ™ ( 5abrcrd r ^b r cr d ) , 
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and integrating over the angles yields 



P = ^ c ^Q a S }Qf d \ta b {r c r d ) - (r a r h r c r d )). 



The angular integrals 

If 1 

( r c r d) = J r c r d dtt = -S cd 



{r a rbr c rd) = ^~ / r anr c rd dtt = (5 ab S cd + 8 ac 5 hd + S ad 8 bc ^j 



and 

lr~rur„ rj \ = „, „ 

47T J 15 

are evaluated in Sec. 5.7.7. Taking these results into account, we obtain 
T = Q^Qf}{5ab{r c r d )-{r a r b r c r d )) 

= iQa^QbJtabScd - ^Q ( a}Qf} (S ab S cd + 5 ac S bd + S ad 5 bc ) 

= iQ(3)Q(3) _ 1q(3)q(3) _ iQ{3)Q(3) _ Iq(3) Q (3) 
Q^ac^tac ^^ac^ac ^t!oot>66 ^g^ac^t-ca 



1 2 

3 ~ 15 



Vac Vac J 



where we have used the property that Q ab is a symmetric, tracefree tensor (so that 
Qca = Qac and Q aa = 0). 

Our final result for the radiated power is therefore 

-Pclcctric-quadrupolc = ^2Q^3 Qab ( u )Qab ( U ) ' (5.7.20) 

where summation over the two repeated indices is understood. In orders of mag- 
nitude, Q ab ~ (pr"c)( r c) ~ Fci and Qab ~ /° r c/*c ~ FcAc ~ jw c - We therefore 
have 

P / 7! ^ ^ 

clcctric-quadrupolc i ^< 



-Pelectric-dipole \ c 

for slowing- moving distributions, the power emitted in clectric-quadrupolc radiation 
is smaller than the power emitted in electric-dipole radiation by a factor of order 
(v c /c) 2 <C 1. We see also that electric-quadrupole radiation is of the same order of 
magnitude as magnetic-dipole radiation. 



5. 7. 6 Total radiated power 

In situations in which all three types of radiation contribute, the total power is the 
sum of individual contributions, 



6irc 



1 



1 



)(3), 



(5.7.21) 



There is no interference between the different terms in Eq. (5.7.16), because the 
angular profiles associated with each term are orthogonal. 
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5. 7. 7 Angular integrations 
In this subsection we establish the formulae 

(r„) = 0, (5.7.22) 

(r a r b ) = l -5 ab , (5.7.23) 

(r a r b r c ) = 0, (5.7.24) 

(r a r b r c r d ) = — (s ab 5 cd + S ac S bd + (5 a Ac), (5.7.25) 

where r a = x a /r are the components of the vector r, and (• • •) = (47r) _1 J (• • •) dCl 
is the average of the quantity (• • •) over a sphere of constant r. 

These relations can be derived by brute-force computation, by simply substi- 
tuting r = (sin 9 cos (f>, sin 9 sin 0, cos 6*) and integrating over (Kl = s\n9d9d<p. But 
there is an easier, prettier way that relies on arguments of symmetry. 

Take the zero result for (r a ). If these integrals were not zero, they would have to 
be equal to some vector a that does not depend on the angles 9 and 4> (because they 
have been integrated over) nor on r (because r appears nowhere on the left-hand 
side of the equation). But there is no such vector available, because there is no 
preferred direction in a homogeneous three-dimensional space. The right-hand side 
must therefore be zero. 

To arrive at the result for (r a r b ) we observe that the right-hand side must be 
a constant tensor that is symmetric in the pair of indices a and b. The only such 
tensor available is S ab , and we must therefore have {r a r b ) = XS ab} for some constant 
A. To determine A we simply take the trace of the left-hand side, 

5 ab (r a n) = (S ab r a r b ) = (1) = 1, 

and the trace of the right-hand side, 

XS ab S ab = 3A. 

Equating the results produces A = 1/3, and this agrees with the result of Eq. (5.7.23). 

The zero result for (r a r b r c ) follows from the fact that we do not have a vector 
that could be combined with S ab to form a symmetric, three-index tensor. 

To derive the result for (r a r b r c r d ) we note that the right-hand side must be a 
four-index tensor that is symmetric in all pairs of indices. We must form this tensor 
with two copies of S ab , and it is easy to see that it must be equal to X'(5 ab S cd + 
S ac S b d + Sadfihc)- To determine A' we evaluate the double trace of the left-hand side, 
which gives 1, and the double trace of the right-hand side, which gives 

X8 ab S cd (5 ab S cd + S ac S bd + S ad S bc ) = A'(3 x 3 + 3 + 3) = 15A'. 

Equating the results produces A' = 1/15, in agreement with Eq. (5.7.24). 



5.8 Pulsar spin-down 

As an application of magnetic-dipole radiation we consider the oblique-rotator model 
of a pulsar, a rotating neutron star that emits pulses of electromagnetic radiation 
at regular intervals. The most famous pulsar is located in the Crab nebula, and the 
pulsar is believed to be the source of energy for this supernova remnant. 

The Crab pulsar is observed to spin down: the frequency of the pulses decreases 
with time. If we associate the pulse frequency with the neutron star's spin frequency, 
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Figure 5.4: The neutron star rotates with angular velocity SI around the z axis. 
The magnetic moment m makes an angle of a with respect to this axis, and it 
rotates along with the star. 

then we can interpret this phenomenon in terms of the pulsar losing rotational 
energy, which is given by 

E rot = ^m 2 , (5.8.1) 

where I is the star's moment of inertia and S7 its angular velocity. If P denotes 
the pulses period (which is not to be confused with the power radiated by electro- 
magnetic waves), then fl — 2ir/P. If the neutron star is modeled as a solid sphere 
of uniform mass density, then I = |M_R 2 , where M is the pulsar's mass and R its 
radius. The loss of rotational energy translates into a decrease of SI, or an increase 
of P: 

E mt = mil = -(2tt) 2 /^. (5.8.2) 
P A 

For the Crab we have the observational values P~4x 10 -13 s/s and P ~ 0.03 s. 
If we take the neutron star to have a mass of 1.4 solar masses and a radius of 12 
km, then the rate of loss of rotational energy is 

E Iot ~ -7 x 10 31 J/s. 

This is comparable to what is required to power the Crab nebula. The energetics 
of the Crab can therefore be explained by the pulsar losing its rotational energy. 

What is responsible for this loss? We can imagine that the neutron star is carry- 
ing a magnetic field, and that the spinning motion of this field produces magnetic- 
dipole radiation. The energy carried off by the radiation will then come at the 
expense of the star's rotational energy. For this model to work we need the field's 
orientation to differ from the star's rotational axis — otherwise we would have a 
stationary situation analogous to the one analyzed in Sec. 3.3. This is the oblique- 
rotator model depicted in Fig. 5.4. 

We imagine that the neutron star maintains a magnetic dipole moment m that 
is oriented at an angle a with respect to the rotation axis. This vector rotates along 
with the star, also with an angular velocity S7. If the rotation axis is aligned with 
the z direction, we can write 

m(t) = toq ( sm a cos x + sin a sin fit y + cos a z) , (5.8.3) 
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where mo is the magnitude of the pulsar's magnetic moment. This can be related 
to the magnitude of the magnetic field on the surface of the neutron star. Because 
the star is located well within the near zone, we can use Eq. (3.4.2) to describe the 
surface field: 

_ fi 3(m ■ r)r - m 

where R is the star's radius. The field is maximum at the magnetic pole, where f 
is aligned with m. We therefore have 

x ~ 4~7r ~FP~ ' (5 ' 8 ' 4) 

a relationship between m , the star's radius, and the maximum surface magnetic 
field. 

The time-changing magnetic moment of Eq. (5.8.3) produces magnetic-dipolc 
radiation. This radiation takes energy away from the star, at a rate given by 
Eq. (5.7.18), 

#rad = ~ a o \m{u)\ , 

D7TC 1 1 

where u — t — r/c is retarded time. Substituting Eq. (5.8.3) gives 

Mo 



or 

2tt 



£ rad = ^(^o^ 2 sin a) , (5.8.5) 

2tt 

3/j.qc 3 



^rad - (B m ^n 2 R 3 sin a) (5.8.6) 



after involving Eq. (5.8.4). 

To produce an equality between £ ra( j and \E ro t\ — 7 x 10 31 J/s, we need a 
magnetic field such that 

B max sina ~ 5 x 10 s T. 

This is very large, but not unreasonable. A main-sequence star typically supports a 
magnetic field of 0.1 T. When the star collapses to form a neutron star, the magnetic 
field is frozen in and the magnetic flux 4ttR 2 B is conserved. Because R decreases 
by a factor of 10 5 during the collapse, B increases by a factor of 10 10 , and fields of 
the order of 10 8 T can easily be achieved for neutron stars. 

5.9 Problems 

1. A point electric dipole is rotating around the z axis at an angular velocity ui, 
so that 

p(t) — po [cos(wt) x + sin(wt) y] , 
where po is a constant. 

a) Calculate the wave-zone electric and magnetic fields of this rotating dipole. 

Express the fields in terms of the unit vectors 9 and <fi. 

b) Calculate (dP/dtt), the angular profile of the radiated power, averaged 

over a complete wave cycle. 

c) Calculate the total radiated power. How does it compare with the power 

radiated by an oscillating dipole? 

2. Two equal charges q are moving with constant angular velocity f2 on a circle 
of radius b. The charges are situated at diametrically opposite points on the 
circle. Calculate, to leading order in a slow-motion approximation, the total 
power radiated by this distribution of charges. 
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3. The current density inside a centre- fed, linear antenna is given by 

j(t, x) — I cos(wi) sm(kl — k\z\)8(x)8(y) z, 

where I is the current's peak amplitude, lu the oscillation frequency I the 
antenna's half-length, and k — oj/c. In this problem we do not impose the 
slow-motion approximation, so that kl is not assumed to be small. 

a) Show that in the wave zone (kr » 1), the vector potential can be expressed 

as 

A( )= *>«(M) 
4tt r 

and find an expression for the vector w(u, 9). Here, u = t—r/c is retarded 
time and 9 is the angle between x and the z axis. 

b) Calculate (dP/dil), the angular profile of the radiated power, averaged 

over a complete wave cycle. For the two special cases kl — n/2 (half- 
wave antenna) and kl = tt (full- wave antenna), provide parametric plots 
of this angular profile (like those displayed in Jackson's Figure 9.7). Com- 
pare these profiles with the sin 2 9 dependence obtained for electric-dipolc 
radiation. 

c) From the general expression for (dP/dQ,) derived in part b), obtain the 

short-antenna limit kl <C 1. Make sure that your result matches the one 
derived in the text. 

4. An oscillating current lit) — Ioe~ Jt sin(wi) flows in a circular loop of radius a 
placed in the equatorial plane of the coordinate system (r, 9, (p). The current 
is started at t = 0, and it decays exponentially due to dissipation within the 
loop; we assume that 7 <C w, so that the time scale for dissipation is much 
longer than the oscillation period. The current density is given by 

j(t, x) = ^ S(r - a)S(9 - it/2) 0. 

Calculate (P), the total power radiated by the loop, averaged over an os- 
cillation period. Calculate also the total amount of energy radiated by the 
loop. 

5. A nearly spherical surface described by 

r = R(t, 9) = a(l + e cos 9 cos uit), 

where a, uj, and e are constants, supports a nearly uniform distribution of 
charge. It is assumed that 6<1, and the charge density is given by 

p(t,x) = ^5{r-R(t,6)), 

where q is the total charge. 

Calculate (P), the electromagnetic power emitted by this distribution of 
charge, averaged over a complete wave cycle. You may work to leading order 
in the slow-motion approximation. You may also work to leading order in e. 

6. A ring of radius b is placed in the x-y plane, and the origin of the coordinate 
system is placed at the centre of the ring. Initially, the ring is at rest and it 
supports a linear charge density proportional to sin0, where <j) is the angle 
from the a;-axis. The ring is then made to rotate with a uniform angular 
velocity 0, so that the linear charge density becomes proportional to sin((/> — 
£lt). The ring's volume charge density is described by 

p(t, r, 9, 4>) = ^S(r - b)6(9 - f ) shu> - fit). 
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(a) Calculate, to leading order in a slow-motion approximation, the total 
power radiated by the ring. At what angle 9 is most of the radiation 
emitted? 

(b) Find an exact expression for $(t, 0), the scalar potential at the centre of 
the ring. 

(c) Find an exact expression for A(t, 0), the vector potential at the centre of 
the ring. 

7. Three charges are located along the z axis: a charge +2q stays at the origin 
z = 0, a charge — q is at z = +a cos ujt, and the final charge — q is at z = 
— acoswi; a and w are constants. It is assumed that auo <C 2nc and that the 
charges simply go through one another when they meet at z = 0. 

For this situation calculate (dP/dil), the angular profile of the radiated power, 
averaged over a complete wave cycle. (Express this in terms of the angle 9 
from the z axis.) Calculate also (P), the total radiated power. What is the 
frequency of the electromagnetic waves emitted by this system of charges? 

8. A current / is suddenly established in an infinite wire that extends along the 
z axis. The current density is given by 

j(t,x)=I6(t)S(x)S(y)z, 

where 9(t) is the step function. Because V • j = 0, it is consistent to set 
p(t,x)=0. 

Calculate, without approximations, the vector potential A{t, x) associated 
with this distribution of current. From your answer you should be able to 
verify that at late times, the electric and magnetic fields that are produced 
by the distribution of current have the components 

_^2J Mo "My o Mo 21a; 

z ~ 4tt t ' 4^ x 2 + y 2 ' v ~ 4tt x 2 + y 2 ' 

The electric field eventually disappears, and the magnetic field settles down 
to the value it would have for an unchanging current. 
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Chapter 6 
Electrodynamics of point 

charges 



6.1 Lorentz transformations 

[The material presented in this section is also covered in Sees. 11.3, 11.9, and 11.10 
of Jackson's text.] 

In this final chapter we will look at the electrodynamics of charge and current 
distributions for which the characteristic speed v c is close to the speed of light. To 
keep the discussion concrete we will restrict it to the electromagnetic field produced 
by a single point particle in arbitrary (relativistic) motion. To set the stage we 
begin with a review of the Lorentz transformations of special relativity. 

We consider a reference frame S' that moves with constant speed v relative to 
another frame S; we take the motion to take place in the x direction, so that the 
frame's velocity vector is v = vx. The coordinates (t', x' , y', z') of an event seen 
in 5" are related to the coordinates (t, x, y, z) of the same event seen in S by the 
well-known relations 



x' = i(x-vt), y' = y, z' = z, t'=7^t-^a:^. 



(6.1.1) 



where; 



7= h 2/2 - (6-1.2) 



These are the Lorentz transformations between the two reference frames. The re- 
versed transformation is obtained by reversing the sign of v: 



x = 1 {x' + vt l ), y = y', z = z', i = 7 ^'+Ja/). 



(6.1.3) 



We now would like to derive how the electric and magnetic fields transform under 
such a change of reference frame. 

We note first that Lorentz transformations alter the differential operators. For 
example, applying the chain rule on an arbitrary function ip gives 

dip dip dt dip dx dip dy dip dz 
dt' ' ~ ~dtdt' " 4 'fadt' ~dy~dt' ' ~dz~dt n 



or 



dt' \dt dx) 
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after using Eq. (6.1.3). Doing a similar calculation for all other partial derivatives, 
we obtain 

d ( d v d\ d d d d d ( d d 

dx' ^\dx c 2 dt )' dy' dy ' dz' dz ' dt' \ dt dx / 

(6.1.4) 

It follows from this that the wave operator is a Lorentz invariant: 

= i ay d 2 j> d 2 j> o> 2 y> 

V c 2 dt ,i + dx ,2 + dy a + dz a 

i 2 ( d d\fd d \ . 



+ 7 ' [dx + c 2 dt ) [dx + c 2 dt ) ^ + dy 2 + dz 2 



d v d \ ( d v d \ d ip d ip 



c 2 dt 2 dx 2 dy 2 dz 2 
Dip. 



In other words, 



□' = □. (6.1.5) 



To figure out how the electric and magnetic fields transform under a change of 
reference frame, we first work out how p and j should transform. This will tell us 
how the potentials $ and A transform, and from this we will be able to get to the 
fields. 

It is sufficient to consider a simple special case: a point charge q moving in the 
x direction with a speed v (the same speed that enters in the Lorentz transforma- 
tion). As seen from S we have p — qS(x — vt)S(y)5(z), j x — qvS(x — vt)S(y)6(z), 
and j y = j z = 0. As seen from S', however, the particle is not moving, and 
p' = qS(x')S(y')S(z'), while j' x = j' y = j' z = 0. Wc can expect that the Lorentz 
transformation will involve a linear mixture of p and j x , and that the other compo- 
nents of the current density will be left alone. Rewriting p' in terms of the unprimcd 
coordinates gives 

p' = q5[ 1 {x + vt)]5{y)8{z) 

= ^5{x-vt)5{y)5{z). 

7 

This expression for p' involves a factor of 1 /-f instead of the expected 7, and j x has 
not yet made an appearance. But we can write it in the equivalent form 

.2 l„2 



g 7 (l - v 2 /c 2 )S(x - vt)S(y)S(z) 



= 7 



v 



qS(x — vt)S(y)S(z) -qvS(x — vt)S(y)S(z) 



v 



and this is the expected expression of a Lorentz transformation. This shows that 
p transforms as t, and wc further expect that j x will transform as x, that is, j' x — 
l{jx ~ vp). This is confirmed by the fact that the right-hand side is zero for the 
situation considered here. We conclude that under a change of reference frame, the 
charge and current densities transform as 

fx = lUx - vp), j' y = j y , j' z =j z , p' = l(p~ ^jx^j ■ (6.1.6) 



6.1 Lorentz transformations 
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While the derivation of these transformation rules was restricted to a very simple 
situation, it is possible to show that they are in fact completely general. 

To see how the potentials transform we recall the wave equations □(<f>/c 2 ) = 
— p/(eoc 2 ) = — pop and OA — — poj from Sec. 5.1. Because □ is an invariant, these 
equations tell us that $/c 2 must transform as p, and that A must transform as j. 
Thus, 

A' x = 1 [a x -^, A' y = A v , A' Z =A X , <&' = 7 ($ - vA x ) (6.1.7) 

are the Lorentz transformations for the scalar and vector potentials. These are 
designed to leave the form of the wave equations unchanged under a change of 
reference frame. You should check that the Lorenz gauge condition, 

„ . 1 9$ « 
V-A+^— =0, 

& at 

also takes the same form in S' . All the equations satisfied by the potentials are 
therefore preserved under a Lorentz transformation. This implies that Maxwell's 
equations themselves are invariant under a change of reference frame. 

We can now figure out how the fields transform. We simply write down how the 
fields in S' are related to the potentials, 

E' = - 9 ^--V'&, £?' = V' x A', 
dt' 

and we substitute Eqs. (6.1.4) and (6.1.7). For example, 

*' ~ ^(l^s)(^-?*)-^5 + ?i)(*-^) 

dA x <9$ 
dt dx 
= E x . 



Or 



\ dt dy ) \ dx dy 
= 1 (E y -vB z ). 



Another example is 



Yet another is 



B , = 9A^_dA l = dA ± _dAy =Bx 
dy' dz' dy dz 



dAL dA' 



dz' dx' 

d ( v \ ( d v d\ . 
^ z { A *-sV-\o- X + Sdt) A * 

jdA x dA z \ i, / dA z dt> 

\ dz dx J c 2 I dt dz 
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The complete summary of these results is that the electric field transforms as 

E' X =E X , E' y = 7 (E y -vB z ), E' z = 1 {E z +vB y ) (6.1.8) 
under a change of reference frame, while the magnetic field transforms as 

B' X = B X , B ' y = 1 (B y + ^E^ B' z =j(B z -^E y y (6.1.9) 

Two observers in relative motion will therefore disagree as to the value of the electric 
and magnetic fields. 

There are interesting special cases to the foregoing results. Suppose first that 
B = in some reference frame 5*. In another reference frame S' we have E' x = E x , 
E' y = ~fE y , E' z = 1 E Z1 B' x =0,B' y = l{v/c 2 )E z = (v/c 2 )E' z = -(v x E') y /c 2 , and 
B' z = -j(v/c 2 )E y = -{v/c 2 )E' y = -(« x E') z /c 2 . In other words, 

B = -> B' = -\vxE'. (6.1.10) 

c 2 

A pure electric field in 5* is therefore perceived as a mixture of electric and magnetic 
fields in S' . Suppose next that E = in S. Then B' x = B x , B' y = ^B y , B' z = jB z , 
E' x = 0, E' y = - 1V B Z = -vB' z = (vx B') y , and E' z = ~fvB y = vB' y = (v x B') z . 
In other words, 

E = -» E' = vxB'. (6.1.11) 

A pure magnetic field in S is therefore perceived as a mixture of electric and mag- 
netic fields in S'. We conclude that E and B are observer-dependent fields; only 
the electromagnetic field as a whole has an observer-independent existence. 



6.2 Fields of a uniformly moving charge 

[The material presented in this section is also covered in Sec. 11.10 of Jackson's 
text.] 

We can use the Lorentz transformations of Eqs. (6.1.8) and (6.1.9) to calculate 
the electric and magnetic fields of a point charge q that moves with constant speed 
v in the x direction. The charge's position vector is 

r(t) = (vt,b,0), (6.2.1) 

where b is the charge's displacement above the x axis. We want E and B as 
measured at O, the origin of the reference frame S. We shall first do the calculation 
in S' , which moves along with the particle. 

As seen in the reference frame S", the charge's position vector is r' = (0,6,0), 
but the observer at O has a position r' Q = (—vt',0,0), because relative to S', O 
moves in the negative x direction with speed v. In 5" there is only an electric field 
E', and its value at O is given by the standard expression 



E , = g r - r o 



47re \r' - *«'-l 3 " 



We have r' — r' Q = (vt' , b, 0), \r' — r' \ = \Jb 2 + vH' 2 , and the primed components 
of the electric field are 

pi = g V l_ pi = Q b p/ = n ffi OOA 

47re (6 2 +w 2 t' 2 ) 3 / 2 ' y 4^e (b 2 + v 2 t' 2 f/ 2 ' z K ' ' 

As measured in S' , the electric field lines are distributed isotropically around r' , 
the charge's position. 
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Figure 6.1: Electric field of a uniformly moving charge, as measured by an observer 
at the origin of a fixed reference frame. The x component of the field goes through 
zero at t = 0, while the y component is always negative and achieves its maximum 
(in magnitude) at that time. 

We now transform the field to the frame of reference S. The coordinates of O 
are x = y = z = 0, and the time coordinates are therefore related by t' = "ft. The 
unprimed fields are obtained by inverting Eqs. (6.1.8) and (6.1.9), which is done by 
reversing the sign of v. Because B' = we have E x = E' x , E y = jE' v , E z = E' z = 0, 
B X =B' X = 0, B y = -{~fv/c 2 )E' z = 0, and B z = (jv/c 2 )E' y . We therefore obtain 

rp = I l vt F = 9 7fr F =n 

47Te (& 2 + 7 2w2 t 2)3/2 ' V 4^ ( fo 2 + ^2^)3/2 ' 

(6.2.3) 

and 

B x =0, B y =0, B Z = ^E V (6.2.4) 

for the fields measured at O, in the S frame. Plots of E x and E y as functions of 
time are shown in Fig. 6.1. 

The fields can also be expressed in terms of the charge's position vector, r(t) = 
(vt, b, 0). For this purpose we let ip be the angle between the vectors r and v. This 
is defined by 

co^= V -^, (6.2.5) 
vr(t) 

where r = |r| = \/b 2 + v 2 t 2 . We have v ■ r = v 2 t, so cos-0 = vtj\Jb 2 + v 2 t 2 , which 
implies sin-0 = b/y/b 2 + v 2 t 2 . We also have 

„2/j,2 /„,2 , „.2.2\ 



2 f .2 , „,2 .2 V .2 



2 2 A V 2 . 2 A 

]r II g sin if; J, 
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and this can be substituted into Eqs. (6.2.3). The end result is 

E = -T 1 1 ( 6 - 2 - 6 ) 

47T60 7 2[!_ ( w / c )2 sin 2 ^] 3 / 2 r3 

for the electric field at O, and 

B= ^rvx E (6.2.7) 
c 2 

for the magnetic field. 

We see from Eq. (6.2.6) that E is always directed along — r(t), a vector that 
points from the charge to O. The electric field therefore points away from the charge, 
and this means that the electric field lines originate at the position of the moving 
charge. The field lines, however, are not distributed isotropically around the charge; 
this is revealed by the presence of sin 2 tp in the expression for the electric field. When 
ip = 0, that is, when r is directed along v, we have that [1 — (v/c) 2 sin 2 ip]~ 3 / 2 = 1, 
and the electric field is reduced by a relativistic factor of -f~ 2 < 1. But when instead 
ip = 7r/2, so that r is orthogonal to v, then [1 — (v/c) 2 sin 2 ?p}~ 3 / 2 = 7 3 > 1, and the 
field is amplified by a factor 7 > 1. The electric field, as seen in S, is strongest in 
the directions transverse to the motion, and weakest in the longitudinal direction. 

The results derived in this section are at once remarkable and unremarkable. It 
is unremarkable that the field lines should emanate from the position of the charge 
— we have the same picture in the charge's rest frame. And it is unremarkable that 
the field lines should be distorted by a Lorentz contraction in the direction of the 
motion. But it is remarkable that the field registered at O does not involve a time 
delay associated with the propagation of a signal from charge to observer. 

Allowing for such a delay would suggest that the field at time t should point not 
toward the charge at its current position r(t), but instead toward an earlier position 
r(t'). In this picture, a signal originates at the charge when it was at the earlier 
position, and it propagates to O in a time At = t — t' determined by the condition 
that the signal travels at the speed of light; since the distance traveled is 
we must have t — t' = \r(t')\/c, which in principle can be solved for t'. During 
the time in which the signal propagates, the particle moves to its current position, 
but the field registered at O should not yet be aware of this; it should know only 
of the earlier position, and E should therefore point in the direction of r(t'). But 
this is not what we have found: the field does point in the direction of r(t), the 
current position of the charge. Our picture, according to which information should 
be delayed, is curiously wrong. We will return to this issue later on, and investigate 
it more fully. 



6.3 Fields of an arbitrarily moving charge 

[The material presented in this section is also covered in Sees. 6.5 and 14.1 of 
Jackson's text.] 

We now allow the point charge q to move on an arbitrary trajectory r(t). We 
let v(t) = dr/dt be its velocity vector, and a(t) = dv/dt its acceleration vector. 
The charge density is 

p(t,x)=qS(x-r(t)), (6.3.1) 

and the current density is 

j(t,x) = qv{t)5(x-r(t)). (6.3.2) 
To obtain the fields produced by this particle we shall first calculate the potentials. 
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6.3.1 Potentials 



To calculate the potentials associated with the densities of Eqs. (6.3.1) and (6.3.2) 
we employ a technique that was first put to use in Sees. 5.2 and 5.3. For example, 
we express the vector potential as 



A{t,x) 



47T 



Jj(t',x')- 



IX — X' 



in terms of the retarded Green's function obtained in Sec. 1.9. Substituting Eq. (6.3.2) 
yields 



A(t,x)= q ^jv(t')S(x f -r(t')) 
and integrating over d 3 x' returns 



^(t-t'-lx-xVc)^, 



IX — X' 



A(t,x) 



47T 



The same manipulations produce 

q 



S(t-t'-\x-r(t')\/c) 
V(t) ' |aj-r(f)| 

5(t-t'- \x-r{t')\/c) 



47T60 



dt' 



(6.3.3) 



(6.3.4) 



\x-r(t')\ 

for the scalar potential. At this stage it is convenient to introduce the quantities 

H(f) 



R{t') =x- r(t'), R{t') = \x- r(t') 



R(t') 



R(t') 



(6.3.5) 



In terms of these the argument of the (5-function is t — t' — R(t')/c. 

To evaluate the integrals of Eqs. (6.3.3) and (6.3.4) we let s = t' + R(t')/c — t 
be the new variable of integration. We have ds/dt' = 1 + c^dR/ 'dt' , and it is easy 
to verify that dR/dt' = R -v/R = —R ■ v. With this we obtain ds/dt 1 = k, where 



k= 1 -v(t') -R(t')/c. 
The integral for the vector potential is then 

qfio f vS(s) dt' 



(6.3.6) 



A = 



A-k 



R ds 



ds 



qjio f vS(s) 
4tt ./ kR 

4tt 



ds 



kR 



s=0 



Similarly, 



47re n 



kR 



s=0 



The condition s — gives R(t')/c = t — t', or 

\x-r(t')\=c(t-t'), 



(6.3.7) 



and this implicit equation can in principle be solved for t' . The meaning of Eq. (6.3.7) 
is clear: The signal received at position x and time t depends on the state of motion 
of the particle at an earlier time t'; the delay (t — t') is set by the time required by 
light to travel the distance R(t'). We shall call t' the retarded time corresponding 
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to the spacetime point (t, x), and Eq. (6.3.7) will be referred to as the retarded 
condition. 

We shall write the potentials as 



and 



A(t,x) 



Q 

47re 



47T 



1 



kR 



(6.3.8) 



(6.3.9) 



indicating that the quantities within the square brackets must be evaluated at the 
time t' determined by the retarded condition of Eq. (6.3.7). These expressions are 
known as the Lienard- Wiechert potentials. Apart from the presence of k in the de- 
nominator, the potentials take their familiar form from time-independent situations, 
and they implement the notion that the potentials here and now should depend on 
the retarded position and velocity of the particle. Given that the retardation effect 
is present in the Lienard- Wiechert potentials, it is all the more curious that the 
effect was not seen in the preceding section. (Stay tuned.) 
The factor of 1/k accounts for the fact that 



/ 



p(t- \x-x'\/c, x')d 3 x' 



is not equal to the charge q. While we do have J p(t,x') d 3 x' = q, in the previous 
expression we sample the charge density at different times while performing the 
integration, and the result is not q. Instead, 



J P (t -\x- x'\/c, x') d 3 x' = J p(t', x')S(t - t' - \x - x'\/c) d 3 x'dt 

= q J S(t-t' -\x-r{t')\/c)dt' 



and this result clarifies the origin of n = 1 — v ■ R/c. 



6.3.2 Fields 

Having obtained the potentials, we can calculate the fields by straightforward dif- 
ferentiation, 

OA 

E = -— V$, B = V x A. 

dt 

The dependence of A on t comes from its explicit dependence on t' and the retarded 
condition \x — r(t')\ = c(t — t'). Its dependence on x comes from the explicit depen- 
dence of R and k on position, but there is also an implicit dependence contributed 
by the retarded condition. This mixture of explicit and implicit dependence com- 
plicates the task of taking derivatives. But we can simplify things by going back to 
earlier forms of Eqs. (6.3.3) and (6.3.4), 

*^=4k Jw) 6iS)dt ' (6 ' 3 - 10) 

and 
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where s = t' +R(t')/c—t. Their advantage is that t' inside the integrand is no longer 
restricted, so that the dependence of the potentials on t and x is fully explicit. It will 
then be straightforward (if tedious) to differentiate first, and integrate next (thereby 
setting t' by the retarded condition). Before we proceed we recall the definitions of 
Eqs. (6.3.5) and (6.3.6): R(t') =x- r(t') and k = 1 - v{t') ■ R(t')/c = ds/dt'. 
We begin with the scalar potential. We have 



47Te 



6(a) + ^S'(s)Vs 



dt', 



and since s depends on x through R, we can replace Vs by Vii/c. This gives 



g 

47T60 
g 

47T£ 



/ 



v( 

1, 



(>> 

5(s 



+ —S'(a) 



/[Ms) 



cR 

VR .,. 
f ^R & {S 

d (VR\ 

ds \ ckR J 



dt' 



ds 



S(s) ds, 



where we have integrated the second term by parts. The derivative with respect to 
s can now be expressed in terms of d/dt', and we finally obtain 



Q 



47T60 



i) 



1 d fVR\ 



CK 



dt 1 \kR J 



S(s) ds. 



Leaving this result alone for the time being, we move on and calculate dA/dt. We 
have 



8A 

~dt 



47re 6 
q 1 



ds 

at 



dt' 



47T£o C 2 

Li/ 

T£0 C 2 J 



S'(s) ds 



47re 

q 1 
47reo c 2 



nR 

-(-) 
ds \ kR J 



1 d 



(5(s) cZs 
(5(s) ds. 



In the first step we used the fact that the only t dependence is contained in s. In 
the second step we substituted ds/dt = — 1. In the third step we integrated by 
parts. And in the fourth step we replaced the s-derivative by a t'-derivative. 
Collecting these results gives us an expression for the electric field: 



E 



47re 

q 



.1, 

K 
1. 



■a 



1 d (VR\ 

ck dt' \ kR J 
1 d fVR-v/c 
ck dt' \ 



J_d_ 
c Y ~K~dt' 



V 

kRJ 



S(s) ds 



kR 



S(s) ds. 



This can be integrated at once, and we obtain 



E 



q 



47re 



K \ 



+ 



1 d fVR-v/c 
ck dt' \ kR 



This expression badly needs to be simplified. It is easy to establish that Vi? = R, 
so that Vi? -1 = —R/R 2 , and we already have obtained dR/dt' = —v ■ R in the 
preceding subsection. We now set to work on 



d_ 
dt' 



VR-v/c 
kR 



1 / dR a\ R vie d / \ 

= -'— ' dt\ KR h 



kR \ dt' c 



K 2 R 2 
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where a = dv/dt' is the particle's acceleration. We have dR/dt' = d(R/ R) / dt' = 
-v/R + (R/R 2 )(y ■ R) = -[v - (v ■ R)R]/R, and also d(nR)/dt' = d(R - v ■ 
R/c)/dt' = —v R (a ■ R — v ■ v)/c = —v R a R/c + v 2 /c. Plugging these 
results back into our previous expression gives 



cdf 



/VR-v/c\ 
\ kR ) 


1 

kR 2 


v/c 




R - 


v/c 




K 2 


R 2 




1 

kR 2 


v/c 



v/c - (v ■ R/c)R 



a 

kRc 2 



(-v ■ R/c + v 2 /c 2 ^ + 



R — v/c a R 

K 2 R 2 c 2 ^ 



+ k - v/c (l 

k 2 R 2 V 



v 2 /c 2 ) 



1 



k 2 Rc 2 

R-v/c 

kR 2 
1 



k 2 Rc 2 



(R — v/c)(a ■ R) — Ka 

R R v/c f 2 2 \ 
~W + k 2 R 2 \ /C ) ~ kR 2 

(R — v/c)(a ■ R) — na . 



R-v/c 



The first and fourth terms cancel out, and the expression within the square brackets 
is Rx [(R-v/c) xa} = (R- a)(R-v/c) - R (R-v/c)a = (R- a)(R- v/c) - na. 
This gives 



1 d fVR-v/c 
cdfy kJ1 



R R 

+ 



v/c Rx [(R - v/c) x a] 



R 2 ~/ 2 k 2 R 2 



K 2 C 2 R 



Substituting all this into our previous result for the electric field returns 



E(t,x) 



47re n 



R-v/c Rx [(R-v/c) x a] 
j 2 k 3 R 2 + k 3 c 2 R 



(6.3.12) 



ret 



our final expression for the electric field of a moving point charge. The electric 
field is naturally decomposed into a "velocity field" that decays as 1/R 2 and an 
"acceleration field" that decays as 1/R. The velocity field represents that part of 
the electric field that stays bound to the particle as it moves; it is the particle's 
generalized Coulomb field. The acceleration field, on the other hand, represents 
electromagnetic radiation that propagates independently of the charge. 

To obtain the magnetic field we take the curl of Eq. (6.3.11), which gives 



V x A = 



1 



_/ 1 \ P , . Vs x v . , 
V I — x vS(s) H — 5'(s) 



\R 



47T£o C 2 



47TC0 C 

q 1 



R 



VRXV 5>( S ) 



47re c 2 

q 



1, 



ckR 
d ( VR x v 



x v 



so that 



B = 



47T£o C 



q 1 

47T£o C 2 



ds \ ckR 
1 d /VR x v 
ck dt' \ kR 



dt' 

ds 
S(s) ds 
S(s) ds, 



1 d (VR 



en dt' V kR 



We now work on simplifying this expression. We have V(l/R) xv = —(R xv)/R 2 , 
and d(R x v/nR)/dt' can be expressed as 



d (R-v/c) xv 



dt' 



kR 



d_ 
dt' 



R-v/c 
~kR 



R-v/c 

x v H — x a, 

nR 
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which allows us to use results that were generated during the calculation of the 
electric field. We have 



1 d (R-v/c) xv 



cdt> 



kR 



R R—v/c (R- a)(R — v/c) — na 



+ 



R 2 j 2 k 2 R 2 
R-v/c 



k 2 Rc 2 



X V 



ckR 
Rx v 



x a 



+ 



R 2 
Rx a 



+ 



Rxv (R- a)Rx v axv 



ckR 
Rxv 



+ 



J 2 K 2 R 2 

axv 
' kRc 2 
Rxv 



+ 



k 2 Rc 2 



kRc 2 



R 2 j 2 k 2 R 2 



(R ■ a)R xv Rx a 

+ " s4"5 + 



K 2 R C 2 



ckR 



Collecting these results, we arrive at 



or 



B 



B{t,x) = - 



1 



47reo c 2 



47re c 2 



Rxv (R-a)Rxv Rxa 

j 2 k 3 R 2 k 3 Rc 2 CK 2 R 



Rxv (R- a)R x v/c + kR x a 



ret 



J 2 K 3 R 2 



+ 



CK 



3 R 



(6.3.13) 



A little more algebra reveals that this can also be expressed as 



c 1 J 



x E. 



(6.3.14) 



The magnetic field also is decomposed into a velocity field that stays bound to the 
particle, and an acceleration field that represents electromagnetic radiation. 

Our general results for the electric and magnetic fields allow us to resolve the 
issue that was first brought up near the end of Sec. 6.2: where does the electric 
field point? Does it point toward r(t'), the retarded position of the particle, as it 
was argued it might do? Or does it point toward r(t), its current position, as we 
actually found in Sec. 6.2? The answer is that in general, the electric field follows 
neither option. 

Consider our expression (6.3.12) for the bound electric field. (We exclude the 
radiative field from this discussion, as we expect it to behave largely independently 
of the charge.) It shows that the field points in the direction opposite to the vector 
— R + Rv/c evaluated at the retarded time t' . If the observation point is at the 
origin of the coordinates (x = 0), then R= —r and the bound field points in the 
direction opposite to 

r(t>) + (At)v(t'), 
where At is determined by the retarded condition 

At = t-t' = r(t')/c. 

This result is interesting: the field points not in the direction opposite to the charge's 
retarded position r(i'), nor in the direction opposite to the charge's current position 
r(t), but in the direction opposite to the charge's anticipated position r a = r{t') + 
(At)v(t'). This is where the charge would be after a time At if it were moving 
with a constant velocity v(t'); the time interval corresponds to the time required 
for a light signal to travel from the charge's retarded position to the observation 
point. (The situation is described in Fig. 6.2.) So the field brings information to the 
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retarded position 



anticipated position 




current position 



Figure 6.2: Retarded position, current position, and anticipated position of an 
accelerated charge. 



observer that does indeed originate from the retarded position, but the field orients 
itself as though it were anticipating the motion of the particle during the time the 
information propagates from the charge to the observer. Because the extrapolation 
is based only on partial information (the values of r and v at the retarded time, 
and nothing more), it is not fully accurate, and the field misses the actual position 
of the particle. 

Notice that uniform motion is a special case: when the velocity is constant, the 
charge's anticipated position r a coincides with its current position r(t), and the 
field "miraculously" orients itself in this direction. This explains the result found 
in Sec. 6.2: the information needed time to reach the observer, but it was sufficient 
to correctly anticipate the actual motion of the charge. 



6.3.3 Uniform motion 

Let us verify that the general expression of Eq. (6.3.12) reduces to Eq. (6.2.3) in 
the case of uniform motion in the x direction. 

We set r(t) = (vt,b,0), r{t') = (vt',b,0) = r', and v(t) = v(t') = vx. We 
once more take the observation point to be at x = 0. Because a = 0, Eq. (6.3.12) 
reduces to 

q f' + v/c q r' + r'v/c 

47re 7 2 K 3 r' 2 47re j 2 K 3 r' 3 ' 

or 



q ^(r' + r'v/c) 
47T£o (7«t') 3 



We recall that k = 1 + v ■ f'/c. 
The retarded condition is 



(6.3.15) 



r' = v/fc2 + {vt'f = c(t - t'), 
and it gives, after squaring, 

(1 - v 2 /c 2 )t' 2 - 2tt' + t 2 - b 2 /c 2 = 0. (6.3.16) 



This is a quadratic equation for t', which we do not need to solve. 
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The numerator of Eq. (6.3.15) involves the vector r' + r'v/c. Its x component 
is vt' + r'v/c = vt' + c(t — t')v/c = vt, and this vector therefore coincides with r. 
In the denominator we have (-jur 1 ) 2 raised to the power 3/2. This is 

(7k/) 2 = 7 2 (r' + v ■ r'/c) 2 

= 7 2 (r' + v 2 t'/ c ) 2 
= ^^(t-t' + vh'/c 2 ) 2 
= ^(t-t'/j 2 ) 2 

= c 2 ( 1 2 t 2 -2tt' + t' 2 / 1 2 ). 

We now use Eq. (6.3.16) to eliminate t' and obtain 

( 1K r'f - c 2 [b 2 /c 2 + 7 2 t 2 (l-l/ 7 2 )] 
= & 2 + 7 Vt 2 . 

Substituting these results into Eq. (6.3.15) gives 

F = 1 7^ft) (a<^7) 

4TT€ (62+^^2)3/2' I?-*- 1 ') 

the same statement as Eq. (6.2.3). 



6.3.4 Summary 



The scalar and vector potentials of a point charge in arbitrary motion described by 
the position vector r(t) are given by the Lienard-Wiechert expressions 



and 



A(t,x) 



Q 



47re 



47T 



1 



kR 



ret 



where R(t') =x- r(t'), R = \R\, R = R/R, 

K = 1 -v(t') ■ R(t')/c, 

and where t' is determined by the retarded condition 

\x-r(t')\=c{t-t'). 

The electric field is 

q 



and the magnetic field is 



R v/c Rx [(R-v/c) x a] 



7 2 K 3 i? 2 



K 3 C 2 R 



B = \R] xE. 

n L J ret 



(6.3.18) 



(6.3.19) 



(6.3.20) 



(6.3.21) 



(6.3.22) 



(6.3.23) 



Both fields are naturally decomposed into velocity (bound) fields that decay as 
1/i? 2 , and acceleration (radiative) fields that decay as 1/R. 
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current position 



Figure 6.3: Sphere S' centered at the retarded position r(t') of the accelerated 
charge. The sphere has a constant radius R(t') = \x — r(t')\. 

6.4 Radiation from an accelerated charge 

[The material presented in this section is also covered in Sees. 14.2, 14.3, and 14.4 
of Jackson's text.] 

At large distances from the particle — in the wave zone — we can neglect the 
velocity fields and approximate the electric field by 



E(t,x) = -«- 
47re 

and the magnetic field by 



R x [(R-v/c) x a] 



K 3 C 2 R 



(6.4.1) 



We recall that R = x - r(t'), R = \R\, R = R/R, K = 1 - v(t') ■ R(t')/c, and 
that t' is determined by the retarded condition t — t' = R(t')/c. In this section we 
use these expressions to calculate how much power is radiated by the accelerated 
charge. We also determine the angular distribution of the radiated power. 



6.4-1 Angular profile of radiated power 

Because our expressions for the fields refer directly to the charge's retarded position, 
it is a relatively simple task to calculate the power that crosses a sphere S' centered 
at this position (see Fig. 6.3); it would be much more complicated to adopt any 
other reference point. The sphere S" has a constant radius R(t') and its element of 
surface area is da' — R 2 (t')R(t') rffi', where dfl' is an element of solid angle on the 
sphere. 

The energy crossing an element of S' per unit time t is, by definition, S(t') ■ da' , 
where S = ExB /fi is the Poynting vector. But because all quantities are expressed 
in terms of retarded time t' instead of i, it is convenient to convert this "per-unit- 
i-time" flux to a "per-unit-t'-time" measure. We have that 

is the energy crossing an element of S' per unit i'-time, and the conversation factor 
is dt/dt' = 1 + c^dR^/dt' = 1 - v(t') ■ R(t')/c = k. The displayed quantity can 
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therefore be expressed as S(t') ■ R(t')R 2 (t')Kdfl' , which is a definition for dP' = 
(dP'/dQ')dSl'. We have arrived at 

dP' 

^— = kR 2 S-R (6.4.3) 
dlV 

= energy crossing an element of S' 
per unit i'-time and unit solid angle. 

Because it is understood that all quantities refer to the retarded time t', we no 
longer have a need for the redundant notation [ ] rct . Integrating Eq. (6.4.3) over 
dtt' gives the total power crossing S' . Integrating this over dt' (which involves 
displacing S' in response to the particle's motion during integration) returns the 
total energy radiated by the charge. 
The Poynting vector is given by 

S = — E xB = —E x(Rx E) = —\E\ 2 R, 
Mo Moc Moc 

because E is orthogonal to R. This gives 

dP' 1 _ kR 2\ e \2_ 1 I 2 \Rx[(R-v/c)xa]\ 2 



(6.4.4) 



dQ' HqC [i c (47re ) 2 k 5 c 4 

or 

dP' _ n q 2 \Rx [(R-v/c) x a}\ 2 

dSV ~ (4tt) 2 c ^ ' 

The general expression for the angular profile is complicated, and it depends on the 
individual directions of v, a, and R. To gain insight into this it will help to look at 
specific simple situations; we shall do so in the following subsections. To ease the 
notation we shall now agree to drop the primes on dP' /dQ! . 

6.4-2 Slow motion: Larmor formula 

Our first special case will be the limit v <C c. In this slow-motion approximation 
we can neglect v/c in front of R in Eq. (6.4.4), and also set k ~ 1. We then have 

dP /j, q 2 ,~ - |2 

dh * (4^ |i?X(i?Xa)| 

~ ^ \a-{a-R)R\ 2 



(4tt) 2 c 



Mo<7 



2 



(4tt) 2 c 



a\ z -(a-RY 



If 9 denotes the angle between the vectors a and R, then \a\ 2 — (a ■ R) 2 = |a| 2 (l — 
cos 2 9) — \a\ 2 sin 2 9, and we finally obtain 

dP /i / h 2 



dn (47r) 2 c 



(q\a\) sin 2 9. (6.4.5) 



This expression could have been derived on the basis of the electric-dipole ap- 
proximation of Sec. 5.4. According to the results of that section, in particular 
Eq. (5.4.22), the angular profile of radiation produced by a slowly-moving source is 
given by 

dP _ Mo | ..,2 • 2 , 



dn (4tt) 



|p| 2 sin^6, 
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where p is the dipole moment of the charge distribution, and O the angle between 
the vectors r = x/\x\ and p. Here p = qr(t), so that p = qa. And because the 
sphere S' is centered on the particle, the vector R plays the same role as f, so that 
O = 9. The electric-dipole approximation therefore produces the same result as in 
Eq. (6.4.5). 

Integration of dP/dfl over all angles yields 




This result is known as Larmor's formula. 



6.4-3 Linear motion 

In this subsection we no longer restrict the size of v /c, but we assume that v and a 
are pointing in the same direction. This condition might hold at all times, but it is 
sufficient to assume that it holds only momentarily. At this moment the particle is 
moving on a straight line. We now let 9 be the angle between R and the common 
direction of v and a. 

We have that v x a = 0, so Eq. (6.4.4) becomes 

dP _ fi q 2 \R x (R x a)\ 2 
dll ~ (4tt) 2 c «5 ' 

We already know that \Rx (Rxa)\ 2 = \a\ 2 sin 2 9, and we recall that k = l—v-R/c = 
1 — vcos9/c. Substituting these results gives 

dP Mo / i ,n2 sin 2 (9 

dn = wrM a ^ (i-.ccsw ( } 

For v jc <C 1 this reproduces the sin 2 9 behaviour encountered in the preceding 
subsection. For v/c approaching unity, however, the angular profile is strongly 
modified — it tips forward in the direction of v (9 = 0), and the power radiated 
at small angles is very large (see Fig. 6.4). For such values of v/c, the radiation is 
beamed in the forward direction. 
The angular function is 

sin 2 ,„ . „n 

/= ( l-,cosW (6A8) 

Differentiating with respect to \i = cos 9 and putting the result to zero indicates 
that the maximum of the distribution is at the value of [i that solves (3v/c)/j, 2 + 
2fi — 5v/c = 0. So most of the radiation is emitted at angles close to # max , which is 
determined by 

co^= ^™f -. (6.4.9) 

For v/c <C 1, this expression is well approximated by cos # max — |(^/c) — ^§-(v/c) 3 + 
■ ■ ■, and # max ~ 7r/2: most of the radiation is emitted in the directions transverse 
to v. This is the usual nonrelativistic situation. For v/c close to unity, however, so 

that 7^1, Eq. (6.4.9) is well approximated by cos# max ~ 1 — §7~ 2 — ^7 -4 H , 

and 6* max is very small: most of the radiation is emitted in the forward direction. 
For such small angles we have cos6> max = 1 — |#max + ■ • •, and this gives us the 
approximate relation 

# max ^, 7»1- (6-4.10) 
27 
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Figure 6.4: The first plot shows the angular profile of radiation produced by a 
nonrelativistic particle in linear motion; the radiation propagates mostly in the 
directions transverse to the charge's velocity (the horizontal direction in the plot). 
The second plot shows the angular profile of radiation produced by a relativistic 
particle in linear motion; the radiation is beamed in the forward direction with an 
angular width of order 1/(27). Notice the different scales in the two diagrams. 
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This is the statement that when the charge moves with a highly relativistic speed, 
the radiation is strongly beamed in the direction of the velocity vector. 

The total power radiated is obtained by integrating Eq. (6.4.7) with respect to 
dQ = sin 9 dOdS. We have 



Mo 



(47T) 2 C 



(q\a\) 2 2n 



sm 3 0d0 



(1 - wcos0/c) 5 ' 



and the integral evaluates to §7 • The final answer is therefore 



(6.4.11) 



This generalizes the Larmor formula of Eq. (6.4.6). It is instructive to express this 
in terms of the force acting on the particle. This force is F = dp/dt, where p = mjv 
is the particle's relativistic momentum. For linear motion in the x direction we have 
p x = mv/(l - v 2 /c 2 ) 1 / 2 and 



dp x 
dt 



ma 



+ 



{mv){v/c 2 )a 



ma 



{1-V 2 / C 2 ) 1 / 2 (1 _ i; 2/ c 2)3/2 (! _ w 2/ c 2)3/2- 

This result can be expressed as F — m-f 3 a, and Eq. (6.4.11) becomes 

2 



Mo 

67TC 



m 



(6.4.12) 



This is the total power radiated by a charge in instantaneous linear motion. 



6.4-4 Circular motion 

We now take a and v to be orthogonal vectors. This condition might hold at all 
times, but it is sufficient to assume that it holds only momentarily. At this moment 
the particle's motion is circular. For concreteness we let a be directed along the x 
axis, and v points in the z direction. We let the angles {9, <j>) give the direction of 
the vector R relative to this set of axes. We therefore have a = ax, v = vz, and 

R = sin 9 cos </> x + sin 9 sin <j> y + cos 9 z. 

According to this, v ■ R = v cos 9, a ■ R = a sin 9 cos <f), n — 1 — v cos9/c, and 

R x [(R - v/c) x a] = (R ■ a)(R — v/c) — R ■ (R — v/c)a 



{Ft ■ a){R— v/c) — na. 



The squared norm of this vector is 
\Rx [{R-v/c) x a}\ 2 



{R ■ a) 2 {R - v/c) ■ {R-v/c) 
- 2k{R ■ a){R - v/c) ■ a + n 2 a 2 
{R ■ a) 2 {\ - 2R ■ v/c + v 2 /c 2 - 2k) + n 2 a 2 
-{R-a) 2 {l-v 2 /c 2 ) + K 2 a 2 
{R-a) 2 ~ 



2 2 

K a 



1 



K 2 7 2 



Substituting this result into Eq. (6.4.4) gives 



Mo q 2 a 2 



dP _ 
dn ~ (4tt) 2 c k 3 



(R-a) 2 

k 2 7 2 



1 - 
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Figure 6.5: Angular profile of radiation produced by a relativistic charge in circular 
motion. The radiation is beamed in the (horizontal) direction of the velocity vector. 



dP fi q \a\ sin: 



2 , 



7 2 (1 — vcos8/c) 2 



.4.13) 



dfl {Att) 2 c(1-v cos 9/cf 

This angular profile is rather complex, and quantitatively different from the one 
encountered in the preceding subsection; but it still features a strong beaming in 
the forward direction (direction of v) when v is close to the speed of light. This is 
illustrated in Fig. 6.5. 

The total power radiated is obtained by integrating dP/dD, over the angles 9 
and <j>. The result is 

P=^ c { q \a\) 2 1 \ (6.4.14) 

This is a different generalization of the Larmor formula of Eq. (6.4.6). We may also 
express this in terms of the force F acting on the particle. In this case the force 
changes only the direction of the vector v, but not its magnitude. So F = dp/dt = 
d(mjv)/dt = mja, and Eq. (6.4.14) becomes 



p= Mo fq\F\ 
6irc \ m 



Y- (6.4.15) 



This is the total power radiated by a charge in instantaneous circular motion. You 
should compare this with Eq. (6.4.12) and notice the extra factor of j 2 . 

6.4-5 Synchrotron radiation 

We continue our discussion of circular motion and provide additional details. We 
now assume that the charge moves at all times on a circle of radius p at a constant 
angular velocity f2. Its position vector is 

r(t') = p[cos(fit')* + sm(nt')y] , (6.4.16) 

its velocity vector is 

v[t') =v[- sm(fLt')x + cos(Qt')y] , (6.4.17) 
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where v = Vlp, and its acceleration vector is 

a(t') = -vQ[cos(Qt')x + sm(Qt')y] . (6.4.18) 

We wish to calculate the electric field that would be measured by a detector situated 
far away, in the wave zone. We place this detector on the y axis, at position x = ry, 
and we assume that r ^> p. 

The wave-zone electric field is given by Eq. (6.4.1), which we first evaluate as 
a function of retarded time t'. The relation between t' and true time t is given by 
Eq. (6.4.20). It is useful to recall that t' is time measured by a clock that moves 
with the charged particle, while t is time measured by a clock attached to the 
detector. As we shall see, these measurements of time differ substantially when v is 
comparable to c, and this leads to an interesting modulation of the electromagnetic 
wave. 

We use the condition r > p to calculate E approximately. We shall need 
some level of precision in the calculation of t(t'), because it is this function which 
determines the phasing of the electric field. A much cruder approximation will 
suffice when calculating the field's amplitude. 

Working to first order in p/r, it is easy to see that the length of the vector 
R = x r(t') is given by R{t') = r — psinfM'. Equation (6.3.21) then gives 
t = t' + r/c — (p/c) sin fit', or 

Q(t -r/c) = fit' - - sin Sit'. (6.4.19) 

c 

This gives the relationship between t — r/c, the retarded time measured at the 
detector, and t' , the time measured at the particle. From this it follows that 

k= ^* = 1 - - cos fit'. (6.4.20) 
dt' c J 

We see that k will deviate strongly from 1 when v is close to the speed of light and 
when | cosfM'| is close to unity. 

To perform the remaining calculations we allow ourselves to neglect all terms 
linear in p/r. We take R = y, 

R-v/c= {v/c)sm(flt')x- [l - (v/c) cos(ftt')]y, 

and Rx (R—v/c) x a = — vQ(v/c — cos Slt')x. Substituting all this into Eq. (6.4.1) 
gives E = Ex, with 

An e c 2 r [1 -(v/c) cos Qt'] 3 ' ( " ' 

The electric field can be expressed in terms of detector time t by inverting Eq. (6.4.19). 
Its behaviour is displayed in Fig. 6.6 for selected values of 7 = [1 — (u/c) 2 ] -1 / 2 . The 
plots reveal that for 7 moderately larger than unity, the wave deviates strongly from 
a simple sinusoid. The phase modulations which originate from the discrepancy be- 
tween t and t' produce a wave that is characterized by two distinct timescales: a 
long one associated with the rotational frequency ft, and a short one associated 
with the pulse duration. The frequency spectrum of the radiation is presented in 
Fig. 6.7. These plots show that in the relativistic regime, the radiation is carried 
by a large number of harmonics of the fundamental frequency Q; as 7 increases the 
dominant harmonic is shifted to higher frequencies. 

To understand these features it is helpful to consider the extreme relativistic 
regime, in which 7 ^> 1. In this regime the basic mechanism behind the narrowing 
of the pulses is very easy to identify: When v/c — ► 1 the denominator of Eq. (6.4.21) 
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synchrotron radiation: gamma = 1.8, via = 0.83 synchrotron radiation: gamma = 2.0, via = 0.87 



Figure 6.6: Electric field produced by a charged particle in circular motion, as 
a function of detector time t, for selected values of the relativistic 7 factor. The 
plots also show k = dt/dt', whose deviation from 1 measures the magnitude of 
the relativistic time-dilation effect. The relativistic modulation of the wave pulse 
for large values of v/c is clearly apparent. In particular, the plots reveal that the 
radiation is characterized by two distinct time scales: a long one associated with 
the rotational frequency f2, and a short one associated with the pulse duration. 
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Figure 6.7: Fourier transform of the electric field produced by a charged particle 
in circular motion, for selected values of the relativistic 7 factor. The plots reveal 
that the number of relevant harmonics of the fundamental frequency f2 increases 
with 7. In addition, the dominant harmonic is shifted to higher frequencies. 
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approaches zero whenever cos fit' — 1 , which happens whenever fit' is a multiple of 
2ir. Most of the radiation is therefore emitted during a very short portion of the 
particle's orbit, just as the particle crosses the x axis in the positive sense. This 
radiation is strongly beamed in the forward direction, and it travels along the y 
axis toward the detector. During the remaining portion of the orbit little radiation 
is emitted, and the detector must wait for the completion of a full orbit before 
receiving the next pulse. 

To express these ideas mathematically we write v/c= 1 — |7~ 2 + 0(-f~ 4 ) and we 
approximate sin fit' ~ fit' and cos fit' ~ 1 — i(fit') 2 . From Eq. (6.4.20) we obtain 



1 



i + (jfit'f 



which shows that k is small (of order 7 2 <C 1) for a duration of order ( 7 fi)~ 
terms of particle time t'. We also obtain 

v/c - cos fit' ~ - ^ 1 - (l^' Y 

and Eq. (6.4.21) reduces to 

E ~ 1 »" l7 4 l-(7^) 2 (6 4 22) 

tj ~ 47re oC 2r 47 [1 + ( 7 fit') 2 ] 3 ' ( ' 

Again, this shows that in terms of particle time t', the radiation is "on" for a 
duration of order (jfl)^ 1 . 

To see how long the pulse lasts in terms of detector time t, we approximate 
Eq. (6.4.19) as fl(t - r/c) ~ \^~ 2 flt' , which gives 

-/fit' ~ 2 7 3 fi(t - r/c). (6.4.23) 

This relation reveals that the pulse duration, as measured by a clock attached to 
the detector (and as displayed in Fig. 6.6) is given by 

AW ~ (6.4.24) 

It is reduced by a factor 7~ 3 <C 1 relative to the orbital timescale associated with 
the angular velocity fl. As a consequence, the characteristic frequency uj c of the 
radiation is shifted relative to the fundamental frequency fl. It is given by lo c ~ ^fl, 
and the displacement of the spectrum toward higher frequencies when 7 increases 
is seen very clearly in Fig. 6.7. 



6.4-6 Extremely relativistic motion 

To conclude this section we note that at any point on the arbitrary trajectory of a 
point charge, the applied force can be decomposed into a component that is parallel 
to the velocity vector v, and a component F± that is perpendicular to it (Fj_-v = 0). 
When the motion is extremely relativistic, 7 ^> 1 and the total power emitted by 
the particle is well approximated by Eq. (6.4.15), 

p= jM L (q\F±\\ 2 2 _ 

6-kc \ m J 

The contribution from the parallel force is smaller by a factor of I/7 2 <C 1, and 
it can be neglected. In this limit, therefore, the parallel force plays no role in the 
energetics of the radiation. 
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6.5 Radiation reaction 

[The material presented in this section is also covered in Sees. 16.1, 16.2, and 16.3 
of Jackson's text.] 

In Sec. 5.6 we saw that a classical electron orbiting a classical proton emits 
electromagnetic radiation. As a result, the classical hydrogen atom loses energy at 
a rate determined by Larmor's formula, 

dE /J.Q , ,»2 

rfT = - p = -6^ |a| ) ' 

This loss of energy translates into a decrease in the electron's orbital radius. The 
orbit gradually shrinks, and the classical atom is unstable to the emission of elec- 
tromagnetic radiation. 

It is natural to assume that the electron is subjected only to the electric held 
produced by the proton. This picture, however, is incomplete: So long as the 
electron feels only the proton's Coulomb held, its motion must stay circular (or more 
generally, elliptical). The electron's inspiraling motion cannot take place in this 
picture, in contradiction with our previous description. We are forced to conclude 
that the proton's electric held cannot be the only held acting on the electron. The 
electron must also be subjected to its own electric field. It is the electron's self-held 
which provides the force that drives the electron's inspiraling motion. 

This conclusion is difficult to understand, because the electron's own held di- 
verges at the electron's position. How can it then produce a (presumably hnitc) 
radiation-reaction force that drives the inspiral? The answer is that the divergent 
part of the electron's self-held exerts no force — it just travels along with the elec- 
tron. The part of the self-held that is responsible for the radiation reaction is hnitc. 

To hgure out what this part is, let us ask what would happen if instead of 
imposing outgoing-wave boundary conditions on our spherical-wave solutions to 
Maxwell's equations, we were to impose incoming-wave boundary conditions. In 
this case we would have incoming waves bringing energy to the system instead of 
outgoing waves taking energy away. The atom's energy would then increase, and 
the electron would move out instead of in. We would, in fact, have a mirror image 
(more precisely, a time reversal) of the correct picture, and the radiation-reaction 
force would be acting in the opposite direction (pushing the orbit out instead of in) . 

We can hgure out which part of the electron's field is producing the radiation- 
reaction force by examining what happens when we switch from outgoing-wave 
boundary conditions to incoming-wave boundary conditions. We will see that while 
part of the held is not affected by this switch, the remaining part changes sign. 
And we will see that while the unaffected part is inhnite at the electron's position, 
the remaining part is hnite. It is this part of the field that is responsible for the 
radiation reaction. 

Let us examine the vector potential hrst. The physical solution to Maxwell's 
equations is the retarded solution, 

where j is the electron's current density. This describes waves that are properly 
outgoing in the wave zone. These waves remove energy from the system, and the 
radiation reaction drives the electron inward. An unphysical solution to Maxwell's 
equations is the advanced solution, 
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and this describes waves that are incoming in the wave zone. These waves bring 
energy to the system, and the radiation reaction drives the electron outward. To 
isolate the part of the (retarded) vector potential that is responsible for the radia- 
tion reaction, we must examine both solutions simultaneously, and determine what 
changes when we choose one sign over the other. 

We want to evaluate A± near the electron, in the near zone. We can therefore 
take \x — x'\ to be small, and Taylor-expand the current density: 



j(t±\x-x'\/c)=j(t)± 



X 



x'\ dj 



dt 



(*) + • 



Relative to the leading term, the second term is of order (r c /t c )/c = v c /c, and we 
shall take this to be small — we are working under a slow-motion approximation. 
We have 

Ha f j(t,x') 
A± = ^ 



We see that the first term is the same irrespective of the choice of boundary con- 
ditions (outgoing waves versus incoming waves); it cannot, therefore, contribute to 
a radiation-reaction force. The second term, however, changes sign under a change 
of boundary conditions, and the force it contributes will also change sign; this must 
be the vector potential's contribution to the reaction-reaction force. Because the 
retarded solution provides the correct sign, we shall write 



(6.5.1) 



and we will see that the fields this potential generates produce the radiation-reaction 
force that acts on the electron. 

For the scalar potential we have the two solutions 



$±(t,a:) 



1 



47T£o 



p(t± \x-x'\/c,xj_ 



IX — X' 



and we Taylor-expand the charge density: 
p(t ± \x - x'\/c) 



\ x -x'\ 3 d 3 P ^ + 



\x — X 

2c 2 



(Pp 
dt 2 



(*) 



± 



6c 3 dt 3 



This gives 



47re 



p(t,x') Sx , ± ld 



IX — X 



+ 



± 



2^2 at 2 " 
J_ <p_ 

6c 3 dt 3 



c dt 

J p{t,x')\x-x'\d 3 x' 
j p{t,x')\x-x'\ 2 d 3 x' 



p{t 1 x')d 3 x' 



The first and third terms are insensitive to boundary conditions; they have nothing 
to do with radiation reaction. The second and fourth terms, on the other hand, 
change sign under a switch of boundary conditions, and they are associated with 
radiation reaction. The second term, however, vanishes by virtue of charge conser- 
vation, and we conclude that the radiation-reaction part of the scalar potential is 
given by (recall that c~ 2 = eoMo) 



<f> IT (t,x) 



J^ld 3 ^ 
'4nc6dt 3 



J p(t,x')\x 



x 



' d 3 x' + ■■■. 



(6.5.2) 
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Up to now the charge and current densities have been left unspecified, and our 
considerations have been general: Eqs. (6.5.1) and (6.5.2) give the radiation-reaction 
potentials for an arbitrary distribution of change and current. These expressions 
will apply, in particular, to a point charge q if we set p(t,x') = qS(x' — r(t)) and 
j(t,x') = qv(t)5(x' — r(t)), where v = dr/dt. Substituting this into our expressions 
for the potentials returns 



$„{t,x) = - 



Mo q 



Q3 



and 



4vrc 6 dt 3 
Mo 



x — r(t)\ 



d , . 
Ancdt v ; 



A TT (t : X) 

2 = (x — r) ■ (x — r), and differentiating once with respect to t gives 



We have \x — r 

—2(x — r)-v. Differentiating once more gives — 2(x — r)-a+2v-v, and differentiating 
a third time yields — 2(x — r) ■ a + 2v ■ a + Av ■ a = —2(x — r) ■ d + 6t> • a. We finally 
arrive at 



and 



Mo 1 

$ rr (t,X = q 

Aire 



A rr (t, x) 



(x — r) ■ a v ■ a 



Mo 
Attc 



qa. 



(6.5.3) 



(6.5.4) 



Here we have defined a = dv/dt and a = da/dt. These are the radiation- reaction 
potentials acting on a point charge q moving on a trajectory r(t). 
The radiation-reaction part of the electric field is given by 



E r 



dt 



Mo 
Aire 



1 



q\ a 



or 



j— i Mo . 

E TT = - — qa. 
birc 



(6.5.5) 
(6.5.6) 
(6.5.7) 



Because A„ does not actually depend on x, we also have 

B„ = 0. 

The radiation-reaction force is therefore F rr = qE II} or 

F rr = ^q 2 a. 
birc 

This force acts on the point charge, and it originates from the charge's own electric 
field! 

It is easy to verify that the radiation-reaction force does negative work on the 
particle and takes energy away from it. As the charge moves with velocity v the 
force does work at a rate 

W = F lr ■ v = ^q 2 a ■ v. 
one 



We write 



a ■ v 



dt 



(a-v) 



a 



and we average W over a time interval T. This gives 



(W) 



o 



Mo 2 

67TC 

birc 



Wdt 
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T 
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T d , , , 1 
— (a • v) at — — 
dV ' T 



\a\ 2 dt 
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Supposing now that the motion is periodic with period T (as it would be in the case 
of atomic motion), or else that it begins and ends with a vanishing acceleration, we 
may discard the boundary terms and obtain 

{ W) = -^(\af) = -(P), (6.5.8) 

where P is the power radiated by the point charge, as given by the Larmor formula 
of Eq. (6.4.6). The work done by the radiation-reaction force therefore matches the 
energy taken away by electromagnetic radiation, and we have energy conservation. 
(It is noteworthy that this statement holds only on the average.) 

We have a complete resolution of the paradox stated at the beginning of this 
section. The electron's own field, or at least the finite part that is sensitive to the 
choice of boundary conditions, acts on the electron and prevents it from keeping its 
circular motion. (Indeed, it is easy to check that for a circular orbit, the radiation- 
reaction force of Eq. (6.5.7) points in the direction of — v and drives the orbit 
to a smaller radius.) Furthermore, the work done by this radiation-reaction force 
matches the amount of energy taken away by the radiation, as it well should if 
energy is to be conserved in the process. 



6.6 Problems 

1. A point charge q moves on a circle of fixed radius a with an angular frequency 
uj(t) that changes slowly with time. (The time scale over which lu changes is 
much longer than the rotational period.) We take aw to be smaller than c, 
but we do not assume that it is much smaller than c; the charge has a Lorentz 
7-factor that can be large compared with unity. 

The charge emits electromagnetic radiation. It therefore loses energy and 
this causes the angular frequency u to change with time. You are asked to 
derive an equation for dco/dt and to express your result in terms of u, a, and 
other relevant constants. (You may need to know that the charge's relativistic 
energy is E = jmc 2 , where m is the charge's rest mass.) 

2. (This problem is adapted from Jackson's problem 14.12.) A charge q moves 
in simple harmonic motion along the z axis. Its position is given as a function 
of t' by z(t!) = z cos(uit'), where z$ and w are constants. It is assumed that 
(3 = zquj/c is fairly close to unity, so that the motion is relativistic. 

a) Show that the time-averaged angular profile of the radiation is described 

by 

dP\ flo 2 2 1 + (/3/2) 2 cos 2 . 2a 
{qzoio ) T -55 , sin e, 



<Kl/ 327r 2 c vy u ; (1 - /3 2 cos 2 6>)V2 

where 9 is the angle from the z axis. The time average is defined by 
((• • •)) = T- 1 j T (- ■ •) dt', where T = 2tt/lo is the period. 

b) Produce a parameteric plot of the angular function for (3 — 0.95. Comment 
on its appearance in view of the beaming effect discussed in Sec. 6.3.3. 

3. A charge q moves on a trajectory described by z(t') = ^Jb 2 + (ct') 2 , where b 
is a constant. The particle comes in from infinity, it turns around at z = b 
when t' — 0, and it eventually returns to infinity. 

Calculate P, the total power emitted by the particle. Show in particular that 
this does not depend on time t' . 
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4. (This problem is adapted from Jackson's problem 14.7.) A nonrelativistic 
particle of charge q, mass m, and initial speed v is incident on a fixed charge 
Q with an impact parameter b that is large enough to ensure that the particle's 
deflection in the course of the collision is very small. 

Calculate the total energy radiated by the particle. Because the deflection 
angle is so small, it is sufficient to assume that the particle's trajectory is a 
straight line and that its velocity vector is constant. But this approximation 
must be introduced only after you have calculated the particle's acceleration. 



